Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,997,991 G→T 5.1% A163A (GCC→GCA) ‡ tcyL ← cystine ABC transporter membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,997,9910GT 5.1% 336.8 / 4.8 98A163A (GCC→GCA) ‡tcyLcystine ABC transporter membrane subunit
Reads supporting (aligned to +/- strand):  ref base G (48/45);  new base T (2/3);  total (50/48)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.60e-01

Only 100 of 103 total aligned reads displayed.

TGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGTAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATCAA  >  NC_000913/1997842‑1998139
                                                                                                                                                     |                                                                                                                                                    
tGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGTAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTcgcg                                                                                                                                                      <  1:469011/150‑1 (MQ=255)
 gATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGTAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCgg                                                                                                                                                     <  1:1639916/150‑1 (MQ=255)
   ttCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGTAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGcc                                                                                                                                                   >  2:205049/1‑150 (MQ=255)
    tctcAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGTAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCa                                                                                                                                                  <  1:1966477/150‑1 (MQ=255)
     ctcAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGTAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAg                                                                                                                                                 >  1:1693982/1‑150 (MQ=255)
      tcAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGTAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGc                                                                                                                                                >  2:1238559/1‑150 (MQ=255)
       cAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGTAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCg                                                                                                                                               <  2:957160/150‑1 (MQ=255)
       cAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGTAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCg                                                                                                                                               >  2:1374872/1‑150 (MQ=255)
       cAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGTAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCg                                                                                                                                               >  1:120271/1‑150 (MQ=255)
         aaaTGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGTAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCTGGCCGCGAc                                                                                                                                             <  2:162943/150‑1 (MQ=255)
         aaaTGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGTAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCTGGCCGCGAc                                                                                                                                             <  2:1423551/150‑1 (MQ=255)
         aaaTGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGTAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGAc                                                                                                                                             <  2:1146696/150‑1 (MQ=255)
          aaTGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGTAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACg                                                                                                                                            <  2:534335/150‑1 (MQ=255)
          aaTGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGTAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACg                                                                                                                                            <  2:1887688/150‑1 (MQ=255)
               ttCTGCAACGTCGATAACACAGTCGCCATGATCCAGTAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCt                                                                                                                                       <  2:1826766/150‑1 (MQ=255)
                tCTGCAACGTCGATAACACAGTCGCCATGATCCAGTAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCTGGCCGCGAcaatccag                                                                                                                                      <  2:1138493/150‑8 (MQ=255)
                tCTGCAACGTCGATAACACAGTCGCCATGATCCAGTAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCtt                                                                                                                                      <  2:711960/150‑1 (MQ=255)
                tCTGCAACGTCGATAACACAGTCGCCATGATCCAGTAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCtt                                                                                                                                      >  2:1205462/1‑150 (MQ=255)
                 cTGCAACGTCGATAACACAGTCGCCATGATCCAGTAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCttt                                                                                                                                     <  1:1238559/150‑1 (MQ=255)
                   gCAACGTCGATAACACAGTCGCCATGATCCAGTAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTa                                                                                                                                   >  1:45517/1‑150 (MQ=255)
                   gCAACGTCGATAACACAGTCGCCATGATCCAGTAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTa                                                                                                                                   >  2:1025645/1‑150 (MQ=255)
                    cAACGTCGATAACACAGTCGCCATGATCCAGTAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTAc                                                                                                                                  >  1:1436281/1‑150 (MQ=255)
                    cAACGTCGATAACACAGTCGCCATGATCCAGTAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTAc                                                                                                                                  >  1:439104/1‑150 (MQ=255)
                      aCGTCGATAACACAGTCGCCATGATCCAGTAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCa                                                                                                                                >  2:169260/1‑150 (MQ=255)
                       cGTCGATAACACAGTCGCCATGATCCAGTAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAg                                                                                                                               <  2:1405935/150‑1 (MQ=255)
                                cacaGTCGCCATGATCCAGTAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATaaa                                                                                                                      >  1:918761/1‑150 (MQ=255)
                                  caGTCGCCATGATCCAGTAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGc                                                                                                                    <  2:1173234/150‑1 (MQ=255)
                                    gTCGCCATGATCCAGTAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTg                                                                                                                  <  1:663987/150‑1 (MQ=255)
                                       gCCATGATCCAGTAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGTGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTc                                                                                                               <  2:58879/150‑1 (MQ=255)
                                       gCCATGATCCAGTAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTc                                                                                                               <  1:362955/150‑1 (MQ=255)
                                       gCCATGATCCAGTAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTc                                                                                                               <  1:370023/150‑1 (MQ=255)
                                          aTGATCCAGTAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGAc                                                                                                            >  2:1635242/1‑150 (MQ=255)
                                           tGATCCAGTGGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTc                                                                                                               <  2:1094162/146‑1 (MQ=255)
                                           tGATCCAGTGGATGAGAGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGtta                                                                                                               >  1:1094162/1‑145 (MQ=255)
                                           tGATCCAGTAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACa                                                                                                           >  2:1137653/1‑150 (MQ=255)
                                           tGATCCAGTAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACa                                                                                                           >  1:464110/1‑150 (MQ=255)
                                                 aGTAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAgcggcg                                                                                                     >  1:439076/1‑150 (MQ=255)
                                                  gTAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGcggcgg                                                                                                    >  2:1107575/1‑150 (MQ=255)
                                                   tAGATGAGCGAAGCCTCCAGATACATGGTTAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTTACGTAACTGCTCCGGAACCTGGTTTGTCGAGGCCAGAGACGTATCTTTTACCAGGCTAATAAAGCTTTTCGACAGCggcggc                                                                                                   <  2:1013098/150‑1 (MQ=255)
                                                   tAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCggcggc                                                                                                   <  2:1356779/150‑1 (MQ=255)
                                                    aGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCa                                                                                                  >  2:21493/1‑150 (MQ=255)
                                                       tGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGACCAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGTAACg                                                                                               <  1:205049/150‑1 (MQ=255)
                                                       tGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCggcggcgcgt                                                                                               >  2:1462239/1‑146 (MQ=255)
                                                         aGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGcc                                                                                             >  1:990598/1‑150 (MQ=255)
                                                         aGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGcc                                                                                             <  2:120271/150‑1 (MQ=255)
                                                          gCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCca                                                                                            >  2:1080344/1‑150 (MQ=255)
                                                           cGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCcac                                                                                           >  2:1038059/1‑150 (MQ=255)
                                                               gccgccAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACAcgc                                                                                       >  2:666963/1‑150 (MQ=255)
                                                                ccgccAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACAcgcg                                                                                      <  1:1107575/150‑1 (MQ=255)
                                                                ccgccAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACAcgcg                                                                                      >  1:859422/1‑150 (MQ=255)
                                                                 cgccAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACgcgc                                                                                     >  1:1132939/1‑150 (MQ=255)
                                                                  gccAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGcgcc                                                                                    <  1:1720565/150‑1 (MQ=255)
                                                                     aGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCgccgcc                                                                                 <  1:1558907/150‑1 (MQ=255)
                                                                        tACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCtgt                                                                              <  2:487916/150‑1 (MQ=255)
                                                                          cATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTgg                                                                            <  1:1374872/150‑1 (MQ=255)
                                                                             ggTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCaa                                                                         >  2:503641/1‑150 (MQ=255)
                                                                             ggTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCaa                                                                         >  2:159053/1‑150 (MQ=255)
                                                                             ggTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCaa                                                                         >  1:779141/1‑150 (MQ=255)
                                                                               tGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCaaaa                                                                       <  2:830783/150‑1 (MQ=255)
                                                                                gAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAAt                                                                      >  2:290588/1‑150 (MQ=255)
                                                                                gAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAAt                                                                      >  1:1755329/1‑150 (MQ=255)
                                                                                gAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAAt                                                                      >  2:935204/1‑150 (MQ=255)
                                                                                 aaCACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACAGGCGGCGCGTGTGGcgttgcc                                                                     <  1:996694/150‑7 (MQ=255)
                                                                                   cacCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATcgc                                                                   >  1:695584/1‑150 (MQ=255)
                                                                                     ccTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCACCGCCTGTGGCAAAATcgcgc                                                                 >  1:1957303/1‑150 (MQ=255)
                                                                                      cTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATcgcgcg                                                                <  2:45517/150‑1 (MQ=255)
                                                                                         cAGCGTACGCGAGGTAATCAAGTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGTCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCAGACGCGCCGCCTGTGGCAAAATCGCGCGACg                                                             <  2:1816985/150‑1 (MQ=255)
                                                                                          aGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGc                                                            >  1:1342688/1‑150 (MQ=255)
                                                                                           gCGTAAGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCa                                                           >  2:853992/1‑150 (MQ=255)
                                                                                            cGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCAt                                                          <  2:918761/150‑1 (MQ=255)
                                                                                                 gcgAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCTCGCGACGCATGGTCt                                                     <  2:757253/150‑1 (MQ=255)
                                                                                                   gAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGc                                                   <  2:1034171/150‑1 (MQ=255)
                                                                                                       tAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCCCg                                               >  1:1416595/1‑150 (MQ=255)
                                                                                                           cAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTATTAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGt                                           >  1:471336/1‑150 (MQ=255)
                                                                                                            aaCTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTc                                          <  2:1972085/150‑1 (MQ=255)
                                                                                                             aCTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCa                                         >  1:1451491/1‑150 (MQ=255)
                                                                                                                gcgcCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCACCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATAc                                      <  2:1957303/150‑1 (MQ=255)
                                                                                                                        cGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTg                              <  2:302292/150‑1 (MQ=255)
                                                                                                                          gAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACCCGCGCCGCCTGTGGCAAAATCGCGCGCCGCATGGTCTGCCACGGCGTCATACCAATACTGGc                            >  2:1156941/1‑150 (MQ=255)
                                                                                                                          gAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATAcc                                     <  2:1158203/141‑1 (MQ=255)
                                                                                                                          gAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATAcc                                     >  1:1158203/1‑141 (MQ=255)
                                                                                                                           aaCAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCt                           >  2:1253154/1‑150 (MQ=255)
                                                                                                                             cAGCTCCGGCAGCTTGATTGTCGCGGGCTTCGACGTATCTTTTACCAGGCTATTATTGCTGTTCGACAGCGGCGGCAAAGCCACACGCGCCGCCTGTGGCAAAATCGGGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGc                         <  2:1696712/150‑1 (MQ=255)
                                                                                                                             cAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGc                         >  2:1424907/1‑150 (MQ=255)
                                                                                                                              aGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGcc                        <  1:290588/150‑1 (MQ=255)
                                                                                                                                 tCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCt                     <  2:439104/150‑1 (MQ=255)
                                                                                                                                 tCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCt                     >  2:1911251/1‑150 (MQ=255)
                                                                                                                                  ccGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCtt                    <  2:439076/150‑1 (MQ=255)
                                                                                                                                   cGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTc                   <  2:1755329/150‑1 (MQ=255)
                                                                                                                                    ggCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTcc                  <  2:695584/150‑1 (MQ=255)
                                                                                                                                     gCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTCCCAGGCTAATACAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTccc                 <  1:1232993/150‑1 (MQ=255)
                                                                                                                                     gCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTccc                 <  1:935204/150‑1 (MQ=255)
                                                                                                                                      cACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCa                >  2:1019767/1‑150 (MQ=255)
                                                                                                                                         aagatacgTCGCTGGCCGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTg             >  1:325447/8‑150 (MQ=255)
                                                                                                                                          tGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTATTAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGa            <  2:471336/150‑1 (MQ=255)
                                                                                                                                           ggATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGAc           >  1:1631956/1‑150 (MQ=255)
                                                                                                                                            atacgTCGCTGGCCGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGAcc          >  1:542861/5‑150 (MQ=255)
                                                                                                                                                  cgcgGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATc    <  1:1173443/150‑1 (MQ=255)
                                                                                                                                                  cgcgGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATc    <  1:1038059/150‑1 (MQ=255)
                                                                                                                                                    cgGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATCaa  <  1:70011/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGTAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATCAA  >  NC_000913/1997842‑1998139

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: