Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,065,426 A→T 100% intergenic (+134/‑159) yccE → / → agp uncharacterized protein/glucose‑1‑phosphatase/inositol phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,065,4260AT97.8% 288.4 / ‑3.3 88intergenic (+134/‑159)yccE/agpuncharacterized protein/glucose‑1‑phosphatase/inositol phosphatase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base T (36/49);  total (37/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TCACACTTTTTTATTGGTAAATATTGTCTCTGTATTGGTAACGCCGCAGATATTCTGTTTAGCCACAGGTGCAATTATCAGCGGCGTACGCGAGGCAGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTGCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGT  >  NC_000913/1065330‑1065520
                                                                                                |                                                                                              
tcACACTTTTTTATTGGTAAATATTGTCTCTGTATTGGTAACGCCGCAGATATTCTGTTTAGCCACAGGTGCAATTATCAGCGGCGTACGCGAGGCTggg                                                                                             <  2:587335/100‑1 (MQ=255)
 cacaCTTTTTTATTGGTAAATATTGTCTCTGTATTGGTAACGCCGCAGATATTCTGTTTAGCCACAGGTGCAATTATCAGCGGCGTACGCGAGGCTgggg                                                                                            <  1:1367313/100‑1 (MQ=255)
 cacaCTTTTTTATTGGTAAATATTGTCTCTGTATTGGCAACGCCGCAGATATTCTGTTTAGCCACAGGTGCAATTATCAGCGGCGTACGCGAGGCTgggg                                                                                            >  2:1373211/1‑100 (MQ=255)
  acacTTTTTTATTGGTAAATATTGTCTCTGTATTGGTAACGCCGCAGATATTCTGTTTAGCCACAGGTGCAATTATCAGCGGCGTACGCGAGGCTGGGGc                                                                                           <  1:1277405/100‑1 (MQ=255)
   cacTTTTTTATTGGTAAATATTGTCTCTGTATTGGTAACGCCGCAGATATTCTGTTTAGCCACAGGTGCAATTATCAGCGGCGTACGCGAGGCAGGGGCt                                                                                          >  1:1524911/1‑100 (MQ=255)
    acTTTTTTATTGGTAAATATTGTCTCTGTATTGGTAACGCCGCAGATATTCTGTTTAGCCACAGGTGCAATTATCAGCGCCGTACGCGAGGCTGGGGCTa                                                                                         >  2:741186/1‑100 (MQ=255)
       tttttATTGGTAAATATTGTCTCTGTATTGGTAACGCCGCAGATATTCTGTTTAGCCACAGGTGCAATTATCAGCGGCGTACGCGAGGCTGGGGCTAATc                                                                                      >  1:1864861/1‑100 (MQ=255)
          ttATTGGTAAATATTGTCTCTGTATTGGTAACGCCGCAGATATTCTGTTTAGCCACAGGTGCAATTATCAGCGGCGTACGCGAGGCTGGGGCTAATCAgg                                                                                   <  1:311352/100‑1 (MQ=255)
            aTTGGTAAATATTGTCTCTGTATTGGTAACGCCGCAGATATTCTGTTTAGCCACAGGTGCAATTATCAGCGGCGTACGCGAGGCTGGGGCTAATCAGGCa                                                                                 <  1:1856471/100‑1 (MQ=255)
              tGGTAAATATTGTCTCTGTATTGGTAACGCCGCAGATATTCTGTTTAGCCACAGGTGCAATTATCAGCGGCGTACGCGAGGCTGGGGCTAATCAGGCATa                                                                               <  1:1338364/100‑1 (MQ=255)
               ggTAAATATTGTCTCTGTATTGGTAACGCCGCAGATATTCTGTTTAGCCACAGGTGCAATTATCAGCGGCGTACGCGAGGCTGGGGCTAATCAGGCATAg                                                                              >  2:1268007/1‑100 (MQ=255)
               ggTAAATATTGTCTCTGTATTGGTAACGCCGCAGATATTCTGTTTAGCCACAGGTGCAATTATCAGCGGCGTACGCGAGGCTGGGGCTAATCAGGCATAg                                                                              >  2:465915/1‑100 (MQ=255)
               ggTAAATATTGTCTCTGTATTGGTAACGCCGCAGATATTCTGTTTAGCCACAGGTGCAATTATCAGCGGCGTACGCGAGGCAGGGGCTAATCAGGCATAg                                                                              <  2:191080/100‑1 (MQ=255)
                 tAAATACTGTCACTGTATTGGTTACGCCGCAGATATTCTTTTTAGTCACAGGGGTAATTATCAGCGGCATAGTCGAGGATGGGGATAAGCATGCATAGtt                                                                            >  2:1265519/1‑100 (MQ=255)
                   aaTATTGTCTCTGTATTGGTAACGCCGCAGATATTCTGTTTAGCCACAGGTGCAATTATCAGCGGCGTACGCGAGGCTGGGGCTAATCAGGCATAGTTTg                                                                          <  1:944770/100‑1 (MQ=255)
                     tatTGTCTCTGTATTGGTAACGCCGCAGTTATTCTGTTTAGCCACAGGTGCAATTATCAGCGGCGTACGCGAGGCTGGGGCTAATAAGGCATAGTTTg                                                                          >  2:507585/1‑98 (MQ=255)
                     tatTGTCTCTGTATTGGTAACGCCGCAGATATTCTGTTTAGCCACAGGTGCAATTATCAGCGGCGTACGCGAGGCTGGGGCTAATCAGGCATAGTTTGCg                                                                        <  2:223623/100‑1 (MQ=255)
                      atTGTCTCTGTATTGGTAACGCCGCAGATATTCTGTTTAGCCACAGGTGCAATTATCAGCGGCGTACGCGAGGCTGGGGCTAATCAGGCATAGTTTGCGt                                                                       <  1:861922/100‑1 (MQ=255)
                         gTCTCTGTATTGGTAACGCCGCAGATATTCTGTTTAGCCACAGGTGCAATTATCAGCGGCGTACGCGAGGCTGGGGCTAATCAGGCATAGTTTGCGTCaa                                                                    >  1:1525932/1‑100 (MQ=255)
                            tctGTATTGGTAACGCCGCAGATATTCTGTTTAGCCACAGGTGCAATTATCAGCGGCGTACGCGAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAAcc                                                                 <  2:1681471/100‑1 (MQ=255)
                            tctGTATTGGTAACGCCGCAGATATTCTGTTTAGCCACAGGTGCAATTATCAGCGGCGTACGCGAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAAcc                                                                 <  1:1399060/100‑1 (MQ=255)
                             ctGTATTGGTAACGCCGCAGATATTCTGTTTAGCCACAGGTGCAATTATCAGCGGCGTACGCGAGGCTGGGGCTACTCAGGCATAGTTTGCGTCAAACCt                                                                >  2:1971729/1‑100 (MQ=255)
                               gTATTGGTAACGCCGCAGATATTCTGTTTAGCCACAGGTGCAATTATCAGCGGCGTACGCGAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTG                                                              <  1:1216637‑M1/100‑1 (MQ=255)
                               gTATTGGTAACGCCGCAGATATTCTGTTTAGCCACAGGTGCAATTATCAGCGGCGTACGCGAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTG                                                              <  1:222313‑M1/100‑1 (MQ=255)
                                tATTGGTAACGCCGCAGATATTCTGTTTAGCCACAGGTGCAATTATCAGCGGCGTACGCGAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTGC                                                             <  1:1743913‑M1/100‑1 (MQ=255)
                                 cttccgatctccgcagatattctgTTTAGCCACAGGTGCAATTATCAGCGGCGTACGCGAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTGCC                                                            <  1:2323841‑M1/90‑1 (MQ=255)
                                 aTTGGTAACGCCGCAGATATTCTGTTTAGCCACAGGTGCAATTATCAGCGGCGTACGCGAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTGCC                                                            >  2:1141121‑M1/1‑100 (MQ=255)
                                 aTTGGTAACGCCGCAGATATTCTGTTTAGCCACAGGTGCAATTATCAGCGGCGTACGCGAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTGCC                                                            >  1:1842480‑M1/1‑100 (MQ=255)
                                  ttGGTAACGCCGCAGATATTCTGTTTAGCCACAGGTGCAATTATCAGCGGCGTACGCGAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTGCCT                                                           <  1:1857524‑M1/100‑1 (MQ=255)
                                    ggTAACGCCGCAGATATTCTGTTTAGCCACAGGTGCAATTATCAGCGGCGTACGCGAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTGCCTGT                                                         >  2:1971515‑M1/1‑100 (MQ=255)
                                    ggTAACGCCGCAGATATTCTGTTTAGCCACAGGTGCAATTATCAGCGGCGTACGCGAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTGCCTGT                                                         >  1:1656988‑M1/1‑100 (MQ=255)
                                     gTAACGCTGCAGATATTCTGTTTAGCCACAGGTGCAATTATCAGCGGCGTACGCGAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTGCCTGTT                                                        <  2:1822843‑M1/100‑1 (MQ=255)
                                       aaCGCCGCAGATATTCTGTTTAGCCACAGGTGCAATTATCAGCGGCGTACGCGAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTGCCTGTTTT                                                      >  1:470828‑M1/1‑100 (MQ=255)
                                          gccgcAGATATTCTGTTTAGCCACAGGTGCAATTATCAGCGGCGTACGCGAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTGCCTGTTTTTGa                                                   <  2:1006806/100‑1 (MQ=255)
                                           ccgcAGATATTCTGTTTAGCCACAGGTGCAATTATCAGCGGCGTACGCGAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTGCCagatcggaag                                                  >  2:2323827‑M1/1‑90 (MQ=255)
                                           ccgcAGATATTCTGTTTAGCCACAGGTGCAATTATCAGCGGCGTACGCGAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTGCCTGTTTTTcta                                                  >  2:369279‑M1/1‑97 (MQ=255)
                                            cgcAGATATTCTGTTTAGCCACAGGTGCAATTATCAGCGGCGTACGCGAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTGCCTGTTTTTGAAg                                                 <  2:513226/100‑1 (MQ=255)
                                                gATATTCTGTTTAGCCACAGGTGCAATTATCAGCGGCGTACGCGAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTGCCTGTTTTTGAAGATGt                                             >  1:85001/1‑100 (MQ=255)
                                                 atatTCTGTTTAGCCACAGGTGCAATTATCAGCGGCGTACGCGAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTGCCTGTTTTTGAAGATGta                                            <  1:156750/100‑1 (MQ=255)
                                                 atatTCTGTTTAGCCACAGGTGCAATTATCAGCGGCGTACGCGAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTGCCTGTTTTTGAAGATGta                                            <  1:155104/100‑1 (MQ=255)
                                                 atatTCTGTTTAGCCACAGGTGCAATTATCAGCGGCGTACGCGAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTGCCTGTTTTTGAAGATGta                                            <  2:689445/100‑1 (MQ=255)
                                                  tatTCTGTTTAGCCACAGGTGCAATTATCAGCGGCGTACGCGAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTGCCTGTTTTTGAAGATGtat                                           >  1:1896301/1‑100 (MQ=255)
                                                   atTCTGTTTAGCCACAGGTGCAATTATCAGCGGCGTACGCGAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTGCCTGTTTTTGAAGATGtata                                          <  2:1745235/100‑1 (MQ=255)
                                                     tCTGTTTAGCCACAGGTGCAATTATCAGCGGCGTACGCGAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTGCCTGTTTTTGAAGATGtatata                                        >  2:2203214/1‑100 (MQ=255)
                                                     tCTGTTTAGCCACAGGTGCAATTATCAGCGGCGTACGCGAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTGCCTGTTTTTGAAGATGtatata                                        >  2:488539/1‑100 (MQ=255)
                                                     tCTGTTTAGCCACAGGTGCAATTATCAGCGGCGTACGCGAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTGCCTGTTTTTGAAGATGtatata                                        <  1:1197874/100‑1 (MQ=255)
                                                         tttAGCCACAGGTGCAATTATCAGGGGCGTACGCGAGGCTGGGGCTAATCAGGCATAGTTTGTGTCAAATCTTGCCTGTTTTTGAAGATGTATATAGa                                      <  2:1228745/98‑1 (MQ=255)
                                                           tAGCCACAGGTGCAATTATCAGCGGCGTACGCGAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTGCCTGTTTTTGAAGATGTATATAGaaaaa                                  <  2:973953‑M2/100‑6 (MQ=255)
                                                             gCCACAGGTGCAATTATCAGCGGCGTACGCGAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTGCCTGTTTTTctatatacatcttcaaaaaca                                <  1:1390039‑M1/100‑22 (MQ=255)
                                                              ccccatctGCAATTATCAGCGGCTTACGCTAGGCTGGTTTATTCAGGAATCGTTTGCTTTAAATCTTGCTTGTTTTTGAAGATGtatat                                         <  2:1575165/82‑1 (MQ=255)
                                                               cacaGGTGCAATTATCAGCGGCGTACGCGAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTGCCTGTTTTTGAAGATGTATATAGaaaaacagg                              >  1:2076760‑M2/1‑91 (MQ=255)
                                                               cacaGGTGCAATTATCAGCGGCGTACGCGAGGCTAGGGCTAATCAGTCATAGTTTGCGTCAAACCTTGCCTGTTTTTGAAGATGAATATAGaaaaatagg                              <  2:1111676‑M2/100‑10 (MQ=255)
                                                                acaGGTGCAATTATCAGCGGCGTACGCGAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTGCCTGTTTTTctatatacatcttcaaaaacaggc                             <  2:1874738‑M1/100‑25 (MQ=255)
                                                                acaGGTGCAATTATCAGCGGCGTACGCGAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTGCCTGTTTTTGAAGATGTATATAGaaaaacaggc                             >  1:1987589‑M2/1‑90 (MQ=255)
                                                                acaGGTGCAATTATCAGCGGCGTACGCGAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTGCCTGTTTTTGAAGATGTATATAGaaaaacaggc                             >  1:882167‑M2/1‑90 (MQ=255)
                                                                acaGGTGCAAGTATCAGCGGCGTAAGCGAGGCTGGGTCTAATCAGGCATAGTTTGCGTCAAACCTTGCCTGTTTTTGAAAAAAAAAAAAAAAAAACAGGc                             <  2:1622068/100‑1 (MQ=255)
                                                                   ggTGCAATTATCAGCGGCGTACGCGAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTGCCTGTTTTTctatatacatcttcaaaaacaggcaag                          <  1:322820‑M1/100‑28 (MQ=255)
                                                                   ggTGCAATTATCAGCGGCGTACGCGAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTGCCTGTTTTTctatatacatcttcaaaaacaggcaag                          <  2:159176‑M1/100‑28 (MQ=255)
                                                                   ggTGCAATTATCAGCGGCGTACGCGAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTGCCTGTTTTTctatatacatcttcaaaaacaggcaag                          <  1:369283‑M1/100‑28 (MQ=255)
                                                                     tGCAATTATCAGCGGCGTACGCGAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTGCCTGTTTTTctatatacatcttcaaaaacaggcaaggt                        <  2:1313144‑M1/100‑30 (MQ=255)
                                                                     tGCAATTATCAGCGGCGTACGCGAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTGCCTGTTTTTctatatacatcttcaaaaacaggcaaggt                        <  1:1764339‑M1/100‑30 (MQ=255)
                                                                     tGCAATTATCAGCGGCGTACGCGAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTGCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTca                        <  2:1245832/100‑1 (MQ=255)
                                                                      gCAATTATCAGCGGCGTACGCGAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTGCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTcaa                       <  1:1594587/100‑1 (MQ=255)
                                                                       cAATTATCAGCGGCGTACGCGAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTGCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTcaac                      <  1:205783/100‑1 (MQ=255)
                                                                       cAATTATCAGCGGCGTACGCGAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTGCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTcaac                      >  2:471814/1‑100 (MQ=255)
                                                                        aaTTATCAGCGGCGTACGCGAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTGCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTcaaca                     >  2:1886449/1‑100 (MQ=255)
                                                                        aaTTATCAGCGGCGTACGCGAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTGCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTcaaca                     <  1:175443/100‑1 (MQ=255)
                                                                        aaTTATCAGCGGCGTACGCGAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTGCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTcaaca                     <  1:862988/100‑1 (MQ=255)
                                                                         aTTATCAGCGGCGTACGCGAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTGCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTcaacaa                    <  1:1765693/100‑1 (MQ=255)
                                                                         aTTATCAGCGGCGTACGCGAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTGCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTcaacaa                    <  1:967099/100‑1 (MQ=255)
                                                                            aTCAGCGGCGTACGCGAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTGCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTCAACAAGcc                 <  2:1962956/100‑1 (MQ=255)
                                                                            aTCAGCGGCGTACGCGAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTGCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTCAACAAGcc                 <  2:2000393/100‑1 (MQ=255)
                                                                              cAGCGGCGTACGCGAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTGCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTCAACAAGCCAt               >  1:827499/1‑100 (MQ=255)
                                                                              cAGCGGCGTACGCGAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTGCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTCAACAAGCCAt               >  1:900610/1‑100 (MQ=255)
                                                                                gcggcgTACGCGAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTGCCTGTTTTTctatatacatcttcaaaaacaggcaaggtttgacgcaaac             <  1:19385‑M1/100‑41 (MQ=255)
                                                                                  ggcgTACGCGAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTGCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTg           <  1:949935/100‑1 (MQ=255)
                                                                                  ggcgTACGCGAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTGCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTg           >  2:2534175/1‑100 (MQ=255)
                                                                                  ggcgTACGCGAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTGCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTg           >  1:2021480/1‑100 (MQ=255)
                                                                                  ggcgTACGCGAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTGCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTg           >  1:613097/1‑100 (MQ=255)
                                                                                   gcgTACGCGAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTGCCTGTTTTTTTTTTTGTTTATAGAAAAACAGGCGTTCAACAAGCCATTTTGc          >  2:1649794/1‑100 (MQ=255)
                                                                                        cgctAGGCTGGGGCTAATCAGGCATAGTTGGCGTCAAACCTGTCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAAcc     <  1:2137937/100‑1 (MQ=255)
                                                                                        cgcgAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTGCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAAcc     >  2:398149/1‑100 (MQ=255)
                                                                                        cgcgAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTGCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAAcc     <  1:1143713/100‑1 (MQ=255)
                                                                                         gcgAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTGCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCt    >  1:358006/1‑100 (MQ=255)
                                                                                         gcgAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTGCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCt    >  1:400695/1‑100 (MQ=255)
                                                                                          cgAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTGCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTg   >  2:1883171/1‑100 (MQ=255)
                                                                                          cgAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTGCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTg   >  1:1392746/1‑100 (MQ=255)
                                                                                           gAGGCTGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTGCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGt  <  1:1283254/100‑1 (MQ=255)
                                                                                                |                                                                                              
TCACACTTTTTTATTGGTAAATATTGTCTCTGTATTGGTAACGCCGCAGATATTCTGTTTAGCCACAGGTGCAATTATCAGCGGCGTACGCGAGGCAGGGGCTAATCAGGCATAGTTTGCGTCAAACCTTGCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGT  >  NC_000913/1065330‑1065520

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