Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,484,461 A→C 7.8% A497A (GCT→GCG pepA ← multifunctional aminopeptidase A: a cyteinylglycinase, transcription regulator and site‑specific recombination factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,484,4610AC 7.8% 208.7 / 5.1 90A497A (GCT→GCGpepAmultifunctional aminopeptidase A: a cyteinylglycinase, transcription regulator and site‑specific recombination factor
Reads supporting (aligned to +/- strand):  ref base A (36/45);  new base C (6/2);  total (42/48)
Fisher's exact test for biased strand distribution p-value = 1.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.14e-01

CAGGTATACACCGCAGCGCCTGATGCGACGCTGTTGTGTCTTATCAGGCCTACGAGTTCAGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGTACTTACGGGTAAAGCGTGACAGGAAGCAACCTGCGG  >  NC_000913/4484316‑4484606
                                                                                                                                                 |                                                                                                                                                 
ctggTATACACCGCAGCGCCTGATGCGACGCTGTTGTGTCTTATCAGGCCTACGAGTTCAGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAgcg                                                                                                                                                <  2:42556/147‑1 (MQ=255)
  ggTATACACCGCAGCGCCTGATGCGACGCTGTTGTGTCTTATCAGGCCTACGAGTTCAGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGcgcg                                                                                                                                              >  1:858383/1‑149 (MQ=255)
   gTATACACCGCAGCGCCTGATGCGACGCTGTTGTGTCTTATCAGGCCTACGAGTTCAGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCTCCGTTAAACCCAGCGCgg                                                                                                                                             >  2:464828/1‑149 (MQ=255)
   gTATACACCGCAGCGCCTGATGCGACGCTGTTGTGTCTTATCAGGCCTACGAGTTCAGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAACCCCAGCGCgg                                                                                                                                             >  2:630828/1‑149 (MQ=255)
     ataCACCTCAGCGCCTGATGCGACGCTCTTGTGTCTTAGCAGGCCTACGAGTTCAGTGCTGTGGAGGTCGGATTAGGCGTTCACTCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCATCGCGGtt                                                                                                                                           <  2:924837/149‑1 (MQ=255)
     ataCACCGCAGCGCCTGATGCGACGCTGTTGTGTCTTATCAGGCCTACGAGTTCAGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGtt                                                                                                                                           <  1:971670/149‑1 (MQ=255)
      taCACCGCAGCGCCTGATGCGACGCTGTTGTGTGTTATCAGGCCTACGAGTTCAGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGttt                                                                                                                                          <  2:745874/149‑1 (MQ=255)
        cacCGCAGCGCCTGATGCGACGCTGTTGTGTCTTATCAGGCCTACGAGTTCAGTGCTGTGTAGGTCGGATAAGGCGATCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTaa                                                                                                                                        >  2:866393/1‑149 (MQ=255)
        cacCGCAGCGCCTGATGCGACGCTGTTGTGTCTTATCAGGCCTACGAGTTCAGTGCTCAGAAGGTCGGATAAGGCGATCACGCCGCATCCGGCAATAACACCCATGCCTGACGCAATTACTCTTCTCCGTTAAACACAGCCCGGTTTaa                                                                                                                                        >  2:862574/1‑149 (MQ=37)
         cccGCTGCGCCTGATGCGACGCTGTTGTGTCTTATCAGGCCTACGAGTTCAGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAAc                                                                                                                                       <  2:242819/148‑1 (MQ=255)
         acCGCAGCGCCTGATGCGACGCTGTTGTGTCTTATCAGGCCTACGAGTTCAGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTAGCCGTTAAACCCCGCGCGGTTTAAc                                                                                                                                       >  2:1182249/1‑149 (MQ=255)
          ccGCAGCGCCTGATGCGACGCTGTTGTGTCTTATCAGGCCTACGAGTTCAGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCGCGCTGGGTTTAAcg                                                                                                                                      <  2:1084093/149‑2 (MQ=255)
          ccGCAGCGCCTGATGCGACGCTGTTGTGTCTTATCAGGCCTACGAGTTCAGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACa                                                                                                                                      <  2:131320/149‑1 (MQ=255)
           cGCAGCGCCTGATGCGACGCTGTTGTGTCTTATCAGGCCTACGAGTTCAGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCGCGCTGGGTTTAAcgg                                                                                                                                     <  1:406670/149‑3 (MQ=255)
           cGCAGCGCCTGATGCGACGCTGTTGTGTCTTATCAGGCCTACGAGTTCAGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAg                                                                                                                                     >  2:167628/1‑149 (MQ=255)
           cGCAGCGCCTGATGCGACGCTGTTGTGTCTTATCAGGCCTACGAGTTCAGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAg                                                                                                                                     <  1:1146852/149‑1 (MQ=255)
            gCAGCGCCTGATGCGAGGCTGTTGTGTCTTATCAGGCCTACGAGTTCAGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGACTTGCCTGACGCAATTACTCTTCGACGTTAAACCCAGCGCGGTTTAACAgg                                                                                                                                    <  1:706277/149‑1 (MQ=255)
            gCAGCGCCTGATGCGACGCTGTTGTGTCTTATCAGGCCTACGAGTTCAGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGATGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAgg                                                                                                                                    <  1:817955/149‑1 (MQ=255)
               gcgcCTGATGCGACGCTGTTGTGTCTTATCAGGCCTACGAGTTCAGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAAc                                                                                                                                 >  1:310815/1‑149 (MQ=255)
                cgcCTGATGCGACTCTGTTGTGTCTTATAAGGCCTACGAGTTGAGTGCTGTGTAGGTCGGATAATGCGTTCACGCCGCATCCGGCAATAACAGCATTGACAGACGCAATTACTCTTCGCCGTTACACCCAGCGCGGATTAATAGGAACt                                                                                                                                >  2:366027/1‑149 (MQ=255)
                  ccTGATGCGACGCTGTTGTGTCTTATCAGGCCTACGAGTTCAGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTACACCCAGCGCGGTTTAACAGGAACtgt                                                                                                                              <  1:407941/149‑1 (MQ=255)
                   cTGAGGCGACGCTGGTGTGTCTTATCAGGCCTACGAGTTCAGGGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACGCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACtgtg                                                                                                                             <  1:223352/149‑1 (MQ=255)
                      aTGCGACGCTGTTGTGTCTTATCAGGCCTACGAGTTCAGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCa                                                                                                                          >  2:128564/1‑149 (MQ=255)
                       tGCGACGCTGTTGTGTCTTATCAGGCCTACGAGTTCAGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAg                                                                                                                         <  2:1139014/149‑1 (MQ=255)
                        gCGACGCTGTTGTGTCTTATCAGGCCTACGAGTTCAGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGc                                                                                                                        >  2:970466/1‑149 (MQ=255)
                          gACGCTGTTGTGTCTTATCAGGCCTACGACTTCAGTTCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGGGAAACCCAGCGCGGTTTAACAGTATCTGTGCCAGCaa                                                                                                                      <  2:625115/149‑1 (MQ=38)
                               tgttgtGTCTTATCAGGCCTACGAGTTCAGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTa                                                                                                                 >  2:1037606/1‑149 (MQ=255)
                                 ttgtGGCTTATCAGGCCTACGAGTTCAGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATACCAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTAcc                                                                                                               <  1:209488/149‑1 (MQ=255)
                                      tCTTATCAGGCCTACGAGTTCAGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTTCTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGgacg                                                                                                          >  1:306335/1‑149 (MQ=255)
                                      tCTTATCAGGCCTACGAGTTCAGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGgacg                                                                                                          >  2:401006/1‑149 (MQ=255)
                                      tCTTATCAGGCCTACGAGTTCAGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGgacg                                                                                                          >  1:352572/1‑149 (MQ=255)
                                              ggCCTACGAGTTCAGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTgg                                                                                                  >  2:1153717/1‑149 (MQ=255)
                                                  tACGAGTTCAGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCATCAACGCTACCGGACGACCGGTGGCGcc                                                                                              <  2:722108/149‑1 (MQ=255)
                                                  tACGAGTTCAGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGcc                                                                                              >  1:595044/1‑149 (MQ=255)
                                                  tACGACTTCAGTGCTGTGTAGGTCGGAAAAGGCGTGCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACGCGTCGCCGTTCAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCCACCGGTCGTCCGGTAgcgtt                                                                                              <  1:215551/149‑3 (MQ=255)
                                                     gAGTTCAGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGCCGACCGGTGGCGCCttt                                                                                           >  1:53898/1‑149 (MQ=255)
                                                       gTTCAGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCGACGCTACCGGACGACCGGTGGCGCCTTTTg                                                                                         >  2:918814/1‑149 (MQ=255)
                                                        ttCAGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGc                                                                                        <  2:579064/149‑1 (MQ=255)
                                                        ttCAGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGc                                                                                        >  1:281506/1‑149 (MQ=255)
                                                        ttCAGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGc                                                                                        >  1:373462/1‑149 (MQ=255)
                                                        gtCAGTGCTGTGTAGGTCGGATAAGGCTTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGc                                                                                        <  1:399804/148‑1 (MQ=255)
                                                         tCAGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGGATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCt                                                                                       <  1:743955/149‑1 (MQ=255)
                                                          cAGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCtt                                                                                      >  2:851842/1‑149 (MQ=255)
                                                           aGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCttt                                                                                     <  1:848695/149‑1 (MQ=255)
                                                           aGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCACCAACGATACCGTACGCCCGGTGGCGCCATATGCttt                                                                                     >  2:801233/1‑149 (MQ=39)
                                                            gTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGAAATAACAGCCTCGCCTGACGCAATTACTCATCGCCGTTAAACCCTCCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTa                                                                                    <  2:1145697/149‑1 (MQ=255)
                                                             tGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGAGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTAc                                                                                   >  2:286320/1‑149 (MQ=255)
                                                             ggCTGTGGAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTAc                                                                                   <  2:858962/148‑1 (MQ=255)
                                                                       gTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCa                                                                         >  2:67018/1‑149 (MQ=255)
                                                                       gTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCa                                                                         >  1:1155532/1‑149 (MQ=255)
                                                                       gTCGGATAAGGCGTTCACGCCGCAACCGGCAATAACAGCCTTGCCTGAAGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTAGCTTTACCAGAACGCCa                                                                         >  1:131100/1‑149 (MQ=255)
                                                                           gATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCg                                                                     >  2:913109/1‑149 (MQ=255)
                                                                             tAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCCGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGt                                                                   >  2:730883/1‑149 (MQ=255)
                                                                                ggCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATAACTCTTCGCCGTTAAACCCCGCGCGGATTAACAGGAACTGTGCCAGAAACGATACCGGACGACCGGTGCCGCCTTTTGCTTAACCAGAACGCCAGGCCGAAcc                                                                >  1:541573/1‑149 (MQ=255)
                                                                                 gCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAACCCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCg                                                               >  1:1122731/1‑149 (MQ=255)
                                                                                   gTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCCGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGc                                                             >  1:969052/1‑149 (MQ=255)
                                                                                         gccgcATCCGGCAATAACCGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATc                                                       >  1:473408/1‑149 (MQ=255)
                                                                                          ccgcATCCGGCAATATCAGCCTTGCCTGTCTCAATTACTCTGCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACGAGAACGCCAGGCGGTACCGGCGATATcc                                                      <  2:142951/149‑1 (MQ=255)
                                                                                          ccgcATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGTCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATcc                                                      <  1:840452/149‑1 (MQ=255)
                                                                                           cgcATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCa                                                     <  2:699906/149‑1 (MQ=255)
                                                                                             cATCCGGCAATAACAGCGTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAgg                                                   <  1:512465/149‑1 (MQ=255)
                                                                                              aTCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTCAACCCAGCGCGGTTTAACAGGAACTGGGCCAGCAACGCTACCGTACGCCCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGt                                                  <  1:726175/149‑1 (MQ=255)
                                                                                              aTCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGt                                                  <  2:310815/149‑1 (MQ=255)
                                                                                                ccGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTgc                                                <  2:406008/149‑1 (MQ=255)
                                                                                                 cGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTAAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTgcg                                               <  1:971525/149‑1 (MQ=255)
                                                                                                  ggCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAACCCCCGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCCGGCGGTACCGGCGATAGCCAGGTgcgc                                              >  1:365594/1‑149 (MQ=255)
                                                                                                        aaCAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGACAGCAACGCTACCGGACGACCGGTGGCGCCTTGTGGTTTACCAGAACGCCAGGCGGTACCGGCGCTATCCCGGTGCGCCCAggg                                        >  1:697635/1‑147 (MQ=255)
                                                                                                         aCAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTg                                       <  1:866393/149‑1 (MQ=255)
                                                                                                         aCAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGATTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTg                                       <  1:862574/149‑1 (MQ=255)
                                                                                                          cAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGt                                      >  2:1173486/1‑149 (MQ=255)
                                                                                                             ccTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGTACt                                   <  2:514935/149‑1 (MQ=255)
                                                                                                             ccTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGTACt                                   <  1:970466/149‑1 (MQ=255)
                                                                                                                 gCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGTACTTACg                               >  1:785650/1‑149 (MQ=255)
                                                                                                                 gCCTGACGCAATTACTCTTCGACGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGTACTTACg                               <  2:569814/149‑1 (MQ=255)
                                                                                                                    tGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGTACTTACGGGt                            <  1:558887/149‑1 (MQ=255)
                                                                                                                      aCGCAATTACTCATCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGTACTTACGGGTaa                          <  1:699140/149‑1 (MQ=255)
                                                                                                                         cAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGTACTTACGGGTAAAGc                       <  1:67018/149‑1 (MQ=255)
                                                                                                                             tACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGTACTTACGGGTAAGGCGTGa                   >  1:1070073/1‑149 (MQ=255)
                                                                                                                               ctctTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGTACTTACGGGTAAAGCGTGACa                 <  2:1157592/149‑1 (MQ=255)
                                                                                                                                 ctTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGTACTTACGGGTAAAGCGTGACAgg               <  1:1153717/149‑1 (MQ=255)
                                                                                                                                   tCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGTACTTACGGGTAAAGCGTGACAGGaa             <  1:1037606/149‑1 (MQ=255)
                                                                                                                                   tCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGTACTTACGGGTAAAGCGTGACAGGaa             <  2:48134/149‑1 (MQ=255)
                                                                                                                                   tCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGTACTTACGGGTAAAGCGTGACAGGaa             <  1:401006/149‑1 (MQ=255)
                                                                                                                                    cGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGTACTTACGGGTAAAGCGTGACAGGAAg            <  2:373462/149‑1 (MQ=255)
                                                                                                                                      ctgTTAAACCGCGCTGGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGTACTTACGGGTAAAGCGTGACAGGAAGCa          >  2:391683/3‑149 (MQ=255)
                                                                                                                                        gTTAAACCGCGCTGGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGTACTTACGGGTAAAGCGTGACAGGAAGCAAc        >  1:308830/1‑149 (MQ=255)
                                                                                                                                         ttAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGTACTTACGGGTAAAGCGTGACAGGAAGCAAcc       >  2:627841/1‑149 (MQ=255)
                                                                                                                                            aaCCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGTACTTACGGGTAAAGCGTGACAGGAAGCAACCTGc    <  1:393160/149‑1 (MQ=255)
                                                                                                                                             aCCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGTACTTACGGGTAAAGCGTGACAGGAAGCAACCTGCg   >  1:386505/1‑149 (MQ=255)
                                                                                                                                              cccAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGTACTTACGGGTAAAGCGTGACAGGAAGCAACCTGCgg  <  1:402231/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
CAGGTATACACCGCAGCGCCTGATGCGACGCTGTTGTGTCTTATCAGGCCTACGAGTTCAGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGTACTTACGGGTAAAGCGTGACAGGAAGCAACCTGCGG  >  NC_000913/4484316‑4484606

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: