Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,679,639 A→C 6.0% Q28P (CAA→CCA)  ydgI → putative arginine:ornithine antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,679,6390AC 6.0% 205.8 / 5.2 67Q28P (CAA→CCA) ydgIputative arginine:ornithine antiporter
Reads supporting (aligned to +/- strand):  ref base A (35/28);  new base C (2/2);  total (37/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.26e-01

CCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTACCTATGCCCGCGAA  >  NC_000913/1679494‑1679787
                                                                                                                                                 |                                                                                                                                                    
ccATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAat                                                                                                                                                  <  1:755858/150‑1 (MQ=255)
        aTGCGTTTTGCCCTAATATTTATTCTTCTTATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATggcggc                                                                                                                                          <  1:620855/150‑1 (MQ=255)
        aTGCGTTTTGCCCTAATATTTATTCTTCTTATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATggcggc                                                                                                                                          >  1:523882/1‑150 (MQ=255)
          gCGTTTTGCCCTAATATTTATTCTTCTTATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAg                                                                                                                                        >  1:759014/1‑150 (MQ=255)
          gCGTTTTGCCCTAATATTTATTCTTCTTATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAg                                                                                                                                        <  1:727763/150‑1 (MQ=255)
           cGTTTTGCCCTAATATTTATTCTTCTTATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGt                                                                                                                                       <  2:1372405/150‑1 (MQ=255)
            gTTTTGCCCTAATATTTATTCTTCTTATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCCATATTTTGCGGCAgac                                                                                                                                      <  1:1163359/150‑3 (MQ=255)
            gTTTTGCCCTAATATTTATTCTTCTTATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCCATATTTTGCGGCAgac                                                                                                                                      <  1:892741/150‑3 (MQ=255)
              tttGCCCTAATATTTATTCTTCTTATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGc                                                                                                                                    <  1:187225/150‑1 (MQ=255)
                    cTAATATTTATTCTTCTTATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGCGACTGAGCGCACTCAACGCGCTGGTATTAAGCTCAATGCTGGGAGCGGCTGTTTTCAGTCTGACGCAAAATATGGATGCAGTTGCCAGccc                                                                                                                              >  2:638186/1‑150 (MQ=255)
                       atatTTATTCTTCTTATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCGCAATGCTGGGCGCGGGTGGTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGgc                                                                                                                           <  2:477526/150‑1 (MQ=255)
                         atttattCTTCTTATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGAAAAATATGGCGGCAGATGACAGACCGgcag                                                                                                                         >  2:878489/1‑150 (MQ=255)
                               tcttctTATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGc                                                                                                                   <  2:1037636/150‑1 (MQ=255)
                                cttcttATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCt                                                                                                                  <  2:759014/150‑1 (MQ=255)
                                       tCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGc                                                                                                           <  1:1216593/150‑1 (MQ=255)
                                        cACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCt                                                                                                          >  2:1192653/1‑150 (MQ=255)
                                         aCGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTg                                                                                                         >  2:1332374/1‑150 (MQ=255)
                                          cGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTgg                                                                                                        >  2:1176289/1‑150 (MQ=255)
                                            ttttAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTgggg                                                                                                      <  1:1332374/150‑1 (MQ=255)
                                                   cACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACt                                                                                               >  2:1089788/1‑150 (MQ=255)
                                                    aCTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTg                                                                                              >  1:72549/1‑150 (MQ=255)
                                                       ggATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCgg                                                                                                                            >  2:658402/1‑117 (MQ=255)
                                                       ggATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCgg                                                                                                                            <  1:658402/117‑1 (MQ=255)
                                                       ggATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGgcg                                                                                           >  1:240876/1‑150 (MQ=255)
                                                       ggATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGgcg                                                                                           >  1:245961/1‑150 (MQ=255)
                                                          tatCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTg                                                                                        <  1:489422/150‑1 (MQ=255)
                                                             cGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCa                                                                                     >  2:1124703/1‑150 (MQ=255)
                                                              gATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCAt                                                                                    >  2:348571/1‑150 (MQ=255)
                                                                   aaaaGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTAt                                                                               <  1:184916/150‑1 (MQ=255)
                                                                        aaaCTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTg                                                                          <  1:878489/150‑1 (MQ=255)
                                                                         aaCTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTgg                                                                         >  2:1337841/1‑150 (MQ=255)
                                                                            tGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCt                                                                      >  1:214569/1‑150 (MQ=255)
                                                                                    aGCGCACTCACCGCGCGGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATg                                                              <  2:674552/150‑1 (MQ=255)
                                                                                     gcgcACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGc                                                             <  2:72549/150‑1 (MQ=255)
                                                                                        cACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGa                                                          >  2:1325085/1‑150 (MQ=255)
                                                                                        cACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGa                                                          <  1:1192653/150‑1 (MQ=255)
                                                                                        cACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGa                                                          >  2:497348/1‑150 (MQ=255)
                                                                                        cACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGa                                                          >  2:803481/1‑150 (MQ=255)
                                                                                          cTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATc                                                        >  2:310523/1‑150 (MQ=255)
                                                                                           tCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATcc                                                       <  1:310523/150‑1 (MQ=255)
                                                                                            cACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCt                                                      >  1:203289/1‑150 (MQ=255)
                                                                                                     ggTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCAt                                             >  1:1002011/1‑150 (MQ=255)
                                                                                                     ggTATTAAGCTCAATGATGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTAAACGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCATTTGCCATGCTCATCCTCACGCGCAt                                             >  2:3872/1‑150 (MQ=255)
                                                                                                       tATTAAGCTCAATGCTGGCCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTc                                           >  2:1292556/1‑150 (MQ=255)
                                                                                                           aaGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTcc                                       >  1:123556/1‑150 (MQ=255)
                                                                                                            aGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCa                                                    >  1:701156/1‑136 (MQ=255)
                                                                                                            aGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCa                                                    <  2:701156/136‑1 (MQ=255)
                                                                                                              cTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGa                                    <  1:638186/150‑1 (MQ=255)
                                                                                                                  aTGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCg                                                                                                             <  2:849100/73‑1 (MQ=255)
                                                                                                                  aTGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCg                                                                                                             >  1:849100/1‑73 (MQ=255)
                                                                                                                     cTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGAc                             >  1:879593/1‑150 (MQ=255)
                                                                                                                      tGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGAcg                            >  2:1203179/1‑150 (MQ=255)
                                                                                                                       gggCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATc                                                                                                              <  1:622467/67‑1 (MQ=255)
                                                                                                                       gggCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATc                                                                                                              >  2:622467/1‑67 (MQ=255)
                                                                                                                          cgcgGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGAcggcgg                        <  1:1089788/150‑1 (MQ=255)
                                                                                                                          cgcgGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGAcggcgg                        <  1:348571/150‑1 (MQ=255)
                                                                                                                          cgcgGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGAcggcgg                        <  1:1124703/150‑1 (MQ=255)
                                                                                                                              ggTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATc                    >  2:1121905/1‑150 (MQ=255)
                                                                                                                                  ttttCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTa                <  2:879593/150‑1 (MQ=255)
                                                                                                                                    ttCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGGCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTAcc              <  1:1325085/150‑1 (MQ=255)
                                                                                                                                    ttCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTAcc              <  1:803481/150‑1 (MQ=255)
                                                                                                                                       caacTGCCGCCATATTTTGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTACCTAt           >  1:82202/4‑150 (MQ=255)
                                                                                                                                       caacTGCCGCCATATTTTGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTACCTAt           >  1:661882/4‑150 (MQ=255)
                                                                                                                                            gCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTACCTATGCCcg      >  1:523310/1‑150 (MQ=255)
                                                                                                                                              cGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTACCTATGCCcgcg    >  2:402796/1‑150 (MQ=255)
                                                                                                                                               gCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTACCTATGCCCGCGa   >  1:86467/1‑150 (MQ=255)
                                                                                                                                                cAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTACCTATGCCCGCGaa  >  2:952273/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTACCTATGCCCGCGAA  >  NC_000913/1679494‑1679787

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