Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,071,616 G→A 7.2% G47G (GGC→GGT)  yihS ← sulfoquinovose isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,071,6160GA 7.2% 298.2 / 7.9 97G47G (GGC→GGT) yihSsulfoquinovose isomerase
Reads supporting (aligned to +/- strand):  ref base G (43/47);  new base A (4/3);  total (47/50)
Fisher's exact test for biased strand distribution p-value = 7.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCATACCAGCCGCCATATTTTTTATCACGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGTTTCCTGTTCCAGCCAACGGTTGTGGCTTAGGGTGTTAAACCATTTCATTCCTTTTAT  >  NC_000913/4071469‑4071765
                                                                                                                                                   |                                                                                                                                                     
ggCATACCAGCCGCCATATTTTTTATCACGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTggc                                                                                                                                                     <  1:281249/150‑3 (MQ=255)
ggCATACCAGCCGCCATATTTTTTATCACGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGcc                                                                                                                                                     >  2:153507/1‑150 (MQ=255)
 gCATACCAGCCGCCATATTTTTTATCACGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGccc                                                                                                                                                    <  2:564938/150‑1 (MQ=255)
  cATACCAGCCGCCATATTTTTTATCACGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCa                                                                                                                                                   <  1:153920/150‑1 (MQ=255)
  cATACCAGCCGCCATATTTTTTATCACGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCa                                                                                                                                                   >  2:270823/1‑150 (MQ=255)
    tACCAGCCGCCATATTTTTTATCACGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCCACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCAtc                                                                                                                                                 >  1:289997/1‑150 (MQ=255)
    aaCCAGCCGCCATATTTTTTATCACGCAGTGCGCCGTTCATGGGTTTGATAGCGGGATCAACGATCGAGTCTGGGCCAGGTCGACCCATCGCCGGAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACATATTGGTGGACAtc                                                                                                                                                 <  2:120952/149‑1 (MQ=255)
      ccAGCCGCCATATTTTTTATCACGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCAtctc                                                                                                                                               <  2:1597674/150‑1 (MQ=255)
         gccgccATATTTTTTATCACGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGCACCCAtctgtgg                                                                                                                                            <  2:1344361/150‑5 (MQ=255)
          ccgccATATTTTTTATCACGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGCACCCAtctgtgga                                                                                                                                           <  1:638599/150‑6 (MQ=255)
            gccATATTTTTTATCACGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCttt                                                                                                                                         <  2:45561/150‑1 (MQ=255)
              cATATTTTTTATCACGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGa                                                                                                                                       <  2:761520/150‑1 (MQ=255)
               atatTTTTTATCACGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGAt                                                                                                                                      >  1:127620/1‑150 (MQ=255)
                tatTTTTTATCACGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTACGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGAtt                                                                                                                                     >  2:1505046/1‑150 (MQ=255)
                tatTTTTTATCACGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGAtt                                                                                                                                     >  2:596859/1‑150 (MQ=255)
                 atTTTTTATCACGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCCGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCACCATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGCttt                                                                                                                                    >  2:862432/1‑150 (MQ=255)
                 atTTTTTATCACGCAGTGCGACGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGACAGGTCGCCCAATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGAttt                                                                                                                                    >  1:816527/1‑150 (MQ=255)
                    ttttATCACGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATTAACTTGCAACATACGAGCGGTGATCCACAGATTGTTGCTCATCTCTTATTTGATTTGcc                                                                                                                                 <  2:1467738/150‑1 (MQ=255)
                    ttttATCACGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGcc                                                                                                                                 <  2:1525368/150‑1 (MQ=255)
                    ttttATCACGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGcc                                                                                                                                 <  2:647648/150‑1 (MQ=255)
                     tttATCACGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGccc                                                                                                                                <  1:744659/150‑1 (MQ=255)
                         tCACGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTa                                                                                                                            >  1:200953/1‑150 (MQ=255)
                           aCGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTAtt                                                                                                                          <  1:1531532/150‑1 (MQ=255)
                           aCGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTAtt                                                                                                                          >  2:161050/1‑150 (MQ=255)
                            cGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTg                                                                                                                         <  1:1452673/150‑1 (MQ=255)
                                ctGCGCCGTTCATGGCTTTCATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTTATCCACAGATGGGTGCCCATCTCTTGTTTGATTTGCCCTTTATTGCCTa                                                                                                                     <  2:283045/149‑1 (MQ=255)
                                   cgcCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAAcc                                                                                                                  <  2:847871/150‑1 (MQ=255)
                                          cATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCaaa                                                                                                           <  2:1261385/150‑1 (MQ=255)
                                             ggCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAAcc                                                                                                        >  1:616514/1‑150 (MQ=255)
                                               cTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAg                                                                                                      >  2:261365/1‑150 (MQ=255)
                                                  tGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGCCCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCg                                                                                                   >  2:687072/1‑150 (MQ=255)
                                                   gATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGATCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCgg                                                                                                  <  1:722487/150‑1 (MQ=255)
                                                   gATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCgg                                                                                                  <  1:976196/150‑1 (MQ=255)
                                                    aTACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGc                                                                                                 <  2:1384597/150‑1 (MQ=255)
                                                    aTACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGc                                                                                                 <  2:1390263/150‑1 (MQ=255)
                                                     tACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCAATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCa                                                                                                <  2:133081/150‑1 (MQ=255)
                                                      aCCGTTATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCAc                                                                                               >  2:197735/1‑150 (MQ=255)
                                                      aCCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCAc                                                                                               <  1:1099645/150‑1 (MQ=255)
                                                       ccGAGATCCACCAACGAGTAAGAGCCAGGTCGACCCATCGCCGCAGCAACGGAAT‑AACGTGCCACAAACGAGCGGTGATCCACAGATCGGTGACCATCTCTTCTTTGATTTGACCTTTATTGCCTAACCAGCCAAAACAAGTCGGCACTa                                                                                             >  2:1338134/1‑150 (MQ=255)
                                                              caaccaacGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTACCCAGCCAAAACCAGTCGGCACTACGGAAt                                                                                       >  1:1049692/1‑150 (MQ=255)
                                                              caaccaacGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAAt                                                                                       >  2:1333590/1‑150 (MQ=255)
                                                                 ccaacGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAAtttt                                                                                    >  2:1056437/1‑150 (MQ=255)
                                                                 ccaacGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAAtttt                                                                                    >  2:663098/1‑150 (MQ=255)
                                                                 ccaacGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAAtttt                                                                                    >  1:1299044/1‑150 (MQ=255)
                                                                   aacGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTa                                                                                  <  2:301704/150‑1 (MQ=255)
                                                                     cGAGTAAGCGCCAGGTCGCCCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCCCTACGGAATTTTTAcc                                                                                >  2:906102/1‑150 (MQ=255)
                                                                       aGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCaa                                                                              <  2:652883/150‑1 (MQ=255)
                                                                        gTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCaaa                                                                             >  1:1560564/1‑150 (MQ=255)
                                                                           aGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATc                                                                          <  2:200953/150‑1 (MQ=255)
                                                                             cgcCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCaa                                                                        >  1:758106/1‑150 (MQ=255)
                                                                              gcCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCACTTTATTGCCTAACCAGCCAAAACCAGTCGGCGCTACGGAATTTTTACAAAAATCaaa                                                                       >  2:553497/1‑150 (MQ=255)
                                                                               ccAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAg                                                                      <  2:1115491/150‑1 (MQ=255)
                                                                               ccAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAg                                                                      <  2:580619/150‑1 (MQ=255)
                                                                                    tCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGcg                                                                 <  2:274525/150‑1 (MQ=255)
                                                                                     cGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGcgg                                                                >  2:640016/1‑150 (MQ=255)
                                                                                       aCCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGcggtc                                                              >  2:573794/1‑150 (MQ=255)
                                                                                          cATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGcggtcggt                                                           >  1:1609467/1‑150 (MQ=255)
                                                                                           aTCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGtt                                                          <  2:816527/150‑1 (MQ=255)
                                                                                           aTCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGtt                                                          >  2:1547982/1‑150 (MQ=255)
                                                                                              gccgcAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGTTTcc                                                       >  2:973981/1‑150 (MQ=255)
                                                                                                cgcAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCATTTATTGCATAACCAGCAAAAACCAGTCGTCACTACGGAATTTTTACCAAAATAAAAGATGAGGTCGGTTTCCTg                                                     >  1:402630/1‑150 (MQ=255)
                                                                                                  cagcaACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACAAAAATCAAAGATGCGGTCGGTTTCCTGtt                                                   <  2:127620/150‑1 (MQ=255)
                                                                                                    gcaACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGTTTCCTGTTcc                                                 <  1:395248/150‑1 (MQ=255)
                                                                                                    gcaACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGTTTCCTGTTcc                                                 <  2:123757/150‑1 (MQ=255)
                                                                                                        cGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGTTTCCTGTTCCAGcc                                             <  1:1333590/150‑1 (MQ=255)
                                                                                                         ggAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGTTTCCTGTTCCAGCCa                                            >  1:988531/1‑150 (MQ=255)
                                                                                                          gAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGTTTCCTGTTCCAGCCaa                                           >  2:928298/1‑150 (MQ=255)
                                                                                                           aaTAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGTTTCCTGTTCCAGCCAAc                                          <  1:153507/150‑1 (MQ=255)
                                                                                                             tAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGTTTCCTGTTCCAGCCAACgg                                        >  1:1536026/1‑150 (MQ=255)
                                                                                                             tAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGTTTCCTGTTCCAGCCAACgg                                        >  1:1635923/1‑150 (MQ=255)
                                                                                                                 cGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGTTTCCTGTTCCAGCCAACGGTtgt                                    >  1:430421/1‑150 (MQ=255)
                                                                                                                   tGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTCATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGTTTCCTGTTCCAGCCAACGGTTGTgg                                  >  2:95840/1‑150 (MQ=255)
                                                                                                                    gCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGTTTCCTGTTCCAGCCAACGGTTGTGGc                                 >  2:180810/1‑150 (MQ=255)
                                                                                                                        cATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGTTTCCTGTTCCAGCCAACGGTTGTGGCTTAg                             <  2:758106/150‑1 (MQ=255)
                                                                                                                          tACGAGCGGTGATCCACTGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGGTTCCTGTTCCAGCCAACGGTTGTGGCTTAggg                           <  1:591431/150‑1 (MQ=255)
                                                                                                                          tACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGTTTCCTGTTCCAGCCAACGGTTGTGGCTTAggg                           <  2:224305/150‑1 (MQ=255)
                                                                                                                           aCGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGTTTCCTGTTCCAGCCAACGGTTGTGGCTTAGGgt                          <  1:596859/150‑1 (MQ=255)
                                                                                                                               gCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGTTTCCTGTTCCAGCCAACGGTTGTGGCTTAGGGTGTTa                      <  1:261365/150‑1 (MQ=255)
                                                                                                                                cGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGTTTCCTGTTCCAGCCAACGGTTGTGGCTTAGGGTGTTaa                     <  2:513022/150‑1 (MQ=255)
                                                                                                                                  gTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGTTTCCTGTTCCAGCCAACGGTTGTGGCTTAGGGTGTTAAAc                   >  2:1091085/1‑150 (MQ=255)
                                                                                                                                  gTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTAACCAAAATCAAAGATGCGGTCGGTTTCCTGTTCCAGCCAACGGTTGTGGCTTAGGGTGTTAAAc                   <  2:973674/150‑1 (MQ=255)
                                                                                                                                   tGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGTTTCCTGTTCCAGCCAACGGTTGTGGCTTAGGGTGTTAAAcc                  <  2:1609467/150‑1 (MQ=255)
                                                                                                                                    gATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGTTTCCTGTTCCAGCCAACGGTTGTGGCTTAGGGTGTTAAACCa                 >  1:1136279/1‑150 (MQ=255)
                                                                                                                                     aagaagaGATGGGCACCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGTTTCCTGTTCCAGCCAACGGTTGTGGCTTAGGGTGTTAAACCAt                >  1:152760/7‑150 (MQ=255)
                                                                                                                                     aagaagaGATGGGCACCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGTTTCCTGTTCCAGCCAACGGTTGTGGCTTAGGGTGTTAAACCAt                >  2:403653/7‑150 (MQ=255)
                                                                                                                                      tCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGTTTCCTGTTCCAGCCAACGGTTGTGGCTTAGGGTGTTAAACCAtt               >  1:451511/1‑150 (MQ=255)
                                                                                                                                      agaagaGATGGGCACCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGTTTCCTGTTCCAGCCAACGGTTGTGGCTTAGGGTGTTAAACCAtt               >  2:99475/6‑150 (MQ=255)
                                                                                                                                       gaagaGATGGGCACCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGTTTCCTGTTCCAGCCAACGGTTGTGGCTTAGGGTGTTAAACCAttt              >  1:1443396/5‑150 (MQ=255)
                                                                                                                                       ccACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGTTTCCTGTTCCAGCCAACGGTTGTGGCTTAGGGTGTTAAACCAttt              >  1:924251/1‑150 (MQ=255)
                                                                                                                                       ccACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGTTTCCTGTTCCAGCCAACGGTTGTGGCTTAGGGTGTTAAACCAttt              >  1:539197/1‑150 (MQ=255)
                                                                                                                                        cacaGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGTTTCCTGTTCCAGCCAACGGTTGTGGCTTAGGGTGTTAAACCATTTc             <  1:197735/150‑1 (MQ=255)
                                                                                                                                          caGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGTTTCCTGTTCCAGCCAACGGTTGTGGCTTAGGGTGTTAAACCATTTCAt           <  2:988531/150‑1 (MQ=255)
                                                                                                                                            gATGGGCACCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGTTTCCTGTTCCAGCCAACGGTTGTGGCTTAGGGTGTTAAACCATTTCATTc         <  2:1443396/150‑1 (MQ=255)
                                                                                                                                                ggTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGTTTCCTGTTCCAGCCAACGGTTGTGGCTTAGGGTGTTAAACCATTTCATTCCttt     >  2:630315/1‑150 (MQ=255)
                                                                                                                                                ggTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGTTTCCTGTTCCAGCCAACGGTTGTGGCTTAGGGTGTTAAACCATTTCATTCCttt     <  2:1461903/150‑1 (MQ=255)
                                                                                                                                                ggTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGTTTCCTGTTCCAGCCAACGGTTGTGGCTTAGGGTGTTAAACCATTTCATTCCttt     <  1:1091085/150‑1 (MQ=255)
                                                                                                                                                 gTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGTTTCCTGTTCCAGCCAACGGTTGTGGCTTAGGGTGTTAAACCATTTCATTCCtttt    <  1:417799/150‑1 (MQ=255)
                                                                                                                                                  tGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGTTTCCTGTTCCAGCCAACGGTTGTGGCTTAGGGTGTTAAACCATTTCATTCCTTTTa   <  1:1056437/150‑1 (MQ=255)
                                                                                                                                                   gCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGTTTCCTGTTCCAGCCAACGGTTGTGGCTTAGGGTGTTAAACCATTTCATTCCTTTTAt  >  1:1530529/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GGCATACCAGCCGCCATATTTTTTATCACGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGTTTCCTGTTCCAGCCAACGGTTGTGGCTTAGGGTGTTAAACCATTTCATTCCTTTTAT  >  NC_000913/4071469‑4071765

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: ‑