Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,966,876 A→C 5.2% 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 5.2% 905.9 / 32.1 264T154P (ACT→CCT) rhotranscription termination factor
Reads supporting (aligned to +/- strand):  ref base A (116/127);  new base C (10/4);  total (129/135)
Fisher's exact test for biased strand distribution p-value = 1.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 267 total aligned reads displayed.

CGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTG  >  NC_000913/3966727‑3967025
                                                                                                                                                     |                                                                                                                                                     
cgAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTa                                                                                                                                                       >  1:1438984/1‑150 (MQ=255)
 gAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTAc                                                                                                                                                      <  2:877589/150‑1 (MQ=255)
 gAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTAc                                                                                                                                                      <  2:2682482/150‑1 (MQ=255)
 gAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTAc                                                                                                                                                      <  1:1701910/150‑1 (MQ=255)
 gAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTAc                                                                                                                                                      <  1:1697463/150‑1 (MQ=255)
  aaaGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACt                                                                                                                                                     >  2:1522442/1‑150 (MQ=255)
     gaagGTGACCGCTATTTTGCGCTGCTGAATGTTAACGACGTTAACTTCGACAAACCTGAACACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGaa                                                                                                                                                  <  1:2014008/150‑1 (MQ=255)
     gaagGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGACAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGaa                                                                                                                                                  <  2:2413401/150‑1 (MQ=255)
            aaCGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAATGGTTCTACTGAAGATTTaa                                                                                                                                           <  2:2744177/150‑1 (MQ=255)
              cGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACt                                                                                                                                         <  2:3709187/150‑1 (MQ=255)
               gatATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCATCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTg                                                                                                                                        >  2:3803236/3‑150 (MQ=255)
               gCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTg                                                                                                                                        >  2:449305/1‑150 (MQ=255)
                cTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGCTTTAACTGc                                                                                                                                       <  2:1579697/150‑1 (MQ=255)
                  aTTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTc                                                                                                                                     <  2:3428775/150‑1 (MQ=255)
                               gAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGACCGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCt                                                                                                                        >  2:592583/1‑150 (MQ=255)
                               gAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCt                                                                                                                        <  1:1579650/150‑1 (MQ=255)
                               gAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCt                                                                                                                        >  2:3354682/1‑150 (MQ=255)
                               gAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCt                                                                                                                        >  2:1954448/1‑150 (MQ=255)
                                  aGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGc                                                                                                                     >  1:2532501/1‑150 (MQ=255)
                                   gTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTCCTGAAGATTTAACTGCTAGCGTACTGGATCTGGCa                                                                                                                    >  1:2702254/1‑150 (MQ=255)
                                   gTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCa                                                                                                                    >  1:3331851/1‑150 (MQ=255)
                                   gTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCa                                                                                                                    <  2:2217634/150‑1 (MQ=255)
                                    ttAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTCCTGAAGATTTAACTGCTCGCGTACTGGATCTGGCAt                                                                                                                   >  1:1320838/1‑150 (MQ=255)
                                    ttAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCAt                                                                                                                   >  1:1366752/1‑150 (MQ=255)
                                    ttAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCAt                                                                                                                   >  1:3271923/1‑150 (MQ=255)
                                      aaCGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCa                                                                                                                 <  2:2381832/150‑1 (MQ=255)
                                      aaCGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCa                                                                                                                 <  2:1092880/150‑1 (MQ=255)
                                      aaCGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCa                                                                                                                 <  1:14648/150‑1 (MQ=255)
                                         gAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTAGTCTGCGTATGGAACGTGGTAACGGTTCTAATGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCt                                                                                                              >  1:2872682/1‑150 (MQ=255)
                                            gTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATc                                                                                                           <  2:3137739/150‑1 (MQ=255)
                                               aaCTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGt                                                                                                        <  1:3566250/150‑1 (MQ=255)
                                                aCTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTc                                                                                                       <  1:955028/150‑1 (MQ=255)
                                                  ttCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGt                                                                                                     >  1:4106673/1‑150 (MQ=255)
                                                     gACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTATCTGCTCGCGTACTGGATCTGGCATCTCCTATCGGTCGTGGt                                                                                                  >  2:2013965/1‑150 (MQ=255)
                                                        aaaCCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAg                                                                                               <  1:3141010/150‑1 (MQ=255)
                                                        aaaCCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAg                                                                                               <  1:1076259/150‑1 (MQ=255)
                                                         aaCCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGc                                                                                              <  1:1625506/150‑1 (MQ=255)
                                                           ccTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGCTTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGt                                                                                            <  2:3741648/150‑1 (MQ=255)
                                                             tGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTGCTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTgg                                                                                          <  1:1533665/150‑1 (MQ=255)
                                                                   cGCCCGCACCAAAATCCTCTTTGAGAACTGAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGAt                                                                                    <  2:1608305/150‑1 (MQ=255)
                                                                    gCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGAtt                                                                                   <  1:702801/150‑1 (MQ=255)
                                                                     cccGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATtg                                                                                  <  1:3021420/150‑1 (MQ=255)
                                                                      ccGCAACATAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATtgt                                                                                 <  2:3322988/150‑1 (MQ=255)
                                                                      ccGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATtgt                                                                                 >  2:198424/1‑150 (MQ=255)
                                                                      ccGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATtgt                                                                                 >  1:3784294/1‑150 (MQ=255)
                                                                       cGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATtgtg                                                                                <  1:3951247/150‑1 (MQ=255)
                                                                       cGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATtgtg                                                                                <  2:385547/150‑1 (MQ=255)
                                                                        gCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGCTCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTgg                                                                               <  2:2327217/150‑1 (MQ=255)
                                                                            caaAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCAcc                                                                           >  2:1416251/1‑150 (MQ=255)
                                                                              aaaTCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTCCTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCAccgc                                                                         <  2:1995705/150‑1 (MQ=255)
                                                                              aaaTCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCAccgc                                                                         >  2:280298/1‑150 (MQ=255)
                                                                               aaTCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCAccgcc                                                                        >  1:4042663/1‑150 (MQ=255)
                                                                               aaTCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCAccgcc                                                                        <  1:811797/150‑1 (MQ=255)
                                                                                aTCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCAccgccg                                                                       <  2:3637878/150‑1 (MQ=255)
                                                                                aTCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCAccgccg                                                                       >  1:3978675/1‑150 (MQ=255)
                                                                                aTCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCAccgccg                                                                       >  1:3521771/1‑150 (MQ=255)
                                                                                    tctTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTCCTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTTGCACCGCCGAAAg                                                                   >  2:2830112/1‑150 (MQ=255)
                                                                                      tttGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGcc                                                                 <  2:2357001/150‑1 (MQ=255)
                                                                                       ttGAGAACTTAACCCCGCTTCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGAGTGTGGCACCGCCGAAAGCCg                                                                <  1:891226/150‑1 (MQ=255)
                                                                                         gagaACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGt                                                              >  2:1186174/1‑150 (MQ=255)
                                                                                           gaACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTGCTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTaa                                                            >  1:4079304/1‑150 (MQ=255)
                                                                                           gaACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTGCTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTaa                                                            >  1:2044321/1‑150 (MQ=255)
                                                                                             aCTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTaaaa                                                          <  2:2459441/150‑1 (MQ=255)
                                                                                             aCTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTaaaa                                                          <  2:1634627/150‑1 (MQ=255)
                                                                                             aCTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTaaaa                                                          <  2:3010547/150‑1 (MQ=255)
                                                                                               ttAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAAcc                                                        >  1:21821/1‑150 (MQ=255)
                                                                                                    cccGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGACTGATTGTGGCACCGCCGAAAGCCGGTAAAACCAtgct                                                   >  1:2143122/1‑150 (MQ=255)
                                                                                                     ccGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTCCTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCAtgctg                                                  >  1:2551123/1‑150 (MQ=255)
                                                                                                               cAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACa                                        <  2:834788/150‑1 (MQ=255)
                                                                                                               cAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACa                                        <  2:1450059/150‑1 (MQ=255)
                                                                                                               cAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACa                                        <  1:739499/150‑1 (MQ=255)
                                                                                                                 aaCTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACAtt                                      >  1:2969445/1‑150 (MQ=255)
                                                                                                                  aCTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTg                                     >  2:1946173/1‑150 (MQ=255)
                                                                                                                   ctctCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGc                                    <  1:1153398/150‑1 (MQ=255)
                                                                                                                    tctcGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCt                                   <  1:3716011/150‑1 (MQ=255)
                                                                                                                    tctcGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCt                                   >  1:3838328/1‑150 (MQ=255)
                                                                                                                        gTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCaga                               <  2:1929580/150‑1 (MQ=255)
                                                                                                                        gTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCaga                               <  1:3006770/150‑1 (MQ=255)
                                                                                                                              gTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTg                         <  1:15760/150‑1 (MQ=255)
                                                                                                                                aTGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCt                       >  1:2207559/1‑150 (MQ=255)
                                                                                                                                 tGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCtt                      >  1:281/1‑150 (MQ=255)
                                                                                                                                 tGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCtt                      >  1:1942734/1‑150 (MQ=255)
                                                                                                                                   gAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTAc                    <  1:3438479/150‑1 (MQ=255)
                                                                                                                                      cGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAac                 >  2:1788423/1‑150 (MQ=255)
                                                                                                                                      cGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAac                 >  2:1166854/1‑150 (MQ=255)
                                                                                                                                        tGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAacca               >  1:1289244/1‑150 (MQ=255)
                                                                                                                                        tGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAacca               <  2:2716190/150‑1 (MQ=255)
                                                                                                                                         ggTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAaccac              <  2:1166370/150‑1 (MQ=255)
                                                                                                                                           tAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCAccc            >  2:354042/1‑150 (MQ=255)
                                                                                                                                           tAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCAccc            >  1:1246798/1‑150 (MQ=255)
                                                                                                                                                gTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGAtt       >  2:3172427/1‑150 (MQ=255)
                                                                                                                                                 ttCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATtg      >  2:1427670/1‑150 (MQ=255)
                                                                                                                                                 ttCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATtg      >  2:2519883/1‑150 (MQ=255)
                                                                                                                                                  tCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATtgt     <  2:3114448/150‑1 (MQ=255)
                                                                                                                                                   cTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGGATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATtgtg    <  2:3581929/150‑1 (MQ=255)
                                                                                                                                                   cTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATtgtg    >  1:3380781/1‑150 (MQ=255)
                                                                                                                                                    tACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATtgtgt   <  1:449305/150‑1 (MQ=255)
                                                                                                                                                    tACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATtgtgt   <  2:3873881/150‑1 (MQ=255)
                                                                                                                                                    tACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATtgtgt   >  2:4062136/1‑150 (MQ=255)
                                                                                                                                                     aCTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATtgtgtg  >  1:2332665/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTG  >  NC_000913/3966727‑3967025

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