Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,214,433 A→T 5.5% N52Y (AAT→TAT) ‡ metA → homoserine O‑succinyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,214,4330AT 5.5% 1053.9 / 30.0 320N52Y (AAT→TAT) ‡metAhomoserine O‑succinyltransferase
Reads supporting (aligned to +/- strand):  ref base A (140/161);  new base T (11/7);  total (151/169)
Fisher's exact test for biased strand distribution p-value = 3.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

Only 100 of 331 total aligned reads displayed.

GATTCGTGTGCCGGACGAGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGAACAACTTCTACTGTAACTTTGAAGATATTCAGGATCAGAACTTTGAC  >  NC_000913/4214285‑4214582
                                                                                                                                                    |                                                                                                                                                     
 aTTCGTGTGCCGGACGAGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAAt                                                                                                                                                     <  2:1908407/150‑1 (MQ=255)
         gCCGGACGAGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGATTTTCAGTTTca                                                                                                                                             <  2:1798386/150‑2 (MQ=255)
            ggACGAGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTgcg                                                                                                                                          <  2:328266/150‑1 (MQ=255)
             gACGAGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGATTTTCAGTTTcaatct                                                                                                                                         <  1:799452/150‑6 (MQ=255)
             gACGAGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTgcgc                                                                                                                                         >  2:682121/1‑150 (MQ=255)
                gAGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCTGTTTCTGCGCCTg                                                                                                                                      <  2:555851/150‑1 (MQ=255)
                gAGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTg                                                                                                                                      >  2:1415395/1‑150 (MQ=255)
                 aGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGc                                                                                                                                     >  1:1434409/1‑150 (MQ=255)
                    tACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGCTCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCttt                                                                                                                                  >  2:390465/1‑150 (MQ=255)
                    tACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCttt                                                                                                                                  >  1:1969315/1‑150 (MQ=255)
                    tACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTCCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAACAACATTGAAACTGAAAATCAGTTTATGCGGCTGCttt                                                                                                                                  >  1:494002/1‑150 (MQ=255)
                     aCCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTc                                                                                                                                 >  1:1390387/1‑150 (MQ=255)
                     aCCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTc                                                                                                                                 <  2:634421/150‑1 (MQ=255)
                        cgccgTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCaaa                                                                                                                              <  1:321595/150‑1 (MQ=255)
                         gccgTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAAc                                                                                                                             <  1:248684/150‑1 (MQ=255)
                              cAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCAcc                                                                                                                        >  1:1270897/1‑150 (MQ=255)
                              cAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCAcc                                                                                                                        >  2:1063706/1‑150 (MQ=255)
                                    cTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCa                                                                                                                  >  1:604447/1‑150 (MQ=255)
                                    cTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCa                                                                                                                  <  1:1891362/150‑1 (MQ=255)
                                     ttGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAg                                                                                                                 <  1:430984/150‑1 (MQ=255)
                                          tGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGa                                                                                                            <  2:1440598/150‑1 (MQ=255)
                                            aagaaAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGata                                                                                                          <  1:173040/150‑1 (MQ=255)
                                               aaaaCGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTc                                                                                                       <  2:80261/150‑1 (MQ=255)
                                                  aCGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGc                                                                                                    <  2:1745848/150‑1 (MQ=255)
                                                    gTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTG                                                                                                  >  2:672919‑M1/1‑150 (MQ=255)
                                                        ttGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGC                                                                                              <  1:962629‑M1/150‑1 (MQ=255)
                                                         tgtgATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCG                                                                                             >  2:124422‑M1/1‑150 (MQ=255)
                                                            gatgaCAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCAt                                                                                          >  1:111036/1‑150 (MQ=255)
                                                             atgaCAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATc                                                                                         <  1:682121/150‑1 (MQ=255)
                                                             atgaCAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATc                                                                                         <  1:1666581/150‑1 (MQ=255)
                                                              tgaCAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCg                                                                                        >  1:442712/1‑150 (MQ=255)
                                                               gaCAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGa                                                                                       >  1:542348/1‑150 (MQ=255)
                                                                 cAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCAcc                                                                                                                        <  2:146857/115‑1 (MQ=255)
                                                                 cAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGAtt                                                                                     >  2:1219055/1‑150 (MQ=255)
                                                                 cAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGAtt                                                                                     >  1:1778901/1‑150 (MQ=255)
                                                                  accTTCTCGGGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGTAGAAGATTGACACTGAACAGCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTc                                                                                    <  2:1945853/148‑1 (MQ=255)
                                                                   aCTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTcc                                                                                   >  2:1109184/1‑150 (MQ=255)
                                                                   aCTTCTAGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGCTTCTGATCCTTAACCTGATGCAGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTcc                                                                                   >  2:1772631/1‑150 (MQ=255)
                                                                    cTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTccc                                                                                  >  2:578046/1‑150 (MQ=255)
                                                                    cTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCT                                                                                                   >  1:944101‑M1/1‑133 (MQ=255)
                                                                    cTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCT                                                                                                   <  2:944101‑M1/133‑1 (MQ=255)
                                                                       ctcGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTg                                                                               <  2:245876/150‑1 (MQ=255)
                                                                         cGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGaa                                                                             <  1:1607379/150‑1 (MQ=255)
                                                                             cGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCgc                                                                         <  2:1433962‑M2/150‑3 (MQ=255)
                                                                              gTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCgcg                                                                        >  1:1033868‑M2/1‑147 (MQ=255)
                                                                                cTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCgcgca                                                                      <  1:1408258‑M2/150‑6 (MQ=255)
                                                                                 tGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCgcgca                                                                      <  2:431592‑M2/149‑6 (MQ=255)
                                                                                  ggTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCgattcacgggaatcgatgcgcaac                                                                    <  2:1379594‑M1/150‑25 (MQ=255)
                                                                                  ggTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCgattcacgggaatcgatgcgcaac                                                                    <  1:1127876‑M1/150‑25 (MQ=255)
                                                                                  ggTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTCCGCCTGCTTTCAAACTCACCTTTGCAGCTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCgcgcaac                                                                    <  2:868205‑M2/150‑8 (MQ=255)
                                                                                   gTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGTCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCgcgcaaca                                                                   >  2:1779915‑M2/1‑142 (MQ=255)
                                                                                    tCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCgattcacgggaatcgatgcgcaacag                                                                  <  1:22038‑M1/150‑27 (MQ=255)
                                                                                     cAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCgattcacgggaatcgatgcgcaacagc                                                                 <  1:896989‑M1/150‑28 (MQ=255)
                                                                                     cAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCgattcacgggaatcgatgcgcaacagc                                                                 <  2:1554615‑M1/150‑28 (MQ=255)
                                                                                     cAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACg                                                                 >  2:1005251/1‑150 (MQ=255)
                                                                                     cAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACg                                                                 >  1:1932560/1‑150 (MQ=255)
                                                                                     cAGGAAATTCGTCCACTTAAGGGTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCgattcacgggaatcgatgcgcaacagc                                                                 <  1:1006655‑M1/150‑28 (MQ=255)
                                                                                       ggAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCgattcacgggaatcgatgcgcaacagctg                                                               <  1:655658‑M1/150‑30 (MQ=255)
                                                                                       ggAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGcc                                                               <  2:1828281/150‑1 (MQ=255)
                                                                                          aaTTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGc                                                            >  1:1839789/1‑150 (MQ=255)
                                                                                              cGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTCAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTAGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCagag                                                        >  2:358637/1‑150 (MQ=255)
                                                                                               gTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGc                                                       >  1:719786/1‑150 (MQ=255)
                                                                                                 ccACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCgattcacgggaatcgatgcgcaacagctgaatatcgacc                                                     <  1:180638‑M1/150‑40 (MQ=255)
                                                                                                   aCTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCt                                                   <  2:1529581/150‑1 (MQ=255)
                                                                                                    cTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTg                                                  <  1:558286/150‑1 (MQ=255)
                                                                                                    attAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGGCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTg                                                  <  2:1874453/149‑1 (MQ=255)
                                                                                                        aGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGaaca                                              >  1:45398/1‑150 (MQ=255)
                                                                                                         ggTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGc                                                            >  1:1506157/1‑135 (MQ=255)
                                                                                                              tGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGa                                                 <  2:1987787/141‑1 (MQ=255)
                                                                                                               gATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGAACAACTTCTa                                       >  1:1715770/1‑150 (MQ=255)
                                                                                                               gATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGAACAACTTCTa                                       >  2:1622735/1‑150 (MQ=255)
                                                                                                               gATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGAACAACTTCTa                                       >  2:1341388/1‑150 (MQ=255)
                                                                                                               gATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACATTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGAACAACTTCTa                                       >  2:250607/1‑150 (MQ=255)
                                                                                                               gATCCTTAAACTGATGCCGAAGAATATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGAACAACTTCTa                                       >  2:1489466/1‑150 (MQ=255)
                                                                                                                aTCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGAACAACTTCTAc                                      <  1:1063706/150‑1 (MQ=255)
                                                                                                                aTCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGAACAACTTCTAc                                      >  2:1364009/1‑150 (MQ=255)
                                                                                                                  ccTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGAACAACTTCTACTg                                    <  1:1448196/150‑1 (MQ=255)
                                                                                                                  ccTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCG                                                                                             <  2:761118‑M1/93‑1 (MQ=255)
                                                                                                                  ccTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCG                                                                                             >  1:761118‑M1/1‑93 (MQ=255)
                                                                                                                         cTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGAACAACTTCTACTGTAACttt                             <  1:297164/150‑1 (MQ=255)
                                                                                                                            ttGCCGAATAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCAGCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGAACAACTTCTACTGTAACTTTGaa                          <  2:847393/149‑1 (MQ=255)
                                                                                                                            aTGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGAACAACTTCTACTGTAACTTTGaa                          <  1:362919/150‑1 (MQ=255)
                                                                                                                             tGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGAACAACTTCTACTGTAACTTTGAAg                         <  1:1373417/150‑1 (MQ=255)
                                                                                                                               ccGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGAACAACTTCTACTGTAACTTTGAAg                         <  2:119689/148‑1 (MQ=255)
                                                                                                                               ccGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGAACAACTTCTACTGTAACTTTGAAg                         >  1:119689/1‑148 (MQ=255)
                                                                                                                                  aagaagATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGAACAACTTCTACTGTAACTTTGAAGATAtt                    <  1:1334452/150‑1 (MQ=255)
                                                                                                                                  aagaagATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGAACAACTTCTACTGTAACTTTGAAGATAtt                    >  1:1076915/1‑150 (MQ=255)
                                                                                                                                   agaagaTTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGAACAACTTCTACTGTAACTTTGAAGATATTc                   <  2:1969315/150‑1 (MQ=255)
                                                                                                                                   agaagaTTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGAACAACTTCTACTGTAACTTTGAAGATATTc                   >  2:1313174/1‑150 (MQ=255)
                                                                                                                                       gaTTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGAACAACTTCTACTGTAACTTTGAAGATATTCAGGa               >  2:368355/1‑150 (MQ=255)
                                                                                                                                       cgcagAAACTGATTTTCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGAACAACTTCTACTGTAACTTTGAAGATATTCAGGa               >  1:262306/5‑150 (MQ=255)
                                                                                                                                        gcagAAACTGATTTTCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGAACAACTTCTACTGTAACTTTGAAGATATTCAGGAt              >  1:725620/4‑150 (MQ=255)
                                                                                                                                        aTTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGAACAACTTCTACTGTAACTTTGAAGATATTCAGGAt              >  1:709404/1‑150 (MQ=255)
                                                                                                                                         cagAAACTGATTTTCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGAACAACTTCTACTGTAACTTTGAAGATATTCAGGATc             >  2:909110/3‑150 (MQ=255)
                                                                                                                                           gAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGAACAACTTCTACTGTAACTTTGAAGATATTCAGGATCAg           >  1:1471343/1‑150 (MQ=255)
                                                                                                                                              aCTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGAACAACTTCTACTGTAACTTTg                            >  1:864291/1‑130 (MQ=255)
                                                                                                                                                  aaaaTCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGAACAACTTCTACTGTAACTTTGAAGATATTCAGGATCAGAACTTTg    <  2:502882/150‑1 (MQ=255)
                                                                                                                                                  aaaaTCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGAACAACTTCTACTGTAACTTTGAAGATATTCAGGATCAGAACTTTg    <  2:1041235/150‑1 (MQ=255)
                                                                                                                                                    aaTCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGAACAACTTCTACTGTAACTTTGAAGATATTCAGGATCAGAACTTTGAc  <  1:1313174/150‑1 (MQ=255)
                                                                                                                                                    aaTCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGAACAACTTCTACTGTAACTTTGAAGATATTCAGGATCAGAACTTTGAc  <  2:172659/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GATTCGTGTGCCGGACGAGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGAACAACTTCTACTGTAACTTTGAAGATATTCAGGATCAGAACTTTGAC  >  NC_000913/4214285‑4214582

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