Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,462,678 T→A 6.2% intergenic (+32/‑334) fadL → / → yfdF long‑chain fatty acid outer membrane channel/bacteriophage T2 receptor/protein YfdF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,462,6780TA 6.2% 351.7 / 6.8 108intergenic (+32/‑334)fadL/yfdFlong‑chain fatty acid outer membrane channel/bacteriophage T2 receptor/protein YfdF
Reads supporting (aligned to +/- strand):  ref base T (57/44);  new base A (2/5);  total (59/49)
Fisher's exact test for biased strand distribution p-value = 2.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 108 total aligned reads displayed.

CGTTGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCA  >  NC_000913/2462532‑2462826
                                                                                                                                                  |                                                                                                                                                    
cGTTGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGc                                                                                                                                                   <  2:19077/150‑1 (MQ=255)
 gTTGGTGTTTCTTATATGCACGGTGAGAGCGTGAAAATTAACGTAGGCCGATAGCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTTACTAGGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCa                                                                                                                                                  <  2:598631/150‑1 (MQ=255)
     gtgtTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTg                                                                                                                                              <  1:510895/150‑1 (MQ=255)
     gtgtTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTg                                                                                                                                              >  2:879273/1‑150 (MQ=255)
     gtgtTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTg                                                                                                                                              >  2:884566/1‑150 (MQ=255)
       gtTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGAc                                                                                                                                            <  1:1242844/150‑1 (MQ=255)
         ttgttATATGCACGGTCATATCGTGAAAATTAACGTAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCTGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCTATGCAGGTGACtt                                                                                                                                          <  2:544386/147‑1 (MQ=255)
         ttCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCTATGCAGGTGACtt                                                                                                                                          <  2:1164835/150‑1 (MQ=255)
          tCTTATATGCACGGTCATAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACttt                                                                                                                                         <  1:251473/150‑1 (MQ=255)
          tCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCTATGCAGGTGACttt                                                                                                                                         <  1:778596/150‑1 (MQ=255)
          tCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACttt                                                                                                                                         >  1:184820/1‑150 (MQ=255)
          tCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACttt                                                                                                                                         <  2:530563/150‑1 (MQ=255)
           cTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTa                                                                                                                                        <  2:1208449/150‑1 (MQ=255)
             tataTGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAAc                                                                                                                                      >  1:1221166/1‑150 (MQ=255)
              atatGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCTATGCAGGTGACTTTatcc                                                                                                                                     <  1:706813/150‑4 (MQ=255)
               tatGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTc                                                                                                                                    >  1:243973/1‑150 (MQ=255)
               tatGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTAGATAAAGTCACCTGCATAGCAGGTGACTTTAACTc                                                                                                                                    >  2:342628/1‑150 (MQ=255)
                atGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTcc                                                                                                                                   <  1:1135531/150‑1 (MQ=255)
                 tGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTccc                                                                                                                                  >  1:1018356/1‑150 (MQ=255)
                  gCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTcccc                                                                                                                                 >  1:88558/1‑150 (MQ=255)
                      ggTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGAACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCAAAGCAGGTGACTTTAACTCCCCCACt                                                                                                                             <  1:1356079/150‑1 (MQ=255)
                        gcAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTc                                                                                                                           <  1:680267/149‑1 (MQ=255)
                           gagCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCAcc                                                                                                                        <  1:1344066/150‑1 (MQ=255)
                            agCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCg                                                                                                                       <  1:1006908/150‑1 (MQ=255)
                              cGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGaa                                                                                                                     >  2:649880/1‑150 (MQ=255)
                                   aaaTTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAagtag                                                                                                                <  1:1229350/150‑1 (MQ=255)
                                   aaaTTAACGAAGACACATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAagtag                                                                                                                >  2:147402/1‑150 (MQ=255)
                                    aaTTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTAAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAagtagt                                                                                                               >  1:113777/1‑150 (MQ=255)
                                     aTTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGaa                                                                                                                     <  2:125324/143‑1 (MQ=255)
                                     aTTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGaa                                                                                                                     >  1:125324/1‑143 (MQ=255)
                                           gAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTAACCGAAAGTAGTGTCCTCa                                                                                                        >  1:714860/1‑150 (MQ=255)
                                             aGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCAtt                                                                                                      <  1:633013/150‑1 (MQ=255)
                                               gCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGc                                                                                                    <  1:649880/150‑1 (MQ=255)
                                                cccATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCt                                                                                                   >  1:776822/1‑150 (MQ=255)
                                                    tACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTAcc                                                                                               >  1:659506/1‑150 (MQ=255)
                                                       cAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTAcctcc                                                                                            <  2:1534216/150‑1 (MQ=255)
                                                       cAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTAcctcc                                                                                            >  1:1449667/1‑150 (MQ=255)
                                                       cAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTCGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTAcctcc                                                                                            <  1:342628/150‑1 (MQ=255)
                                                         gTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTg                                                                                          >  2:454464/1‑150 (MQ=255)
                                                         gTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTg                                                                                          <  2:1009602/150‑1 (MQ=255)
                                                          ttCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGa                                                                                         <  1:1013413/150‑1 (MQ=255)
                                                           tCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGAAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAg                                                                                        >  1:1514358/1‑150 (MQ=255)
                                                              aGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGttt                                                                                     >  1:1112778/1‑150 (MQ=255)
                                                              aGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGttt                                                                                     <  1:1413251/150‑1 (MQ=255)
                                                                 cTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGc                                                                                  <  2:433011/150‑1 (MQ=255)
                                                                 cTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGc                                                                                  <  1:144539/150‑1 (MQ=255)
                                                                  tGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCa                                                                                 <  1:138509/150‑1 (MQ=255)
                                                                  tGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCa                                                                                 >  1:1377935/1‑150 (MQ=255)
                                                                   gAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCaa                                                                                >  2:261980/1‑150 (MQ=255)
                                                                    aaGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCaaa                                                                               <  1:209492/150‑1 (MQ=255)
                                                                     aGGTAAAGCCTGGCTGTTCGGTACTTACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAAc                                                                              <  1:261980/150‑1 (MQ=255)
                                                                      ggTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACa                                                                             >  2:370526/1‑150 (MQ=255)
                                                                          aaGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCt                                                                         >  2:621607/1‑150 (MQ=255)
                                                                                  cTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAg                                                                 >  1:1329196/1‑150 (MQ=255)
                                                                                        ggTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGcc                                                           >  1:276956/1‑150 (MQ=255)
                                                                                         gTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGccc                                                          <  1:219627/150‑1 (MQ=255)
                                                                                         gTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGccc                                                          <  2:243973/150‑1 (MQ=255)
                                                                                           aCTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCtt                                                        >  1:1288532/1‑150 (MQ=255)
                                                                                           aCTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTACTGAGTTTTGAAAACAGACTGATGGCAGCTTGCCCtt                                                        >  2:64008/1‑150 (MQ=255)
                                                                                            cTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCttt                                                       >  2:1243777/1‑150 (MQ=255)
                                                                                              aaCTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCa                                                     <  1:535638/150‑1 (MQ=255)
                                                                                              aaCTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCa                                                     <  1:381396/150‑1 (MQ=255)
                                                                                                cTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCaaa                                                   >  2:822379/1‑150 (MQ=255)
                                                                                                   tAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCa                                                >  1:118705/1‑150 (MQ=255)
                                                                                                   ggcgaaCGCGTTATCAGAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCa                                                >  1:825405/6‑150 (MQ=255)
                                                                                                    aaCTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCaa                                               <  2:1100828/150‑1 (MQ=255)
                                                                                                     cgaaCGCGTTATCAGAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTAAAATCAAt                                              >  1:559378/4‑150 (MQ=255)
                                                                                                      cTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATa                                             >  1:119340/1‑150 (MQ=255)
                                                                                                       tACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATaa                                            >  2:1264265/1‑150 (MQ=255)
                                                                                                        aCGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGATAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAg                                           >  2:1476259/1‑150 (MQ=255)
                                                                                                        aCGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAg                                           >  1:786938/1‑150 (MQ=255)
                                                                                                         cgcgTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAAATCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGc                                          >  2:369501/1‑150 (MQ=255)
                                                                                                          gcgTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCg                                         >  1:1395557/1‑150 (MQ=255)
                                                                                                                     aCGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAACTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGa                              >  1:1312452/1‑150 (MQ=255)
                                                                                                                        cgTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACaa                           >  2:1233392/1‑150 (MQ=255)
                                                                                                                        cgTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACaa                           >  2:1477334/1‑150 (MQ=255)
                                                                                                                        cgTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACaa                           >  2:63266/1‑150 (MQ=255)
                                                                                                                         gTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACaaa                          <  2:1018356/150‑1 (MQ=255)
                                                                                                                           tCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAAcc                        <  2:866853/150‑1 (MQ=255)
                                                                                                                            cGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAATCGGTATCCACAGACAAACCt                       >  2:403981/1‑150 (MQ=255)
                                                                                                                             gCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTg                      <  2:1395557/150‑1 (MQ=255)
                                                                                                                               cTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCGGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGaa                    >  1:306034/1‑150 (MQ=255)
                                                                                                                               cTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGaa                    <  1:884566/150‑1 (MQ=255)
                                                                                                                               cTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGaa                    <  1:879273/150‑1 (MQ=255)
                                                                                                                                 ggATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGCCAAACCTGAAgg                  >  1:583212/1‑150 (MQ=255)
                                                                                                                                 ggATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAgg                  <  2:466678/150‑1 (MQ=255)
                                                                                                                                 ggATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAgg                  >  1:291389/1‑150 (MQ=255)
                                                                                                                                   ttAAAGTCACCTGCTATGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGa                >  1:560761/2‑150 (MQ=255)
                                                                                                                                   aTAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGa                <  2:1221166/150‑1 (MQ=255)
                                                                                                                                    tAAAGTCACCTGCTATGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGaa               >  1:699527/1‑150 (MQ=255)
                                                                                                                                    tAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGaa               >  1:91891/1‑150 (MQ=255)
                                                                                                                                     aaaGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGaaa              <  1:1477334/150‑1 (MQ=255)
                                                                                                                                      aaGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCGCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAg             >  1:913813/1‑150 (MQ=255)
                                                                                                                                       aGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAgg            >  2:431760/1‑150 (MQ=255)
                                                                                                                                          cACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCat         >  2:1533813/1‑150 (MQ=255)
                                                                                                                                            ccTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCatat       >  1:346473/1‑150 (MQ=255)
                                                                                                                                            ccTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCatat       >  2:1490906/1‑150 (MQ=255)
                                                                                                                                            ccTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCatat       >  1:340996/1‑150 (MQ=255)
                                                                                                                                              tGCATAGCAGGTGACATTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATAGCTTACCTCCTGAGTTATGCAAACAGCCTGTAGGCAGCTTGCCCTATCAAATCAATAAGCGGTATCCACAGACAAACCAGAAGGGAAAAGCATAttt     >  2:8155/1‑150 (MQ=255)
                                                                                                                                                 aTAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCa  <  1:994054/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CGTTGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCA  >  NC_000913/2462532‑2462826

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: