Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,997,993 C→G 5.5% A163P (GCC→CCC) ‡ 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,9930CG 5.5% 348.1 / 6.8 102A163P (GCC→CCC) ‡tcyLcystine ABC transporter membrane subunit
Reads supporting (aligned to +/- strand):  ref base C (48/48);  new base G (2/4);  total (50/52)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 104 total aligned reads displayed.

TTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGTAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATCAATAG  >  NC_000913/1997845‑1998142
                                                                                                                                                    |                                                                                                                                                     
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:1887688/150‑1 (MQ=255)
       aaTGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGTAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACg                                                                                                                                               <  2:534335/150‑1 (MQ=255)
             tCTGCAACGTCGATAACACAGTCGCCATGATCCAGTAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCTGGCCGCGAcaatccag                                                                                                                                         <  2:1138493/150‑8 (MQ=255)
             tCTGCAACGTCGATAACACAGTCGCCATGATCCAGTAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCtt                                                                                                                                         >  2:1205462/1‑150 (MQ=255)
             tCTGCAACGTCGATAACACAGTCGCCATGATCCAGTAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCtt                                                                                                                                         <  2:711960/150‑1 (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:439104/1‑150 (MQ=255)
                 cAACGTCGATAACACAGTCGCCATGATCCAGTAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTAc                                                                                                                                     >  1:1436281/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)
                                        tGATCCAGTGGATGAGAGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGtta                                                                                                                  >  1:1094162/1‑145 (MQ=255)
                                        tGATCCAGTAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACa                                                                                                              >  1:464110/1‑150 (MQ=255)
                                        tGATCCAGTAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACa                                                                                                              >  2:1137653/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)
                                                    tGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACg                                                                                                  >  2:1558907/1‑150 (MQ=255)
                                                      aGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGcc                                                                                                <  2:120271/150‑1 (MQ=255)
                                                      aGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGcc                                                                                                >  1:990598/1‑150 (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:159053/1‑150 (MQ=255)
                                                                          ggTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCaa                                                                            >  2:503641/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                                        >  1:1158203/1‑141 (MQ=255)
                                                                                                                       gAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATAcc                                        <  2:1158203/141‑1 (MQ=255)
                                                                                                                        aaCAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCt                              >  2:1253154/1‑150 (MQ=255)
                                                                                                                          cAGCTCCGGCAGCTTGATTGTCGCGGGCTTCGACGTATCTTTTACCAGGCTATTATTGCTGTTCGACAGCGGCGGCAAAGCCACACGCGCCGCCTGTGGCAAAATCGGGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGc                            <  2:1696712/150‑1 (MQ=255)
                                                                                                                          cAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGc                            >  2:1424907/1‑150 (MQ=255)
                                                                                                                              tCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCt                        >  2:1911251/1‑150 (MQ=255)
                                                                                                                              tCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCt                        <  2:439104/150‑1 (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)
                                                                                                                                             gTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGCTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTa         <  1:624360/150‑1 (MQ=255)
                                                                                                                                               cgcgGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATc       <  1:1173443/150‑1 (MQ=255)
                                                                                                                                               cgcgGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATc       <  1:1038059/150‑1 (MQ=255)
                                                                                                                                                 cgGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATCaa     <  1:70011/150‑1 (MQ=255)
                                                                                                                                                   gCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATCAATa   <  1:1635242/150‑1 (MQ=255)
                                                                                                                                                   gCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATCAATa   >  1:140682/1‑150 (MQ=255)
                                                                                                                                                    ccAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTACCGCTTCCCACTGACCTTTATCAATag  <  1:1205462/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGTAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATCAATAG  >  NC_000913/1997845‑1998142

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: