Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,466,138 G→A 7.1% pseudogene (747/2513 nt) ydbA → pseudogene, autotransporter homolog; interrupted by IS2 and IS30

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,466,1380GA 7.1% 324.4 / 7.2 131pseudogene (747/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base G (46/68);  new base A (4/5);  total (53/78)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 134 total aligned reads displayed.

GGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCAT  >  NC_000913/1465990‑1466287
                                                                                                                                                    |                                                                                                                                                     
gggTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGTGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCTTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGa                                                                                                                                                      <  1:1248470/150‑1 (MQ=255)
 ggTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaa                                                                                                                                                     >  2:1368701/1‑150 (MQ=255)
 ggTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaa                                                                                                                                                     >  1:879677/1‑150 (MQ=255)
   tGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAg                                                                                                                                                   >  1:1561841/1‑150 (MQ=255)
       cgcAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGAt                                                                                                                                               <  1:737305/150‑1 (MQ=255)
        gcAGATGTCTCAACGAAGTGGCAATAAAATCAAAACGGTGAGCTTGTTATTACGGGTTATAATGACACAGTAAACAACAATTGCAAAAACACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCCACGGTAATAACGGGAAAGTGAtt                                                                                                                                              <  1:348504/150‑1 (MQ=25)
          aGATGTCGCAAAGAAGTGGCAATATAATCACTACTGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAATGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCa                                                                                                                                            <  2:763750/150‑1 (MQ=21)
           gATGGCGCAACGAAGGCGCAATATAATCACAACGGTCATCTTGTTATTACGGGGGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAg                                                                                                                                           <  2:144178/150‑1 (MQ=25)
            attgcgcAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAgg                                                                                                                                          <  2:387779/147‑1 (MQ=255)
              ggCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGAc                                                                                                                                        <  1:1982660/150‑1 (MQ=255)
               gcgcAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACg                                                                                                                                       >  1:1659333/1‑150 (MQ=255)
                 gcAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGAAAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGc                                                                                                                                     <  2:130650/150‑1 (MQ=255)
                    accaaGTGCCAATATAATCACAACGGTGAGATAGTTATTACGGGTGATAATGCCACAGTAAACAAAAATGGAAAAACCACCGTTTACGGAAAAGATTCCACCGGTACTGAAATCAACTCTAATAACCGTAAAGTGATTCAGCACGGCGaa                                                                                                                                  >  2:639613/4‑149 (MQ=255)
                    aCGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGAt                                                                                                                                  >  2:1083836/1‑150 (MQ=255)
                      gAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAACGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCt                                                                                                                                >  2:1629128/1‑150 (MQ=255)
                        aGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTgg                                                                                                                              <  1:1404485/150‑1 (MQ=255)
                         gTGGCAATATAATCACAACGGTGAGCTTGGTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGa                                                                                                                             <  2:1380055/150‑1 (MQ=255)
                           ggCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGTTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATg                                                                                                                           <  1:1863493/150‑1 (MQ=255)
                           ggCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTGAGGACGGCGATCTGGATg                                                                                                                           <  2:371200/150‑1 (MQ=255)
                           ggCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAAAGGTAATAACGGGAAAGAGATTCAGGACGGCGATCTGGAtt                                                                                                                           >  2:515912/1‑149 (MQ=255)
                             cAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTTATTCAGGACGGCGATCTGGATGTc                                                                                                                         >  2:2020777/1‑150 (MQ=255)
                                  taATCACAACGGTGAGCTTGTTATTACGGGTGATTATGCCACAGTAAACAACAATGGCAAAACCACCGTGGACGGTAAAGATTCCACCGGTACGGGAATCAACGGTTATAACGGGATAGTGAGTCAGGACGGCGAGCTGGTGGTCAgcgg                                                                                                                    >  2:1913072/1‑150 (MQ=255)
                                  taATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAgcgg                                                                                                                    >  1:1339321/1‑150 (MQ=255)
                                  taATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAgcgg                                                                                                                    >  2:1242586/1‑150 (MQ=255)
                                   aaTCAAAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAACCCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAgcggc                                                                                                                   <  1:806536/150‑1 (MQ=255)
                                      cactACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAATATTCCACCGGTACGGAAATCACCGGTAATAACGGGAAAGTGATGCAGGCCGGCGATCTGGATGTCAgccgcggc                                                                                                                >  2:362436/1‑150 (MQ=21)
                                        caACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGcggcggcgg                                                                                                              <  1:574700/150‑1 (MQ=255)
                                        caACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGcggcggcgg                                                                                                              >  2:883218/1‑150 (MQ=255)
                                           cGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCTATGAAGCTAAAGATTCAGCAGGAAAGGAAAACACAGATAAAAACTGCAATGTGATGCAGGACAGAGATCTGGATGTGAGCGGCGGCGGTcc                                                                                                           >  1:100587/1‑149 (MQ=255)
                                            ggTGAGCTTGTTATTACGGGTGATAATGACACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAAAAACGGAAAAGTTATTCAGGACGGCGATCTGGATTTCAGCGGCCGCGGTCAc                                                                                                          >  2:805057/1‑150 (MQ=17)
                                             gTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACg                                                                                                         <  1:1326865/150‑1 (MQ=255)
                                              tGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACgg                                                                                                        <  1:1456096/150‑1 (MQ=255)
                                               gAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGt                                                                                                       <  2:1541292/150‑1 (MQ=255)
                                                 gCTTGTTATTACGGGTGATAATGCCACAGTGAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGTAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTAt                                                                                                     >  1:831470/1‑150 (MQ=255)
                                                     gTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGat                                                                                                 <  2:1356381/150‑1 (MQ=255)
                                                      ttattaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAACGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGata                                                                                                <  1:2020777/150‑1 (MQ=255)
                                                      ttattaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGata                                                                                                <  1:640267/150‑1 (MQ=255)
                                                      ttattaCGGGGGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGata                                                                                                <  1:1111966/150‑1 (MQ=255)
                                                       tattaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGatat                                                                                               >  1:1989450/1‑150 (MQ=255)
                                                        attaGGGGTTATAATGCCAGAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAATATTCCACCTGTACGGAAATCAACGGTAATATCTGGAAAGTGATTCAGGACTGCGATCTGGATGTCAGCGGCGGCGGTGACGGTATTGATATc                                                                                              <  2:115233/150‑1 (MQ=17)
                                                        attaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATc                                                                                              <  1:1119893/150‑1 (MQ=255)
                                                        attaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATc                                                                                              <  1:1455217/150‑1 (MQ=255)
                                                         ttaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCa                                                                                             >  1:83220/1‑150 (MQ=255)
                                                         ttaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCa                                                                                             >  1:794902/1‑150 (MQ=255)
                                                           aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCACCGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCAcc                                                                                           >  1:190121/1‑150 (MQ=255)
                                                            cGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCg                                                                                          >  2:1725984/1‑150 (MQ=255)
                                                             gggTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCgg                                                                                         <  1:570139/150‑1 (MQ=255)
                                                               gTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGTAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCggcg                                                                                       <  2:1989450/150‑3 (MQ=255)
                                                                     aTGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGa                                                                                                                                   >  1:455977/1‑100 (MQ=255)
                                                                          acaGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACTGTATTGATATCACCGGAGACAGCGCGACg                                                                            >  1:1903547/1‑150 (MQ=255)
                                                                          acaGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACg                                                                            <  1:859453/150‑1 (MQ=255)
                                                                             gTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTg                                                                         <  2:501395/150‑1 (MQ=255)
                                                                                  caacaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATc                                                                                                                                 >  2:1277588/1‑89 (MQ=255)
                                                                                   aacaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAAc                                                                   <  2:1374714/150‑1 (MQ=255)
                                                                                   aacaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAAc                                                                   >  2:522639/1‑150 (MQ=255)
                                                                                   aacaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAAc                                                                   <  1:522639/150‑1 (MQ=12)
                                                                                     caaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACAGTATTGATATc                                                                                              >  1:617451/1‑121 (MQ=255)
                                                                                     caaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACAGTATTGATATc                                                                                              <  2:617451/121‑1 (MQ=12)
                                                                                       aTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAgcgca                                                                               <  2:2067549/134‑2 (MQ=255)
                                                                                         ggCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGc                                                             <  2:2048566/150‑1 (MQ=255)
                                                                                                ccaccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGa                                                      >  1:1910334/1‑150 (MQ=255)
                                                                                                 caccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGTAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATc                                                                                              <  2:2088797/109‑1 (MQ=12)
                                                                                                 caccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATc                                                                                              >  1:2088797/1‑109 (MQ=25)
                                                                                                   ccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGAccc                                                   >  2:922555/1‑149 (MQ=255)
                                                                                                   ccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCg                                                   <  1:395434/150‑1 (MQ=255)
                                                                                                   ccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCg                                                   >  2:1839538/1‑150 (MQ=255)
                                                                                                   ccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGAGCTGGATGTCAGCGGCGGCGGGCACGGTATTGATCTCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCg                                                   <  1:1107251/150‑1 (MQ=255)
                                                                                                   ccGGTGACGGTAAACATTCCACCGGCACCGAAATCAACGGTAATAACGGGAACGGGACGCAGGACGGCGAGCTGGATGTCAGCGGCGGCGGTCACGGTATGGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCg                                                   <  2:1520260/150‑1 (MQ=14)
                                                                                                    cGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCgg                                                                                         >  2:1386509/1‑111 (MQ=25)
                                                                                                           ggTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGCCAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGAt                                           >  1:88578/1‑150 (MQ=255)
                                                                                                               aaGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCgg                                       <  1:1817446/150‑1 (MQ=40)
                                                                                                                aGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGa                                      >  2:1223025/1‑150 (MQ=40)
                                                                                                                 gATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGGCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAg                                     <  2:1718980/150‑1 (MQ=255)
                                                                                                                 gATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATGCCGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAg                                     <  2:1856211/150‑1 (MQ=255)
                                                                                                                  cttCCCCCGGTACGGAAAGCAACGGTAATAACGGGAAAGTGAATCAGGACGCCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGGGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGt                                    <  2:1339673/149‑1 (MQ=11)
                                                                                                                  aTTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGt                                    <  1:1801821/150‑1 (MQ=39)
                                                                                                                        ccGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCAtcg                              >  1:1253980/1‑148 (MQ=38)
                                                                                                                        ccGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGAGACAGCGCAACGGTGGATAACAAGGGc                                                             >  1:1631706/1‑119 (MQ=12)
                                                                                                                        ccGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCAATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATgg                              <  2:1582491/150‑1 (MQ=255)
                                                                                                                         cGGTACGGAAATCAACGGTAATACCGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATggg                             <  2:74594/150‑1 (MQ=255)
                                                                                                                         cGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATggg                             <  1:666774/150‑1 (MQ=255)
                                                                                                                         cGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGGGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATggg                             <  1:33815/150‑1 (MQ=255)
                                                                                                                         cGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTAGTGATATCACCGCTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATggg                             <  1:103696/150‑1 (MQ=255)
                                                                                                                          ggTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGt                            <  1:1368701/150‑1 (MQ=255)
                                                                                                                              cGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGAGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATcc                        >  2:656634/1‑150 (MQ=21)
                                                                                                                                gAAATCAACGGTAATTACGGGAAAGGGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAg                      <  1:1306251/150‑1 (MQ=255)
                                                                                                                                gAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAgc                                                                                  >  2:1604614/1‑90 (MQ=18)
                                                                                                                                gAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAgc                                                                                  <  1:1604614/90‑1 (MQ=18)
                                                                                                                                   aTCAACGGTAATCACGGGACAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATc                   <  1:1310484/150‑1 (MQ=255)
                                                                                                                                    tCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCg                  <  1:1439866/150‑1 (MQ=255)
                                                                                                                                     cAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTg                                                                                                   <  2:1958969/68‑1 (MQ=40)
                                                                                                                                      aaCGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTgg                                                                                                                              >  1:340020/1‑40 (MQ=255)
                                                                                                                                       aCGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCggcg                                                                                       >  2:1038319/1‑76 (MQ=25)
                                                                                                                                       aCGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCggcg                                                                                       <  1:1038319/78‑3 (MQ=25)
                                                                                                                                           taataaCGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGat                                                                                                 >  1:996381/1‑64 (MQ=37)
                                                                                                                                             ataaCGGGAAAGTGATGCAGGACGGCGATCTGGACGTCAGCGGCGGCGGTCACGGTAGTGATCTCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACa         <  1:1256877/150‑1 (MQ=255)
                                                                                                                                              taaCGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCAACGGt                                                                          >  1:1717196/1‑84 (MQ=18)
                                                                                                                                              taaCGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCAACGGt                                                                          <  2:1717196/84‑1 (MQ=18)
                                                                                                                                                aCGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAgg      >  2:1053429/1‑150 (MQ=255)
                                                                                                                                                   ggAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCa   >  1:1497026/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCAT  >  NC_000913/1465990‑1466287

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: