Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,431,155 G→T 6.1% S87S (TCC→TCA tsaC ← threonylcarbamoyl‑AMP synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,431,1550GT 6.1% 326.4 / 5.9 97S87S (TCC→TCAtsaCthreonylcarbamoyl‑AMP synthase
Reads supporting (aligned to +/- strand):  ref base G (50/41);  new base T (2/4);  total (52/45)
Fisher's exact test for biased strand distribution p-value = 4.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

Only 100 of 101 total aligned reads displayed.

GAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCATCACTGCTGTTTCGCTATC  >  NC_000913/3431006‑3431304
                                                                                                                                                     |                                                                                                                                                     
gAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCgg                                                                                                                                                       <  1:94845/150‑1 (MQ=255)
    ccAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGaaa                                                                                                                                                   >  1:327755/1‑150 (MQ=255)
    ccAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGaaa                                                                                                                                                   <  2:68577/150‑1 (MQ=255)
         ggTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATg                                                                                                                                              <  1:707529/150‑1 (MQ=255)
          gTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATgg                                                                                                                                             >  2:92142/1‑150 (MQ=255)
          gTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATgg                                                                                                                                             >  2:97748/1‑150 (MQ=255)
             tACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGttt                                                                                                                                          <  1:6144/150‑1 (MQ=255)
             tACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGttt                                                                                                                                          >  1:204652/1‑150 (MQ=255)
             tACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGttt                                                                                                                                          >  2:613589/1‑150 (MQ=255)
              aCCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTc                                                                                                                                         >  1:87527/1‑150 (MQ=255)
              aCCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTc                                                                                                                                         >  1:908710/1‑150 (MQ=255)
                cATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCAc                                                                                                                                       <  2:204652/150‑1 (MQ=255)
                cATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCAc                                                                                                                                       >  1:1268264/1‑150 (MQ=255)
                   aaGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCt                                                                                                                                    >  2:1496166/1‑150 (MQ=255)
                    aGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGTGAAACCATTTTTTCCCGCTg                                                                                                                                   <  2:882466/150‑1 (MQ=255)
                    aGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTg                                                                                                                                   >  2:1132567/1‑150 (MQ=255)
                    aGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTg                                                                                                                                   <  2:364294/150‑1 (MQ=255)
                    aGCCTGGCACAAAACAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTg                                                                                                                                   <  1:1393573/150‑1 (MQ=255)
                     gCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCTAATCAAAGCGGCCCGTCAACCAGCGCTTTGTTTTCGCAGTCGCGTGAAATACAAAGGTGACAGGACCTGTCCAGCGTGAAACCATTTTTTCCCGCTgg                                                                                                                                  <  2:1063816/150‑2 (MQ=255)
                      ccTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGTGAAACCATTTTTTCCCGCTggc                                                                                                                                 <  1:708478/150‑3 (MQ=255)
                       cTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCATGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGTGAAACCATTTTTTCCCGCTggcc                                                                                                                                <  2:129993/150‑4 (MQ=255)
                       cTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGCCAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCAc                                                                                                                                >  2:1007180/1‑150 (MQ=255)
                        tGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACg                                                                                                                               <  1:1228884/150‑1 (MQ=255)
                           cacaAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCGCgtca                                                                                                                            >  1:1044359/1‑150 (MQ=255)
                           cacaAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACgtca                                                                                                                            >  1:1402392/1‑150 (MQ=255)
                           cacaAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACgtca                                                                                                                            >  2:978966/1‑150 (MQ=255)
                                 gCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAAc                                                                                                                      <  2:977976/150‑1 (MQ=255)
                                 gCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAAc                                                                                                                      <  1:910776/150‑1 (MQ=255)
                                     ccaccaACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATgg                                                                                                                  <  1:195891/150‑1 (MQ=255)
                                        ccaACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGgtgt                                                                                                               >  1:307869/1‑150 (MQ=255)
                                           aCGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACg                                                                                                                               <  2:244809/131‑1 (MQ=255)
                                           aCGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACg                                                                                                                               >  1:244809/1‑131 (MQ=255)
                                           aCGGATGGGCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCTCAGGCGCTGGAAAGACAAAGGTGACAGGTCCTGGCCAGCGGGAAAATATGGTTTCACGCTGCACGTCAGTCAAGATGGTGtcat                                                                                                            <  2:316869/150‑1 (MQ=255)
                                             ggATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGtcatca                                                                                                          <  1:978966/150‑1 (MQ=255)
                                              gATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCaa                                                                                                         >  1:1172821/1‑150 (MQ=255)
                                                    cGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAgg                                                                                                   >  2:93749/1‑150 (MQ=255)
                                                        gACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGttt                                                                                               >  2:284246/1‑150 (MQ=255)
                                                            cGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGTGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAgc                                                                                           <  1:139449/150‑1 (MQ=255)
                                                            cGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAgc                                                                                           <  1:97748/150‑1 (MQ=255)
                                                            cGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAgc                                                                                           <  1:92142/150‑1 (MQ=255)
                                                              tACAGCAAGCGAATCAAAGCGGCCCGTCAAGCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAgctg                                                                                         <  2:1055760/150‑1 (MQ=255)
                                                              tACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAgctg                                                                                         >  1:629459/1‑150 (MQ=255)
                                                               aCAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGTAAATATGGTTTCACGCTGCACGTCAGTCAACATGGGGTCCTCAATATAGGGTTTAAgctgc                                                                                        >  2:1389295/1‑150 (MQ=255)
                                                               aCAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAgctgc                                                                                        >  2:250683/1‑150 (MQ=255)
                                                                cAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAATACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAgctgct                                                                                       <  2:307869/150‑1 (MQ=255)
                                                                  gCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCg                                                                                     >  2:680434/1‑150 (MQ=255)
                                                                  gCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATCGTGTCATCAATATAGGGTTTAAGCTGCACg                                                                                     >  1:608166/1‑150 (MQ=255)
                                                                     aGCGAATCAAAGCGGCCCGTCAACCAGCGCGCTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGCCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCCTaa                                                                                  >  2:1289820/1‑150 (MQ=255)
                                                                      gCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCTGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTGAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAAt                                                                                 <  2:284311/150‑1 (MQ=255)
                                                                       cGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAAtt                                                                                >  1:912456/1‑150 (MQ=255)
                                                                          aTCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATAtgc                                                                             <  2:808702/150‑1 (MQ=255)
                                                                             aaaGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCTGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTtgctgc                                                                          <  2:216688/150‑1 (MQ=255)
                                                                             aaaGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTtgctgc                                                                          >  2:1497875/1‑150 (MQ=255)
                                                                              aaGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCg                                                                         <  1:131421/150‑1 (MQ=255)
                                                                                  ggCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTa                                                                     >  2:671336/1‑150 (MQ=255)
                                                                                  ggCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTa                                                                     <  1:1282653/150‑1 (MQ=255)
                                                                                     ccGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTaaaa                                                                  >  1:416338/1‑150 (MQ=255)
                                                                                            ccAGCGCGCTGTTGTCCCAGGCGCGGGAAATACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCACTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGccc                                                           >  2:1129491/1‑150 (MQ=255)
                                                                                              aGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGAGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCt                                                         >  1:7992/1‑150 (MQ=255)
                                                                                                      gttgtCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAAcc                                                 >  2:1488482/1‑150 (MQ=255)
                                                                                                       ttgtCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCg                                                <  1:680434/150‑1 (MQ=255)
                                                                                                       ttgtCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCg                                                <  1:1497875/150‑1 (MQ=255)
                                                                                                         gtCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGa                                              >  2:518727/1‑150 (MQ=255)
                                                                                                           cGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCg                                                <  2:1136321/146‑1 (MQ=255)
                                                                                                           cGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCg                                                >  1:1136321/1‑146 (MQ=255)
                                                                                                            gCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGc                                           >  1:375965/1‑150 (MQ=255)
                                                                                                            gCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGc                                           <  1:749130/150‑1 (MQ=255)
                                                                                                               ggCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGt                                        <  2:1269925/150‑1 (MQ=255)
                                                                                                                 cgcgGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTAt                                                      <  1:485419/134‑1 (MQ=255)
                                                                                                                 cgcgGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTAt                                                      >  2:485419/1‑134 (MQ=255)
                                                                                                                  gcgGGAAAGACAAAGGGTACAGGACCTGGCCAGCTGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGtttt                                     <  2:973601/150‑1 (MQ=255)
                                                                                                                   cgGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTTCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTa                                    <  1:1213572/150‑1 (MQ=255)
                                                                                                                     ggAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAAc                                  <  1:93749/150‑1 (MQ=255)
                                                                                                                     ggAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAAc                                  <  2:1402392/150‑1 (MQ=255)
                                                                                                                      gAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACt                                 >  1:39972/1‑150 (MQ=255)
                                                                                                                       aaaGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAATATCAGCCCCTTATCAACCGGACGCTGTTTTAACTc                                <  1:1047662/150‑1 (MQ=255)
                                                                                                                       aaaGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTc                                <  2:935372/150‑1 (MQ=255)
                                                                                                                        aaGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTca                               >  2:489049/1‑149 (MQ=255)
                                                                                                                          gACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCaa                             >  2:1027445/1‑150 (MQ=255)
                                                                                                                            cAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACa                           >  1:1116925/1‑150 (MQ=255)
                                                                                                                            cAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACa                           >  2:733478/1‑150 (MQ=255)
                                                                                                                                  tGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCa                     <  1:875738/150‑1 (MQ=255)
                                                                                                                                  tGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCa                     <  2:1005871/150‑1 (MQ=255)
                                                                                                                                   gACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCAt                    >  2:805146/1‑150 (MQ=255)
                                                                                                                                   gACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCAt                    <  1:613589/150‑1 (MQ=255)
                                                                                                                                      aGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCATCAc                 >  1:1134486/1‑150 (MQ=255)
                                                                                                                                       ggACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCGCCTTATCAACCGGACGCTGTTTTCACTCCGACAGTCGCATCagt                >  2:1053930/1‑148 (MQ=255)
                                                                                                                                       ggACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCATCAct                >  2:313982/1‑150 (MQ=255)
                                                                                                                                          ccTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCATCActgct             >  2:1526230/1‑150 (MQ=255)
                                                                                                                                           cTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCATCActgctg            >  1:584741/1‑150 (MQ=255)
                                                                                                                                            tGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATTGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCATCACTGCTGt           >  2:1350968/1‑150 (MQ=255)
                                                                                                                                            cgtgcAGCGTGAAACCATTTTTTCCCGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCATCACTGCTGt           >  1:775977/5‑150 (MQ=255)
                                                                                                                                             ggCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCATCACTGCTGtt          >  1:516254/1‑150 (MQ=255)
                                                                                                                                              tgcAGCGTGAAACCATTTTTTCCCGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCATCACTGCTGttt         >  1:549029/3‑150 (MQ=255)
                                                                                                                                               ccAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCATCACTGCTGTTTc        >  2:944783/1‑150 (MQ=255)
                                                                                                                                                cAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCATCACTGCTGTTTCg       <  2:908710/150‑1 (MQ=255)
                                                                                                                                                cAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCATCACTGCTGTTTCg       <  2:327755/150‑1 (MQ=255)
                                                                                                                                                     ggAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAAc                                  <  2:1312290/118‑1 (MQ=255)
                                                                                                                                                     ggAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAAc                                  >  1:1312290/1‑118 (MQ=255)
                                                                                                                                                     ggAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCATCACTGCTGTTTCGCTATc  >  2:311791/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCATCACTGCTGTTTCGCTATC  >  NC_000913/3431006‑3431304

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: