Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,679,641 A→T 5.9% N29Y (AAT→TAT)  ydgI → putative arginine:ornithine antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,679,6410AT 5.9% 211.4 / 2.1 67N29Y (AAT→TAT) ydgIputative arginine:ornithine antiporter
Reads supporting (aligned to +/- strand):  ref base A (35/28);  new base T (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 = 1.00e+00

CCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTACCTATGCCCGCGAAGGG  >  NC_000913/1679494‑1679790
                                                                                                                                                   |                                                                                                                                                     
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:892741/150‑3 (MQ=255)
            gTTTTGCCCTAATATTTATTCTTCTTATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCCATATTTTGCGGCAgac                                                                                                                                         <  1:1163359/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                                                                                                                               <  1:658402/117‑1 (MQ=255)
                                                       ggATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCgg                                                                                                                               >  2:658402/1‑117 (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                                                             <  1:1192653/150‑1 (MQ=255)
                                                                                        cACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGa                                                             >  2:1325085/1‑150 (MQ=255)
                                                                                        cACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGa                                                             >  2:803481/1‑150 (MQ=255)
                                                                                        cACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGa                                                             >  2:497348/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                                                                                                                 >  2:622467/1‑67 (MQ=255)
                                                                                                                       gggCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATc                                                                                                                 <  1:622467/67‑1 (MQ=255)
                                                                                                                          cgcgGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGAcggcgg                           <  1:348571/150‑1 (MQ=255)
                                                                                                                          cgcgGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGAcggcgg                           <  1:1124703/150‑1 (MQ=255)
                                                                                                                          cgcgGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGAcggcgg                           <  1:1089788/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)
                                                                                                                                                   aaTATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTACCTATGCCCGCGAAggg  <  1:1292556/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTACCTATGCCCGCGAAGGG  >  NC_000913/1679494‑1679790

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: