Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,825,938 A→C 5.1% intergenic (‑313/+17) spy ← / ← astE ATP‑independent periplasmic chaperone/succinylglutamate desuccinylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,825,9380AC 5.1% 359.3 / 6.9 112intergenic (‑313/+17)spy/astEATP‑independent periplasmic chaperone/succinylglutamate desuccinylase
Reads supporting (aligned to +/- strand):  ref base A (51/55);  new base C (4/2);  total (55/57)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 113 total aligned reads displayed.

TGATGGAGAGTGATGAGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTC  >  NC_000913/1825791‑1826087
                                                                                                                                                   |                                                                                                                                                     
tGATGGAGAGTGATGAGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGAtt                                                                                                                                                     <  1:1046382/150‑1 (MQ=255)
tGATGGAGAGTGATGAGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGAtt                                                                                                                                                     <  1:772253/150‑1 (MQ=255)
   tGGAGAGTTATGAGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCt                                                                                                                                                  <  2:720436/150‑1 (MQ=255)
        agTGATGAGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAAATCATTCTGCggg                                                                                                                                             <  2:76577/150‑1 (MQ=255)
         gTGATGAGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGt                                                                                                                                            <  2:1001598/150‑1 (MQ=255)
             tgaGAGAATATTACAACCCGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTCATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCTGAATGATTTCTGCGGGTAAGt                                                                                                                                        <  2:404319/150‑1 (MQ=255)
              gagagaATATTACAACACGATGATTTTGCAGACATTATGAAGAACTATAACGCATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTa                                                                                                                                       >  2:312891/1‑150 (MQ=255)
               agagaATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTAt                                                                                                                                      >  1:486457/1‑150 (MQ=255)
                 agaATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTa                                                                                                                                    <  1:1322862/150‑1 (MQ=255)
                         acaacaCGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAATGAAACTTAAAACCCCCAGAATGATTTATCCGGGTAAGTATTACGTTAttt                                                                                                                            >  2:914533/1‑150 (MQ=255)
                            acacGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATTAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTAtttttt                                                                                                                         <  1:1306164/150‑1 (MQ=255)
                            acacGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTAtttttt                                                                                                                         >  2:1345560/1‑150 (MQ=255)
                            acacGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTAtttttt                                                                                                                         >  2:273931/1‑150 (MQ=255)
                            acacGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTAtttttt                                                                                                                         <  2:1242632/150‑1 (MQ=255)
                            acacGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTAtttttt                                                                                                                         >  2:1287154/1‑150 (MQ=255)
                             cacGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCTGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTc                                                                                                                        <  2:503464/150‑1 (MQ=255)
                             cacGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTc                                                                                                                        <  2:1148574/150‑1 (MQ=255)
                              acGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGTTTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCg                                                                                                                       <  2:875370/150‑1 (MQ=255)
                              acGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCg                                                                                                                       <  1:235768/150‑1 (MQ=255)
                               cGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGa                                                                                                                      >  1:982949/1‑150 (MQ=255)
                                 atgatTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGc                                                                                                                    <  1:1000507/150‑1 (MQ=255)
                                 atgatTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAACCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGc                                                                                                                    <  2:786563/150‑1 (MQ=255)
                                  tgatTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCa                                                                                                                   >  2:873017/1‑150 (MQ=255)
                                    atTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCAtt                                                                                                                 <  1:1287154/150‑1 (MQ=255)
                                    atTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCAt                                                                                                                  <  2:1087169/149‑1 (MQ=255)
                                    atTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCAt                                                                                                                  >  1:1087169/1‑149 (MQ=255)
                                     ttttTCAGAGATTATGAAGAACTATACCGGATGAATTGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTa                                                                                                                <  2:643538/150‑1 (MQ=255)
                                         gCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATcc                                                                                                            >  1:467363/1‑150 (MQ=255)
                                            gagaTTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCcgc                                                                                                         <  1:1028760/150‑1 (MQ=255)
                                            gagaTTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCcgc                                                                                                         <  1:616315/150‑1 (MQ=255)
                                              gaTTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTGATCCcgcgc                                                                                                       >  1:1314949/1‑150 (MQ=255)
                                               aTTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCcgcgcg                                                                                                      <  1:587158/150‑1 (MQ=255)
                                                  aTGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTaa                                                                                                   >  1:1506234/1‑150 (MQ=255)
                                                    gaagaaCTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATc                                                                                                 >  2:1348154/1‑150 (MQ=255)
                                                          cTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACg                                                                                           <  1:1125430/150‑1 (MQ=255)
                                                          cTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACg                                                                                           >  1:1265657/1‑150 (MQ=255)
                                                              aCCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTAc                                                                                       >  2:1181460/1‑150 (MQ=255)
                                                                 ggATGACTGGTGATAAATAAAGCAAAAAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGAGAGTAATCCCAAAGCTACCaa                                                                                    >  2:1033793/1‑150 (MQ=255)
                                                                  gATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAAc                                                                                   <  1:473168/150‑1 (MQ=255)
                                                                     gACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGa                                                                                >  1:237019/1‑150 (MQ=255)
                                                                     gACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGa                                                                                <  2:402623/150‑1 (MQ=255)
                                                                     gACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGa                                                                                <  2:1272573/150‑1 (MQ=255)
                                                                       cTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCTGGTAAGTATTAGATTATTTATTCGAGCATTAATCACCCGAGTAATCCCAACGATACAAACGGAtt                                                                              >  1:1541505/1‑150 (MQ=255)
                                                                          gTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAgg                                                                           <  1:324427/150‑1 (MQ=255)
                                                                           tGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAggg                                                                          <  1:1382971/150‑1 (MQ=255)
                                                                            gATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGa                                                                         >  2:1170226/1‑150 (MQ=255)
                                                                            gATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGa                                                                         >  1:622856/1‑150 (MQ=255)
                                                                               aaataaaGCAAATTACCAGGATTAATCTGTATTATTTTATAAGAAAGCAACTTAATACCCGAAGAATGATTTCTGCGGGTATGTATTAGCTTATTTTTTCGAGCATTAATCACGCGCGTAATCCCAACGCTTCCAACGGATTAGGGAATa                                                                      >  2:1528682/1‑150 (MQ=255)
                                                                               aaataaaGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATa                                                                      >  2:207736/1‑150 (MQ=255)
                                                                                aataaaGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAAc                                                                                            <  2:783554/127‑1 (MQ=255)
                                                                                aataaaGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAAc                                                                                            >  1:783554/1‑127 (MQ=255)
                                                                                 ataaaGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAac                                                                    <  2:1275882/150‑1 (MQ=255)
                                                                                 ataaaGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAac                                                                    >  2:1631327/1‑150 (MQ=255)
                                                                                  taaaGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAaca                                                                   >  1:1619640/1‑150 (MQ=255)
                                                                                  taaaGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAaca                                                                   <  1:664500/150‑1 (MQ=255)
                                                                                    aaGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAAcaca                                                                 <  2:813345/150‑1 (MQ=255)
                                                                                      gCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATa                                                               >  2:291776/1‑150 (MQ=255)
                                                                                       cAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATAc                                                              <  2:1213544/150‑1 (MQ=255)
                                                                                        aaaTAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATAct                                                             >  1:82895/1‑150 (MQ=255)
                                                                                        aaaTAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATAct                                                             <  1:671596/150‑1 (MQ=255)
                                                                                         aaTAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATActc                                                            <  1:1276430/150‑1 (MQ=255)
                                                                                            aaCCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTAc                                                         >  2:1388213/1‑150 (MQ=255)
                                                                                             aCCAGGATTAAACTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGCTTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATAAATCCCGCGGGTAATCCCACCGATACCAACGGATTAGGGATTAACACATACTCTAca                                                        >  2:1087330/1‑150 (MQ=255)
                                                                                              ccAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTAcat                                                       >  1:559114/1‑150 (MQ=255)
                                                                                               cAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTAcatc                                                      <  1:1631327/150‑1 (MQ=255)
                                                                                               cAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTAcata                                                      >  1:85404/1‑149 (MQ=255)
                                                                                                aGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTAcatca                                                     >  1:1216921/1‑150 (MQ=255)
                                                                                                  taTTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATg                                                   <  1:1033793/149‑1 (MQ=255)
                                                                                                    ttAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATggg                                                 >  1:773893/1‑150 (MQ=255)
                                                                                                     tAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGt                                                <  2:920889/150‑1 (MQ=255)
                                                                                                         cTGTATTAATTTATAAGAAAGCAACTTAAAACACGCAGAATGATATCTGCGGGTAAGTAATACCTAAAATTTTCGAGAATAAATCCCGCCCGTAATCCCAAACCAACCAACGCATTACGGAATAACAAATACTCAACATCATgcctaacg                                            >  2:527269/1‑143 (MQ=255)
                                                                                                          tGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACgg                                           <  1:812131/150‑1 (MQ=255)
                                                                                                          tGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACgg                                           >  2:606161/1‑150 (MQ=255)
                                                                                                          tGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACgg                                           <  1:273931/150‑1 (MQ=255)
                                                                                                          tGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACgg                                           <  1:914533/150‑1 (MQ=255)
                                                                                                           gTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGt                                          >  1:516639/1‑150 (MQ=255)
                                                                                                            tATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTa                                         <  1:1139317/150‑1 (MQ=255)
                                                                                                             aTTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTaa                                        >  2:1633986/1‑150 (MQ=255)
                                                                                                             aTTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTaa                                        <  1:276147/150‑1 (MQ=255)
                                                                                                               tAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTaaaa                                      <  2:1216921/150‑1 (MQ=255)
                                                                                                                aaTTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAAc                                     >  1:1493338/1‑150 (MQ=255)
                                                                                                                aaTTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAAc                                     >  2:8696/1‑150 (MQ=255)
                                                                                                                 aTTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTAATTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACg                                    >  2:1588245/1‑150 (MQ=255)
                                                                                                                  tttATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAGAACGt                                   >  1:333527/1‑150 (MQ=255)
                                                                                                                      taAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCct                               >  1:469530/1‑150 (MQ=255)
                                                                                                                         gAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCtctc                            <  2:703537/150‑1 (MQ=255)
                                                                                                                            aGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGt                         <  2:589197/150‑1 (MQ=255)
                                                                                                                                aCTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTg                     >  2:915377/1‑150 (MQ=255)
                                                                                                                                 cTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTgc                    <  1:283851/150‑1 (MQ=255)
                                                                                                                                 cTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTgc                    <  2:26736/150‑1 (MQ=255)
                                                                                                                                  ttAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTgcg                   <  1:959394/150‑1 (MQ=255)
                                                                                                                                   acttACCCGCAGAAATCATTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTgcgc                  >  1:1042984/4‑150 (MQ=255)
                                                                                                                                    aaTACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGcc                 <  2:516639/150‑1 (MQ=255)
                                                                                                                                        cccGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCa             >  2:1056708/1‑150 (MQ=255)
                                                                                                                                        cccGCAGAAATCATTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCa             >  1:491658/1‑150 (MQ=255)
                                                                                                                                            cAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGt         >  1:1545462/1‑150 (MQ=255)
                                                                                                                                                aTGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCt     >  1:1502394/1‑150 (MQ=255)
                                                                                                                                                 tGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACg                                    <  2:1339316/118‑1 (MQ=255)
                                                                                                                                                 tGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACg                                    >  1:1339316/1‑118 (MQ=255)
                                                                                                                                                   aTTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTtc  <  2:486457/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TGATGGAGAGTGATGAGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTC  >  NC_000913/1825791‑1826087

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