Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,215,993 T→A 6.8% K167* (AAA→TAA)  yqjH ← putative siderophore interacting protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,215,9930TA 6.8% 332.9 / 8.5 114K167* (AAA→TAA) yqjHputative siderophore interacting protein
Reads supporting (aligned to +/- strand):  ref base T (42/64);  new base A (2/6);  total (44/70)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 115 total aligned reads displayed.

ATCATCGGCAGGGATTTGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTATCGCCCGGTTGC  >  NC_000913/3215844‑3216138
                                                                                                                                                     |                                                                                                                                                 
atcatcGGCAGGGATTTGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCAAACCAGCGCACTAACTTGCGGttt                                                                                                                                                   <  2:88968/150‑1 (MQ=255)
 tcatcGGCAGGGATTTGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTa                                                                                                                                                  >  1:1050875/1‑150 (MQ=255)
 tcagcGGCCGGGATTTGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGGGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTa                                                                                                                                                  <  2:316626/150‑1 (MQ=255)
   atcGGCAGGGATTTGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAAc                                                                                                                                                <  1:887322/150‑1 (MQ=255)
   atcGGCAGGGATTTGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAAc                                                                                                                                                <  1:895562/150‑1 (MQ=255)
    tcGGCAGGGATTTGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACg                                                                                                                                               <  2:452491/150‑1 (MQ=255)
    tcGGCAGGGATTTGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCACGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACg                                                                                                                                               <  2:830177/150‑1 (MQ=255)
     cGGCAGGGATTTTCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACgg                                                                                                                                              <  2:327397/150‑1 (MQ=255)
         aGGGATTTGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAg                                                                                                                                          <  2:585135/150‑1 (MQ=255)
         aGGGATTTGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAg                                                                                                                                          <  2:409028/150‑1 (MQ=255)
          gggATTTGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTATCCAGCGCACTATCTGTCGGTTTAACGGCAAGt                                                                                                                                         >  1:307002/1‑150 (MQ=255)
          gggATTTGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGt                                                                                                                                         >  1:23623/1‑150 (MQ=255)
          gggATTTGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGt                                                                                                                                         >  1:185428/1‑150 (MQ=255)
             aTTTGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTg                                                                                                                                      <  1:1469796/150‑1 (MQ=255)
              tttGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGGTTAACGGCAAGTTTGc                                                                                                                                     <  2:881738/150‑1 (MQ=255)
              tttGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCAGCATGTGCCAGCCATTCGATATTAAAACCATCAAGGGGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGc                                                                                                                                     <  1:1499915/150‑1 (MQ=255)
               ttGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCt                                                                                                                                    <  1:972625/150‑1 (MQ=255)
                tGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGTGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACAAACTTGCCGTTAAACCGCAAGttagtgc                                                                                                                                   <  2:1701311/150‑6 (MQ=255)
                tGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTc                                                                                                                                   >  2:1114271/1‑150 (MQ=255)
                 gCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCa                                                                                                                                  >  1:134617/1‑150 (MQ=255)
                  cATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCaa                                                                                                                                 >  1:626676/1‑150 (MQ=255)
                  cATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGAGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCCGTTAAACCGCAAGttagtgcgc                                                                                                                                 <  2:597836/150‑8 (MQ=255)
                    tCTGCGCCAGACTCGCATCTACCGCCTTCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTTACTTGCGGTTTAACGGCAAGTTTGCTCAACg                                                                                                                               <  1:1300457/150‑1 (MQ=255)
                    tCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACg                                                                                                                               <  1:282850/150‑1 (MQ=255)
                    tCTGCGCCAGACGCGCATCTACCGCCTGCGCATCATCGGCCAGCCAATCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACg                                                                                                                               <  1:1079011/150‑1 (MQ=255)
                     cTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAACCCATCAAGGTGCGCGAGATAATCCGGCCAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCTAGTTTGCTCAACGt                                                                                                                              >  1:290420/1‑150 (MQ=255)
                     cTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGt                                                                                                                              >  1:954531/1‑150 (MQ=255)
                       gcgcCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGTCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCGGACCGttt                                                                                                                            >  2:450756/1‑150 (MQ=255)
                        cgcCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTc                                                                                                                           <  1:297407/150‑1 (MQ=255)
                        cgcCAGACGCGCATCTACCGCCTGCTCATCATGGGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTc                                                                                                                           <  1:857068/150‑1 (MQ=255)
                         gcCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTcc                                                                                                                          <  2:1688444/150‑1 (MQ=255)
                         gcCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTcc                                                                                                                          <  2:1214668/150‑1 (MQ=255)
                         gcCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTcc                                                                                                                          >  1:1125436/1‑150 (MQ=255)
                          ccAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCa                                                                                                                         <  1:724620/150‑1 (MQ=255)
                           cAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAg                                                                                                                        <  2:1576076/150‑1 (MQ=255)
                                gcgcATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAggcggc                                                                                                                   <  1:768064/150‑1 (MQ=255)
                                 cgcATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAggcggcg                                                                                                                  >  1:954860/1‑150 (MQ=255)
                                 cgcATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAggcggcg                                                                                                                  <  1:795008/150‑1 (MQ=255)
                                     tCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCCCGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGgcgc                                                                                                              >  1:135921/1‑150 (MQ=255)
                                      cTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAGCCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCa                                                                                                             >  1:1412272/1‑150 (MQ=255)
                                       tACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCaa                                                                                                            >  2:1089040/1‑150 (MQ=255)
                                       tACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCaa                                                                                                            >  2:243454/1‑150 (MQ=255)
                                        aCCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAAt                                                                                                           <  1:1163587/150‑1 (MQ=255)
                                           gCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCa                                                                                                        <  2:1564348/150‑1 (MQ=255)
                                             cTGCTCATCATGTGCCAGCCATTCGATATTAACACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTAAACGGCAAGTTTGCTCAACGGTTCCAGGCGGCGGCGCAATGCAgg                                                                                                      <  2:84343/150‑1 (MQ=255)
                                              tGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTAAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGc                                                                                                     >  1:1371329/1‑150 (MQ=255)
                                                cTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCAt                                                                                                   >  1:1511997/1‑150 (MQ=255)
                                                       tgtgCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGAt                                                                                            <  2:134617/150‑1 (MQ=255)
                                                        gtgCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATTATCCTGCCAGGCGTTATCACGCACGCTATCCAGCTCACTAACTGGCGGTTTAACGGCAAGGGTGCTCAACGATTCCAGGCGGCGGCGCAATGGATGAATTCCGGAtt                                                                                           >  2:277401/1‑150 (MQ=255)
                                                          gccagccaTTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCa                                                                                         >  2:65520/1‑150 (MQ=255)
                                                            cagccaTTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATc                                                                                       <  2:1709297/150‑1 (MQ=255)
                                                             cgccaTTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGTTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCg                                                                                      <  1:1644391/149‑1 (MQ=255)
                                                              gccaTTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGc                                                                                     <  1:1320267/150‑1 (MQ=255)
                                                              cccaTTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATTCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGc                                                                                     <  1:454441/149‑1 (MQ=255)
                                                                caTTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAg                                                                                   <  2:626676/150‑1 (MQ=255)
                                                                caTTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAg                                                                                   <  1:1089040/150‑1 (MQ=255)
                                                                  ttCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGAc                                                                                 <  1:525844/150‑1 (MQ=255)
                                                                  ttCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGGTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGAc                                                                                 <  1:1544555/150‑1 (MQ=255)
                                                                    cGATATTATAACCATCAAGGTCCGAGAGATAAACCTGACAGGCGTTATCCCGCACGCTAACCAGCGCTCTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGAGGCGCAATGCAGGCATTCCGGATTCATCGCAGACAt                                                                               <  1:273369/150‑1 (MQ=255)
                                                                     gATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATa                                                                              >  2:1456935/1‑150 (MQ=255)
                                                                     gATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACAAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATa                                                                              <  2:1710733/150‑1 (MQ=255)
                                                                       tatTAAAACCATCAAGGTGCGCGAGATTATCCTGACAGGCGTTAACCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGAATCAACGCAGACATACa                                                                            >  2:203429/1‑150 (MQ=255)
                                                                       tatTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCATGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACa                                                                            <  1:404136/150‑1 (MQ=255)
                                                                         ttAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGc                                                                          <  2:954531/150‑1 (MQ=255)
                                                                          tAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCt                                                                         >  1:1337494/1‑150 (MQ=255)
                                                                           aaaaCCATCAAGGTGCGCGAGATAATCCTGACAGGCTTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTg                                                                        <  1:137919/150‑1 (MQ=255)
                                                                            aaaCCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGACATTCCGGATTCATCGCAGACATACAGCTGa                                                                       >  1:1702337/1‑150 (MQ=255)
                                                                               ccATAAAGGTGCGCGAGATAATCCTGACAGGCGATATCCCGCACGCTAACCAGCGCACTAACTTGCGGTATAACGGCTAGTTTGCTCAACGTATCCAGGCGGCGGCGCAATCAAGGCATTCCGGATTCTTCGCAGACATACAGCGGAttc                                                                    >  2:1141961/1‑148 (MQ=255)
                                                                                cATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATAcg                                                                   <  2:274319/150‑1 (MQ=255)
                                                                                  tCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATAcgcg                                                                 <  1:646366/150‑1 (MQ=255)
                                                                                     aGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTaa                                                              <  2:606040/150‑1 (MQ=255)
                                                                                      ggTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAAt                                                             >  1:770097/1‑150 (MQ=255)
                                                                                          cgcgAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTc                                                         >  1:316728/1‑150 (MQ=255)
                                                                                           gcgAGATAATCCTGACAGGCGGTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGGTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTcc                                                        <  2:1267401/150‑1 (MQ=255)
                                                                                              agaTAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCGAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACa                                                                                <  2:240121/123‑1 (MQ=255)
                                                                                              agaTAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGc                                                     >  1:1543298/1‑150 (MQ=255)
                                                                                               gaTAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGca                                                    <  2:185428/150‑1 (MQ=255)
                                                                                                   aTCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTATCGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGCCTTACAGCTGATACGCGTAATCTTCCGGcaccac                                                >  2:323778/1‑150 (MQ=255)
                                                                                                     ccTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCAccacca                                              >  1:451090/1‑150 (MQ=255)
                                                                                                        gACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCg                                           >  2:1194441/1‑150 (MQ=255)
                                                                                                           aGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAAc                                        >  1:382875/1‑150 (MQ=255)
                                                                                                           aGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAAc                                        <  1:1030140/150‑1 (MQ=255)
                                                                                                                 tATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACcgcgcg                                  <  2:1543298/150‑1 (MQ=255)
                                                                                                                  aacccGCACGCTAACCAGCGCACTATCTTGCGGGTTAACGGCAAGATTGCTCAACGTATCCAGGCTGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCgg                                 <  1:154392/148‑1 (MQ=255)
                                                                                                                    cccGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGAc                               <  1:968280/150‑1 (MQ=255)
                                                                                                                       gCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTg                            <  2:619192/150‑1 (MQ=255)
                                                                                                                          cGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCa                         >  1:374036/1‑150 (MQ=255)
                                                                                                                           gCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCAc                        <  1:429966/150‑1 (MQ=255)
                                                                                                                            cTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCAcc                       >  2:450361/1‑150 (MQ=255)
                                                                                                                               cccAGCGCACTAACATGCGGTGTCACGGCACGGTTGCTCAACGGTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTa                    <  2:290420/149‑1 (MQ=255)
                                                                                                                                 cAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAg                  <  2:1125436/150‑1 (MQ=255)
                                                                                                                                  aGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCTGCGGCGCAATGCAGGCATTCCGGCTTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGt                 <  2:1422571/150‑1 (MQ=255)
                                                                                                                                  aGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGt                 <  2:1711825/150‑1 (MQ=255)
                                                                                                                                      agcaaaCTTGCCGTTAAACCGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCAATCCGGATTCATCGCAGACATACAGCGGATACGCGTATTCTTCCGGCACCACCAGCGAACCGCGCGGACCTCCCGCCGTAAGTTTAt             >  1:203063/5‑150 (MQ=255)
                                                                                                                                          aaCTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTATCGcc         <  2:620255/150‑1 (MQ=255)
                                                                                                                                            cTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTATCGCCCg       >  1:1729165/1‑150 (MQ=255)
                                                                                                                                            cTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTATCGCCCg       <  2:902600/150‑1 (MQ=255)
                                                                                                                                            cTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGAATCATCGAAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGCCCCTGCCACCGTAAGTTTATCGCCCg       >  2:282810/1‑150 (MQ=255)
                                                                                                                                                cGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCCCCGTAAGGTTATCGCCCGGTTg   >  1:313103/1‑150 (MQ=255)
                                                                                                                                                 ggTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTATCGCCCGGTTgc  >  1:1658582/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
ATCATCGGCAGGGATTTGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTATCGCCCGGTTGC  >  NC_000913/3215844‑3216138

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