Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,966,876 A→C 8.1% T154P (ACT→CCT)  rho → transcription termination factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,966,8760AC 8.1% 317.7 / 13.9 109T154P (ACT→CCT) rhotranscription termination factor
Reads supporting (aligned to +/- strand):  ref base A (44/56);  new base C (4/5);  total (48/61)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01

Only 100 of 111 total aligned reads displayed.

CGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTG  >  NC_000913/3966727‑3967025
                                                                                                                                                     |                                                                                                                                                     
cgAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTa                                                                                                                                                       >  2:1544675/1‑150 (MQ=255)
 gAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTAc                                                                                                                                                      >  1:945953/1‑150 (MQ=255)
 gAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTAc                                                                                                                                                      <  1:1140710/150‑1 (MQ=255)
 gAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTAc                                                                                                                                                      >  1:1322637/1‑150 (MQ=255)
 gAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTAc                                                                                                                                                      <  1:756130/150‑1 (MQ=255)
    agaagGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGa                                                                                                                                                   <  1:836822/150‑1 (MQ=255)
    agaagGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACACAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGa                                                                                                                                                   >  1:208083/1‑150 (MQ=255)
     gaagGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTCCTGaa                                                                                                                                                  <  1:379639/150‑1 (MQ=255)
     gaagGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGaa                                                                                                                                                  >  1:1402390/1‑150 (MQ=255)
        ggTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGAt                                                                                                                                               >  1:201500/1‑150 (MQ=255)
        ggTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGAt                                                                                                                                               <  1:962923/150‑1 (MQ=255)
          tGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGAttt                                                                                                                                             <  1:56402/150‑1 (MQ=255)
           gAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTCCTGAAGATTTa                                                                                                                                            >  1:552141/1‑150 (MQ=255)
             aCGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAAc                                                                                                                                          <  1:317048/150‑1 (MQ=255)
                cTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGc                                                                                                                                       >  1:566892/1‑150 (MQ=255)
                 tATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCt                                                                                                                                      <  1:1332501/150‑1 (MQ=255)
                   ttttGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTcg                                                                                                                                    <  2:1026385/150‑1 (MQ=255)
                     ttGCGCTGCTGAAAGTTAACGAAGATAACTTCGACAAACCTGAACACGCCCGCAACAATATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTATCGGTTCTCCTGAAGATTTAACTGCTcgcg                                                                                                                                  <  1:1196551/150‑1 (MQ=255)
                        cgcTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATGCTCTTTGAGAACTTAACCCCGCTGCACGCAAGCTCTCGGCTGCGTATGGAACGTGGTAACGGTGCTACTGAAGATTTAACTGCTCGCGTag                                                                                                                               >  1:282024/1‑149 (MQ=255)
                        cgcTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTATGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGAGTATGGAACTTGTTATCTGATCTCCTGAAGATTTAACAGCTCGCGTaa                                                                                                                               >  1:278142/1‑149 (MQ=255)
                           tgctgAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTgg                                                                                                                            <  2:1402390/150‑1 (MQ=255)
                            tctgAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCTCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGa                                                                                                                           <  1:1568745/149‑1 (MQ=255)
                             ctgAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGAt                                                                                                                          <  1:587469/150‑1 (MQ=255)
                                 aaGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTATCGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTgg                                                                                                                      >  1:1436280/1‑150 (MQ=255)
                                 aaGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTgg                                                                                                                      >  1:583012/1‑150 (MQ=255)
                                    ttAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTCACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCAt                                                                                                                   <  1:1067136/150‑1 (MQ=255)
                                     tAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATc                                                                                                                  <  2:444478/150‑1 (MQ=255)
                                       aCGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAATCTCTCCTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCAc                                                                                                                >  2:218568/1‑150 (MQ=255)
                                        cGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCAcc                                                                                                               <  1:1296799/150‑1 (MQ=255)
                                        cGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCAcc                                                                                                               <  2:899792/150‑1 (MQ=255)
                                            gTTAACTTCGACAAACCTGAAAACTCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATc                                                                                                           <  1:1025435/150‑1 (MQ=255)
                                            gTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATc                                                                                                           >  1:223358/1‑150 (MQ=255)
                                                aCTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTc                                                                                                       >  1:476153/1‑150 (MQ=255)
                                                 cTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCg                                                                                                      <  2:948778/150‑1 (MQ=255)
                                                 cTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCg                                                                                                      <  2:1157882/150‑1 (MQ=255)
                                                  ttCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGt                                                                                                     <  1:45230/150‑1 (MQ=255)
                                                   tCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACTTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTg                                                                                                    >  2:1749702/1‑150 (MQ=255)
                                                     gACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGt                                                                                                  >  2:17198/1‑150 (MQ=255)
                                                       cAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCa                                                                                                <  2:833310/150‑1 (MQ=255)
                                                        aaaCCTGAACACGCCCGCAACAAAATGCTCTTTGAGAACTTAACCCCGCTGCACGCAACCTCTCGGCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCTTCACCTATCGGTCGTGGTCAg                                                                                               <  2:1852505/150‑1 (MQ=255)
                                                         aaCCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGc                                                                                              >  1:1594070/1‑150 (MQ=255)
                                                         aaCCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACAGCACGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGc                                                                                              >  2:1905086/1‑150 (MQ=255)
                                                          aCCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCg                                                                                             <  2:826000/150‑1 (MQ=255)
                                                            cTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTg                                                                                           >  2:1081527/1‑150 (MQ=255)
                                                             tGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCAGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTgg                                                                                          <  2:668082/150‑1 (MQ=255)
                                                             tGAAAACGCCCGCAACAAAATCCACATTGAGAACTTCACCCCGCTGCACGCAAACTCTCTTCTGCGAATGGAACGTGGTAACGGTTCTACTGAAGATTAAACTGCGCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTgg                                                                                          <  2:279347/150‑1 (MQ=255)
                                                                aaaCGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCt                                                                                       <  1:1380343/150‑1 (MQ=255)
                                                                aaaCGCAAGCAACAAAATCCTCTTTGAGAACTTAAACCCGCTGCAAGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGATGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGGGGTCt                                                                                       >  2:47657/1‑150 (MQ=255)
                                                                   cGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGAt                                                                                    <  1:845766/150‑1 (MQ=255)
                                                                    gCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGttt                                                                                   >  2:818040/1‑150 (MQ=255)
                                                                     cccGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATtg                                                                                  >  2:500938/1‑150 (MQ=255)
                                                                       cGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTTGTCTGCGTATGGAACGTGGTAACGGTTCTCCTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATtgtg                                                                                <  1:976400/150‑1 (MQ=255)
                                                                       cGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTCCTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATtgtg                                                                                >  2:1903339/1‑150 (MQ=255)
                                                                       cGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATtgtg                                                                                <  1:1684735/150‑1 (MQ=255)
                                                                        gCAACAAAATCCTCTTTGAGCACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTCGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACAGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTgg                                                                               >  1:677455/1‑150 (MQ=255)
                                                                         caacaaAATCCTCTTTGAGACCTTAACCCCGCTGCACGCTAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGGGGTCTGATTGTGGc                                                                              <  2:1239299/150‑1 (MQ=255)
                                                                          aacaaAATACTCTTTGAGAACTTAACCCCGCTGCACACAAACTCTCGTCTGCGTATGGAACGAGGTAACGGTGCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCa                                                                             <  1:81308/150‑1 (MQ=255)
                                                                             aaaaTCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCAccg                                                                          <  1:1624036/150‑1 (MQ=255)
                                                                               aaTCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCAccgcc                                                                        >  1:1387957/1‑150 (MQ=255)
                                                                               aaTCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCAccgcc                                                                        <  1:1544675/150‑1 (MQ=255)
                                                                                 tCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGa                                                                      <  1:616232/150‑1 (MQ=255)
                                                                                    tctTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAg                                                                   >  1:1751719/1‑150 (MQ=255)
                                                                                      tttGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGcc                                                                 >  2:3678/1‑150 (MQ=255)
                                                                                      tttGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGcc                                                                 >  1:1617165/1‑150 (MQ=255)
                                                                                      tttGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGcc                                                                 >  1:1903457/1‑150 (MQ=255)
                                                                                       ttGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCg                                                                <  2:86337/150‑1 (MQ=255)
                                                                                        tGAGAACTTTACCCCGCTGCACGCAAACTCTCGGCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCgg                                                               <  2:187671/150‑1 (MQ=255)
                                                                                        tGAGAACTTACCCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGGCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCgg                                                               <  1:1778409/150‑1 (MQ=255)
                                                                                        tGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCgg                                                               >  2:643434/1‑150 (MQ=255)
                                                                                          agaACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTa                                                             <  2:235963/150‑1 (MQ=255)
                                                                                          agaACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTa                                                             <  2:241563/150‑1 (MQ=255)
                                                                                                 aaCCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCAt                                                      >  2:895109/1‑150 (MQ=255)
                                                                                                  aCCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCAtg                                                     <  2:1751719/150‑1 (MQ=255)
                                                                                                      cGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCAtgctgc                                                 <  2:1482224/150‑1 (MQ=255)
                                                                                                         tGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCAtgctgctgc                                              <  1:1275178/150‑1 (MQ=255)
                                                                                                          gCACGCAAACTCTCGGCTGCGTATGGAACGTGGTAACGGTGCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCa                                             <  2:1129201/150‑1 (MQ=255)
                                                                                                           cACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAg                                            <  2:201500/150‑1 (MQ=255)
                                                                                                           cACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGa                                           <  2:1559085/150‑1 (MQ=255)
                                                                                                           cACGCAAACTCTAGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTCCTCGCGTACTGGATATGGCATCACCTATCGGTCGTGATCAGAGTGGTCTGATAGTGGCAACGCCGAAAGCCGGAAAAACAATGCTTCTGCAg                                            >  1:1338642/1‑150 (MQ=255)
                                                                                                                   ctctCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGc                                    >  2:400367/1‑150 (MQ=255)
                                                                                                                   ctctCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATATAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCCCCGAAATCCGGTAAAACCATTCTGCTGCAGAACATTGc                                    >  1:901370/1‑150 (MQ=255)
                                                                                                                       cGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCag                                >  1:1112273/1‑150 (MQ=255)
                                                                                                                         tCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCagag                              <  1:1626554/150‑1 (MQ=255)
                                                                                                                          cTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGc                             <  2:877050/150‑1 (MQ=255)
                                                                                                                            gCGTATGGAACGTGGTAAAGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGAAGAACATTGCTCAGAGCAt                           >  1:1394593/1‑150 (MQ=255)
                                                                                                                              gTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTg                         <  1:146910/150‑1 (MQ=255)
                                                                                                                                aTGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCt                       <  1:560008/150‑1 (MQ=255)
                                                                                                                                  ggAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGACCATTGCTCAGAGCATTGCTTa                     >  2:1577887/1‑150 (MQ=255)
                                                                                                                                  ggAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTa                     >  1:872366/1‑150 (MQ=255)
                                                                                                                                   gAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCTGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTAc                    <  2:669218/150‑1 (MQ=255)
                                                                                                                                    aaCGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAACGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACa                   >  2:1160600/1‑150 (MQ=255)
                                                                                                                                       gttgTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTTTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAacc                <  2:1137544/147‑1 (MQ=255)
                                                                                                                                       gTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAacc                <  2:1322637/150‑1 (MQ=255)
                                                                                                                                        cggTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAacca               >  2:1251729/2‑150 (MQ=255)
                                                                                                                                          gTAACGGTTCTCCTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGGAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAaccacc             >  2:961850/1‑150 (MQ=255)
                                                                                                                                            aaCGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCg           >  2:510369/1‑150 (MQ=255)
                                                                                                                                               ggTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGGCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGAt        <  1:1411993/150‑1 (MQ=255)
                                                                                                                                                gTTCTCCTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGAtt       <  2:552141/150‑1 (MQ=255)
                                                                                                                                                gTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTTCTTACAACCACCCGGAtt       >  1:1431285/1‑150 (MQ=255)
                                                                                                                                                 ttCTCCTGAAGATTTAACTGCTCGCATACTTGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATtg      <  1:1100332/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTG  >  NC_000913/3966727‑3967025

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