Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,804,987 T→C 5.8% pseudogene (1470/1476 nt) arpB → pseudogene, ankyrin repeats

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,804,9870TC 5.8% 237.9 / 3.5 69pseudogene (1470/1476 nt)arpBpseudogene, ankyrin repeats
Reads supporting (aligned to +/- strand):  ref base T (30/35);  new base C (2/2);  total (32/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

GAACCCCTGCATTATACATCGCCATGAGCAAGGGTAATGAAGACGTCGTGTTATCTTATATATCGACGCTGGGTGCTTTTGCAAAAAAACATTCTTTTAGTCAACATCAGTTATTTACACTATTGGCCGCTAAAAATCATGACAACATGTCAGCTGTTCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGTAATAACAAAAAATCCCTTGATGCCTGTCCCTTTTGTTACACTCCGTTATCACGCACAAGAGATATGCAGGACACTGGTATGCCGACTAAACGCTTTGATAAAAAACACTGGAAGATGGTGGTGGTGCTACTGGCAATCTGTGGCGCTATGTTGTTGCTACGTTGGGCAGCAA  >  NC_000913/1804722‑1805258
                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                               
gAACCCCTGCATTATACATCGCCATGGGCAAGGGTAATGAAGACGTCGTGTTATCTTATATATCGACGCTGGGTGCTTTTGCAAAAAAACATTCTTTTAGTCAACATCAGTTATTTACACTATTGGCCGCTAAAAATCATGACAACATGTCAGCTGTTCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTATAAc                                                                                                                                                                                                                                                                          <  1:310594/273‑1 (MQ=255)
aaaCCCCTTCTTTTTACATCNNCATGANTAAGGGTAATGAAGACGTCTTTTTATCTTATATATCGACGCTGGGTGCTTTTGCAAAAAAACATTCTTTTAGTCAACATCAGTTATTTACACTATTGGCCGCTAAAAATCATGACAACATGTCAGCTGTTCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTATAAc                                                                                                                                                                                                                                                                          <  3:32269/272‑1 (MQ=255)
         cATTATACATCGCCATGAGCAAGGGTAATGAAGACGTCGTGTTATCTTATATATCGACGCTGGGTGCTTTTGCAAAAAAACATTCTTTTAGTCAACATCAGTTATTTACACTATTGGCCGCTAAAAATCATGACAACATGTCAGCTGTTCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTATAACCAGCTAttt                                                                                                                                                                                                                                                                 <  1:121611/273‑1 (MQ=255)
         cATTATACATAGCCATGAGCAAGGGTAATGAAGACGTTGTGTTATCTTATATATCGACGCTTGGTGCTTTTGCAAAAAAACATTCTTTTAGTCAACATCAGTTATTTACACTATTGGCCGCTAAAAATCATGACAACATGTCAGCTGTTCATATATCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTATAACCAGCTAttt                                                                                                                                                                                                                                                                 <  2:330476/273‑1 (MQ=255)
                                          aCGTCGTGTTATCTTATATATCGACGCTGGGTGCTTTTGCAAAAAAACATTCTTTTAGTCAACATCAGTTATTTACACTATTGGCCGCTAAAAATCATGACAACATGTCAGCTGTTCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCa                                                                                                                                                                                                                              >  2:76793/1‑275 (MQ=255)
                                                          tatatCGACGCTGGGTGCTTTTGCAAAAAAACATTCTTTTAGTCAACATCAGTTATTTACACTATTGGCCGCTAAAAATCATGACAACATGTCAGCTGTTCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTATAACCAGCTa                                                                                                                                                                                                                                                                    >  2:695502/1‑221 (MQ=255)
                                                              tCGACGCTGGGTGCTTTTGCAAAAAAACATTCTTTTAGTCAACATCAGTTATTTACACTATTGGCCGCTAAAAATCATGACAACATGTCAGCTGTTCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTATAACCAGCTATTTa                                                                                                                                                                                                                                                                <  4:695499/221‑1 (MQ=255)
                                                               cGACGCTGGGTGCTTTTGCAAAAAAACATTCTTTTAGTCAACATCAGTTATTTACACTATTGGCCGCTAAAAATCATGACAACATGTCAGCTGTTCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTATaa                                                                                                                                                                                                                                                                           >  2:736124/1‑209 (MQ=255)
                                                               cGACGCTGGGTGCTTTTGCAAAAAAACATTCTTTTAGTCAACATCAGTTATTTACACTATTGGCCGCTAAAAATCATGACAACATGTCAGCTGTTCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTATAACCAGCTATTTAc                                                                                                                                                                                                                                                               >  1:322623/1‑221 (MQ=255)
                                                               cGACGCTGGGTGCTTTTGCAAAAAAACATTCTTTTAGTCAACATCAGTTATTTACACTATTGGCCGCTAAAAATCATGACAACATGTCAGCTGTTCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTATAACCAGCTATTTAc                                                                                                                                                                                                                                                               >  2:200207/1‑221 (MQ=255)
                                                                   gCTGGGTGCTTTTGCAAAAAAACATTCTTTTAGTCAACATCAGTTATTTACACTATTGGCCGCTAAAAATCATGACAACATGTCAGCTGTTCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTATAACCAg                                                                                                                                                                                                                                                                       <  4:736121/209‑1 (MQ=255)
                                                                   gCTGGGTGCTTTTGCAAAAAAACATTCTTTTAGTCAACATCAGTTATTTACACTATTGGCCGCTAAAAATCATGACAACATGTCAGCTGTTCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTATAACCAGCTATTTACGata                                                                                                                                                                                                                                                           <  4:200204/221‑1 (MQ=255)
                                                                   gCTGGGTGCTTTTGCAAAAAAACATTCTTTTAGTCAACATCAGTTATTTACACTATTGGCCGCTAAAAATCATGACAACATGTCAGCTGTTCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTATAACCAGCTATTTACGata                                                                                                                                                                                                                                                           <  3:322623/221‑1 (MQ=255)
                                                                           cattttCAAAAAAACATTCTTTGAGTCAACATCAGTTATTTACACTNTTGGCCGCTAAAAAGCATGACAACATGTCAGCTGTTCATATAGCCATTAATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTATAACCAGCTAttt                                                                                                                                                                                                                                                                 <  4:993276/205‑1 (MQ=255)
                                                                                                                      acTATTGGCCGCTAAAAATCATGACAACATGTCAGCTGTTCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTAtt                                                                                                                                                                                                                                                                                >  7:76307/1‑149 (MQ=255)
                                                                                                                           tGGCCGCTAAAAATCATGACAACATGTCAGCTGTTCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTATaa                                                                                                                                                                                                                                                                           <  7:33071/149‑1 (MQ=255)
                                                                                                                               cGCTAAAAATCATGACAACATGTCAGCTGTTCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTg                                                                                                                                                                                                            <  1:125573/208‑1 (MQ=255)
                                                                                                                                  tAAAAATCATGACAACATGTCAGCTGTTCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTta                                                                                                                                                                                                                                                                              >  2:916394/1‑139 (MQ=255)
                                                                                                                                  tAAAAATCATGACAACATGTCAGCTGTTCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGTAATAACAAAAAATCCCTTGATGCCTGTCCCTTTTGTTaca                                                                                                                                      <  4:76791/275‑1 (MQ=255)
                                                                                                                                      aaTCATGACAACATGTCAGCTGTTCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTATAAc                                                                                                                                                                                                                                                                          <  4:916391/139‑1 (MQ=255)
                                                                                                                                        tCATGACAACATGTCAGCTGTTCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTATAACCAGCTATTTACg                                                                                                                                                                                                                                                              >  6:98705/1‑149 (MQ=255)
                                                                                                                                                         cTGTTCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTATAACCAGCTa                                                                                                                                                                                                                                                                    >  2:358852/1‑126 (MQ=255)
                                                                                                                                                            ttCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTATAACCAGCTATTTACGATATATATCGATACATAatt                                                                                                                                                                                                                                          <  8:76307/149‑1 (MQ=255)
                                                                                                                                                             tCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTATAACCAGCTATTTa                                                                                                                                                                                                                                                                <  4:358849/126‑1 (MQ=255)
                                                                                                                                                                      cATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATaaaa                                                                                                                                                                                                                         >  1:116108/1‑156 (MQ=255)
                                                                                                                                                                         tcatcaTAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTACTta                                                                                                                                                                                                                                                                              >  4:378806/1‑100 (MQ=255)
                                                                                                                                                                          catcatAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCa                                                                                                                                                                                                                     <  3:116108/156‑1 (MQ=255)
                                                                                                                                                                            tcatAAACATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGACAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGTAATAACAAAAAATCCCTTGATGCCTGTCCCTTTTGTTACACTCCGTTATCACGCACAAGAGATATGCAGGACACTGGTATGc                                                                                            <  2:257716/275‑1 (MQ=255)
                                                                                                                                                                             catAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTACTTATAAc                                                                                                                                                                                                                                                                          <  2:378809/100‑1 (MQ=255)
                                                                                                                                                                                 aGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCa                                                                                                                                                                                                                     >  7:53179/1‑149 (MQ=255)
                                                                                                                                                                                     ttATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGTAATAACAAAAAATCCCTTGATGCCTGTCCCTTTTGTTACACTCCGTTATCACGCACAAGAGATATGCAGGACACTGGTATGCCGACTaaa                                                                                    >  2:198831/1‑274 (MQ=255)
                                                                                                                                                                                      tataAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGTAATAACAAAAAATCCCTTGATGCCTGTCCCTTTTGTTACACTCCGTTATCACGCACAAGAGCTATGCAGGACACTGTTATGCCGACTAAACg                                                                                  >  2:262161/1‑275 (MQ=255)
                                                                                                                                                                                       ataAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGt                                                                                                                                                                                                             >  1:202495/1‑151 (MQ=255)
                                                                                                                                                                                           aaCTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGAcc                                                                                                                                                                                                         <  3:202495/151‑1 (MQ=255)
                                                                                                                                                                                             cTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGTAATAACAAAAAATCCCTTGATGCCTGTCCCTTTTGTTACACTCCGTTATCACGCACAATAGATATGCAGGACACTGGTATGCCGACTAAACGCTTTGAt                                                                           >  2:241760/1‑275 (MQ=255)
                                                                                                                                                                                              tGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATaa                                                                                                                                                                                                                           >  1:50579/1‑130 (MQ=255)
                                                                                                                                                                                                tAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTATAACCAGCTa                                                                                                                                                                                                                                                                    >  4:286499/1‑87 (MQ=255)
                                                                                                                                                                                                tAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGTAATAACAAAAAATCCCTTGATGCCTGTCCCTTTTGTTACACTCCGTTAACACGCACAAGAGATATGCAGGACACTGGTATGCCGACCAAACCCTTTTATaaa                                                                        >  4:798878/1‑275 (MQ=255)
                                                                                                                                                                                                  gAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACg                                                                                                                                                                                                                       <  3:50579/130‑1 (MQ=255)
                                                                                                                                                                                                   aaaCATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTATAACCAGCTATTTACGATATATATCGATACATaa                                                                                                                                                                                                                                            <  8:12492/108‑1 (MQ=255)
                                                                                                                                                                                                   aaaCATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTATAACCAGCTATTTACGATATATATCGATACATaa                                                                                                                                                                                                                                            >  7:12492/1‑108 (MQ=255)
                                                                                                                                                                                                    aaCATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTATAACCAGCTATTTa                                                                                                                                                                                                                                                                <  2:286502/87‑1 (MQ=255)
                                                                                                                                                                                                        tattatGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGTAATAACAAAAAATCCCTTGATGCCTGTCCCTTTTGTTACACTCCGTTAt                                                                                                                             >  2:443584/1‑214 (MQ=255)
                                                                                                                                                                                                            atGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGTAATAACAAAAAATCCCTTGATGCCTGTCCCTTTTGTTACACTCCGTTATCACg                                                                                                                         <  4:443581/214‑1 (MQ=255)
                                                                                                                                                                                                              gctgctATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCAcc                                                                                                                                                                                        <  5:27819/149‑1 (MQ=255)
                                                                                                                                                                                                                tgctATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCAc                                                                                                                                                                                         >  3:318070/1‑146 (MQ=255)
                                                                                                                                                                                                                 gctATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGa                                                                                                                                                                                     >  6:2133/1‑149 (MQ=255)
                                                                                                                                                                                                                  ctATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTATaa                                                                                                                                                                                                                                                                           >  3:45375/1‑62 (MQ=255)
                                                                                                                                                                                                                  ctATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGTAATAACAAAAAATCCCTTGATGCCTGTCCCTTTTGTTACACTCCGTTATCACGCACAAGAGATATGCAGGACACTGGTATGCCGACTAAACGCTTTGATAAAAAACACTGGAAGAtggtg                                                      >  2:512820/1‑275 (MQ=255)
                                                                                                                                                                                                                  ctATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGTAATAACAAAAAATCCCTTGATGCCTGTCCCTTTTGTTACAATCCGTTATCACGCACAAGAGATATGCAGGACACTGGTATGCCGACTAAACGCTTTGATAAAAAACACTGGAAGAtggt                                                       >  2:217144/1‑274 (MQ=255)
                                                                                                                                                                                                                    aTTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGa                                                                                                                                                                                     <  1:318070/146‑1 (MQ=255)
                                                                                                                                                                                                                     ttAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGTAATAACAAAAAATCCCTTGATGCCTGTCCCTTTTGTTACAATCCGTTATCACGCACAAGAGATATGCAGGACACTGGTATGCCGACTAAACGCTTTGATAAAAAACACTGGAAGAtggtggtg                                                   <  4:217141/275‑1 (MQ=255)
                                                                                                                                                                                                                      tAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTAt                                                                                                                                                                                                                                                                                 <  4:803157/52‑1 (MQ=255)
                                                                                                                                                                                                                      tAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTATAACCAg                                                                                                                                                                                                                                                                       <  1:45375/62‑1 (MQ=255)
                                                                                                                                                                                                                       aaTGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTATAACCAGCTATTTACGatat                                                                                                                                                                                                                                                          >  4:155338/1‑74 (MQ=255)
                                                                                                                                                                                                                           cAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTATAACCAGCTATTTACGatatatat                                                                                                                                                                                                                                                      <  2:155341/74‑1 (MQ=255)
                                                                                                                                                                                                                            aaTCAGCCAACGCCTGAGTCTTAGTGCTGATGAAATAAAGACGTATTTATAACCAGCTATCTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGTAATAACAAAAAATc                                                                                                                                                                >  1:78027/1‑159 (MQ=255)
                                                                                                                                                                                                                             aTCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGTAATAACAAAAAATCCCTTGATGCCTGTCCCTTTTGt                                                                                                                                          >  4:542590/1‑180 (MQ=255)
                                                                                                                                                                                                                              tCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGTAATAACAAAAAATCCCTTGATGCCTg                                                                                                                                                    >  1:334335/1‑169 (MQ=255)
                                                                                                                                                                                                                                atccAAAGCCTGAGTTTTAGNGCTTATGANCTAAAGACGTATTTATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTGTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGTAATAACAAAAAATCCCtt                                                                                                                                                            <  3:78027/157‑1 (MQ=255)
                                                                                                                                                                                                                                 gCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGTAATAACAAAAAATCCCTTGATGCCTGTCCCTTTTGTTaca                                                                                                                                      <  2:542593/180‑1 (MQ=255)
                                                                                                                                                                                                                                  ccAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGTAATAACAAAAAATCCCTTGATGCCTGTccc                                                                                                                                                <  3:334335/169‑1 (MQ=255)
                                                                                                                                                                                                                                    aaaGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTAGAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGTAATAACaaaaaa                                                                                                                                                                  <  7:58289/149‑1 (MQ=255)
                                                                                                                                                                                                                                      aGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTACTTATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGTAATAACaaa                                                                                                                                                                     >  2:134508/1‑144 (MQ=255)
                                                                                                                                                                                                                                          tGAGTTTTAGTGCTGATGAAATAAAGACGTACTTATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGTAATAACAAAAAAt                                                                                                                                                                 <  4:134505/144‑1 (MQ=255)
                                                                                                                                                                                                                                           gAGTTTTAGTGCTGATGAAATAAAGACGTATTTATAACCAGCTATTTACGATATATATCCATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGTAATAACAAAAAATCCCTTg                                                                                                                                                           <  8:46110/149‑1 (MQ=255)
                                                                                                                                                                                                                                                 tAGTGCTGATGAAATAAAGACGTATTTATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGTAATAACAAAAAATCCCTTGATGCCt                                                                                                                                                     <  5:2133/149‑1 (MQ=255)
                                                                                                                                                                                                                                                 tAGTGCTGATGAAATAAAGACGTATTTATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGTAATAACAAAAAATCCCTTGATGCCTGTCCCTTTTGTTACACTCCGTTATCACGCACAAGAGATATGCAGGACACTGGTATGCCGACTAAACGCTTTGATAAAAAACACTGGAAGATGGTGGTGGTGCTACTGGCAATCTGTGGCGCTAtgt                       <  4:198828/275‑1 (MQ=255)
                                                                                                                                                                                                                                                      cTGATGAAATAAAGACGTATTTATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGTAATAACAAAAAATCCCTTGATGCCTGTCCCTTTTGTTACACTCCGtt                                                                                                                               >  2:890294/1‑166 (MQ=255)
                                                                                                                                                                                                                                                        gatgaAATAAAGACGTATTTATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGTAATAACAAAAAATCCCTTGATGCCTGTCCCtt                                                                                                                                              >  7:95170/1‑149 (MQ=255)
                                                                                                                                                                                                                                                          tgaAATAAAGACGTATTTATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGTAATAACAAAAAATCCCTTGATGCCTGTCCCTTTTGTTACACTCCGTTATCa                                                                                                                           <  4:890291/166‑1 (MQ=255)
                                                                                                                                                                                                                                                                      ggttttAAAACCAGCAATTTACGAGATATATCGATACATAATTTAATTTTTTTAATAAAACGCAGAAACAGGGACCATGTGCACAACGCCACACGACGGCGGTAATAACAAAAAATCCCTTGATGCCTGTCCCTTTTGTTACACTCCGTTATCACGCACAAGAGATATGCAGGACACTGGTATGCCGACTAAACGCTTTGATAAAAAACACTGGAAGATGGTGGTGGTGCTACTGGCAATCTGTGGCGCTATGTTGTTGCTACGTTGGGCAGCaa  <  4:262158/272‑1 (MQ=255)
                                                                                                                                                                                                                                                                         tttATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGTAATAACAAAAAATCCCTTGATGCCTGTCCCTTTTGTTACACTCCGTTATCACGCACAAgaga                                                                                                                >  1:92880/1‑162 (MQ=255)
                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                               
GAACCCCTGCATTATACATCGCCATGAGCAAGGGTAATGAAGACGTCGTGTTATCTTATATATCGACGCTGGGTGCTTTTGCAAAAAAACATTCTTTTAGTCAACATCAGTTATTTACACTATTGGCCGCTAAAAATCATGACAACATGTCAGCTGTTCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGTAATAACAAAAAATCCCTTGATGCCTGTCCCTTTTGTTACACTCCGTTATCACGCACAAGAGATATGCAGGACACTGGTATGCCGACTAAACGCTTTGATAAAAAACACTGGAAGATGGTGGTGGTGCTACTGGCAATCTGTGGCGCTATGTTGTTGCTACGTTGGGCAGCAA  >  NC_000913/1804722‑1805258

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: