Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,466,201 T→C 6.6% pseudogene (810/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,2010TC 6.6% 226.8 / 7.0 99pseudogene (810/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base T (24/61);  new base C (3/3);  total (33/66)
Fisher's exact test for biased strand distribution p-value = 3.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.43e-01

GTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTG  >  NC_000913/1466053‑1466349
                                                                                                                                                    |                                                                                                                                                    
gTGATAATGCCACAGTAAACAACACTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGGACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTg                                                                                                                                                     <  1:842994/150‑1 (MQ=255)
     aaTGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGGGACAgc                                                                                                                                                >  1:818930/1‑150 (MQ=255)
          cacaGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGAc                                                                                                                                           >  1:1611654/1‑150 (MQ=255)
           acaGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACTGTATTGATATCACCGGAGACAGCGCGACg                                                                                                                                          >  1:1903547/1‑150 (MQ=255)
           acaGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACg                                                                                                                                          <  1:859453/150‑1 (MQ=255)
              gTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTg                                                                                                                                       <  2:501395/150‑1 (MQ=255)
               tAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGAAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTgg                                                                                                                                      <  2:1438613/150‑1 (MQ=255)
                    aacaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAAc                                                                                                                                 <  2:1374714/150‑1 (MQ=255)
                    aacaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAAc                                                                                                                                 >  2:522639/1‑150 (MQ=255)
                    aacaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAAc                                                                                                                                 <  1:522639/150‑1 (MQ=12)
                        aTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAgg                                                                                                                             <  1:929642/150‑1 (MQ=255)
                        aTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAgcgca                                                                                                                                             <  2:2067549/134‑2 (MQ=255)
                        aTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAgcgca                                                                                                                                             >  1:2067549/1‑133 (MQ=255)
                          ggCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGc                                                                                                                           <  2:2048566/150‑1 (MQ=255)
                             aaaaCCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGTATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCAcc                                                                                                                        <  2:1237018/150‑1 (MQ=255)
                                 ccaccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGa                                                                                                                    >  1:1910334/1‑150 (MQ=255)
                                    ccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGAccc                                                                                                                 >  2:922555/1‑149 (MQ=255)
                                    ccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCg                                                                                                                 >  2:1839538/1‑150 (MQ=255)
                                    ccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCg                                                                                                                 <  1:395434/150‑1 (MQ=255)
                                    ccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGAGCTGGATGTCAGCGGCGGCGGGCACGGTATTGATCTCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCg                                                                                                                 <  1:1107251/150‑1 (MQ=255)
                                    ccGGTGACGGTAAACATTCCACCGGCACCGAAATCAACGGTAATAACGGGAACGGGACGCAGGACGGCGAGCTGGATGTCAGCGGCGGCGGTCACGGTATGGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCg                                                                                                                 <  2:1520260/150‑1 (MQ=14)
                                         gACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCAcc                                                                                                            <  2:1745932/150‑1 (MQ=255)
                                          aCGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCg                                                                                                           <  1:719133/150‑1 (MQ=255)
                                            ggTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGCCAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGAt                                                                                                         >  1:88578/1‑150 (MQ=255)
                                                aaGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCgg                                                                                                     <  1:1817446/150‑1 (MQ=40)
                                                 aGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGa                                                                                                    >  2:1223025/1‑150 (MQ=40)
                                                  gATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAg                                                                                                   <  2:2027316/150‑1 (MQ=39)
                                                  gATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAg                                                                                                   <  2:879677/150‑1 (MQ=39)
                                                  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)
                                                        aCCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCGGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCCCCGAGCCGGAGTCCATg                                                                                             >  1:179168/1‑150 (MQ=255)
                                                         ccGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCAtcg                                                                                            >  1:1253980/1‑148 (MQ=38)
                                                         ccGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGAGACAGCGCAACGGTGGATAACAAGGGc                                                                                                                           >  1:1631706/1‑119 (MQ=12)
                                                         ccGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCAATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATgg                                                                                            <  2:1582491/150‑1 (MQ=255)
                                                          cGGTACGGAAATCACCGGTAATAACGCGAAAGTGATTCAGGACGGCGAGCTGGAGGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATggg                                                                                           <  2:599408/150‑1 (MQ=255)
                                                          cGGTACGGAAATCAACGGTAATACCGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATggg                                                                                           <  2:74594/150‑1 (MQ=255)
                                                          cGGTACGGAAATCAACGGTAATAACGGTAAAGTGATTCAGGACGGCGATCTGGATGTCAGAGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATggg                                                                                           <  1:1042838/150‑1 (MQ=255)
                                                          cGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATggg                                                                                           <  1:597957/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)
                                                              aCGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGAGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATc                                                                                                        >  2:1372024/1‑133 (MQ=17)
                                                               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)
                                                                       aaCGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGAc                                                                              <  2:1659333/150‑1 (MQ=255)
                                                                        aCGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACg                                                                             <  2:2094493/150‑1 (MQ=255)
                                                                              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)
                                                                                    ggAAAGAGATTCAGGACGGCGATCAGGAACTCAGCCGCGGCGGTCACGGTATAGAAATCACCGGTGACAGCGCGACGGTGGATAACAAGGTCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCa                                                                 <  1:1303168/150‑1 (MQ=14)
                                                                                     gAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACTGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCAt                                                                <  1:1483066/150‑1 (MQ=255)
                                                                                        aGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGt                                                             >  2:1941259/1‑150 (MQ=255)
                                                                                            atgcAGGACGGCGATCTGGATGTCAGCGGCGGCGGGCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAAt                                                         <  1:611247/147‑1 (MQ=255)
                                                                                            aTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAAt                                                         >  2:785231/1‑150 (MQ=255)
                                                                                            aTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAAt                                                         <  2:705210/150‑1 (MQ=255)
                                                                                             ttCATGACGGCGATCTGGATGTCAGCGGCGGAGGTCACGGTATTGATATCACCGGTTACAGCGCGACGGTGGATAACAACGGCACCATGACAGTCACCGATCCGGAGTCCATGGGTATCAAGATCGACGGTGACAAGGCCATCGTCAATa                                                        <  1:300090/150‑1 (MQ=14)
                                                                                               cAGGACGGCGATCTGGATGTCAGCGGCGTCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGTATAACAAGCGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAAc                                                      <  1:362436/150‑1 (MQ=255)
                                                                                                aGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACg                                                     <  2:88578/150‑1 (MQ=255)
                                                                                                   ccGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATTACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAg                                                  <  2:1526013/149‑1 (MQ=255)
                                                                                                     ggCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGc                                                <  2:1561841/150‑1 (MQ=255)
                                                                                                        gATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATcgg                                                                                           >  1:570453/1‑104 (MQ=17)
                                                                                                        gATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCgag                                             <  1:432065/150‑1 (MQ=255)
                                                                                                        gATCTGGATGGCAGCGTCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCgag                                             <  2:1278566/150‑1 (MQ=255)
                                                                                                            tGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCa                                         >  1:631346/1‑150 (MQ=255)
                                                                                                            tGGATGGCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCa                                         <  2:1435030/150‑1 (MQ=255)
                                                                                                             ggATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCAc                                        >  2:927759/1‑150 (MQ=255)
                                                                                                             ggATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCAc                                        <  2:1339321/150‑1 (MQ=255)
                                                                                                             ggATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCAc                                        <  1:903565/150‑1 (MQ=255)
                                                                                                              gATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGAGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGt                                                                                                  >  1:853746/1‑91 (MQ=18)
                                                                                                               aTGTCAGCGGCGGCGGTCACGGTATTTATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCa                                      <  1:706290/150‑1 (MQ=255)
                                                                                                               aTGCCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCa                                      <  2:459361/150‑1 (MQ=255)
                                                                                                                 gTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGTCACCATCACCGTCACCGATCCGGAGTCCAGGGGTATCCAGATCGACGGTGACAACGCCAGCGACAATAACGACGGCGAGAGCACCATc                                    >  2:1475370/1‑150 (MQ=14)
                                                                                                                 gTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATc                                    <  1:717928/150‑1 (MQ=255)
                                                                                                                  tCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGGGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCa                                   >  1:1839499/1‑150 (MQ=255)
                                                                                                                    aGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCAcc                                 >  1:1564107/1‑150 (MQ=255)
                                                                                                                      cggcggcggTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCaa                               <  1:742401/150‑1 (MQ=255)
                                                                                                                        gcggcggTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCAGGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGTGAGAGCACCATCACCAACg                             >  1:1307971/1‑150 (MQ=255)
                                                                                                                         cggcggTCACGGTATTGATATCACCGGTGCCAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACgg                            <  2:1928860/150‑1 (MQ=255)
                                                                                                                         cggcggTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACgg                            <  2:1431295/150‑1 (MQ=255)
                                                                                                                          ggcggTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGt                           <  1:1242586/150‑1 (MQ=255)
                                                                                                                            cggTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTgg                         <  1:660754/150‑1 (MQ=255)
                                                                                                                               tCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCgg                                                                                                     >  2:485695/1‑71 (MQ=39)
                                                                                                                                cACGGTTTTGATATCACCGGTGACAGCGCGACTGTGGATAACAATGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTggcacc                     <  1:515912/150‑1 (MQ=255)
                                                                                                                                 aCGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTggcaccg                    <  1:456513/150‑1 (MQ=255)
                                                                                                                                      aTTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTggcaccggcacc               <  1:883218/150‑1 (MQ=255)
                                                                                                                                            atCACCGGTGACAGCGCGACGGTGGATAACTAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGAGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCCCCGGCACGCCGAtt         >  1:1111499/1‑150 (MQ=14)
                                                                                                                                               aCCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAAt      >  2:21873/1‑150 (MQ=255)
                                                                                                                                                  ggTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGt   >  1:1172622/1‑150 (MQ=255)
                                                                                                                                                   gTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtg  <  2:1910334/150‑1 (MQ=255)
                                                                                                                                                   gTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtg  <  1:1083836/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTG  >  NC_000913/1466053‑1466349

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