Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,966,770 A→C 12.4% E118D (GAA→GAC) ‡ rho → transcription termination factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,966,7700AC12.4% 314.7 / 31.8 121E118D (GAA→GAC) ‡rhotranscription termination factor
Reads supporting (aligned to +/- strand):  ref base A (57/49);  new base C (4/11);  total (61/60)
Fisher's exact test for biased strand distribution p-value = 5.79e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.96e-01

Only 100 of 123 total aligned reads displayed.

AGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTAT  >  NC_000913/3966624‑3966919
                                                                                                                                                  |                                                                                                                                                     
cgCTCCTCCCGCGCCGGTCCTGATGACATCTACGATTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCCCTCGTGATACCATCTCTGGTAACATTCGCCCGCCGAAAGAAGGAGACCGCAATTTTGCGCTGCTGAAAGTTAACGAAGtt                                                                                                                                                    <  2:116577/149‑1 (MQ=255)
aGCTCCTACCTCGCCGGTCCAGATGACATCTACGGTTCCCCTATCCAAATCCGCCGTTTCAACCTCCGAACTGGTGATACTATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGACGtt                                                                                                                                                    <  1:1144798/150‑1 (MQ=255)
 gCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTAGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGACGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAACTTa                                                                                                                                                   >  2:68885/1‑150 (MQ=255)
      tACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTc                                                                                                                                              <  2:1462772/150‑1 (MQ=255)
       aCCTCGCCGGTCCTGATGACATCTACGTTTCACCTAGCCAAATCCGCCGATTCAACCTCCGCACTGGTGATACCATCTCTGGGAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGGAAGTTAACTTCGTTAACtttc                                                                                                                                             <  2:104490/150‑3 (MQ=255)
        ccTCTCCGTTCCTGATTACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCAGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCAGAAAGAAGGTTAACGCTATTTTGCGCTGCTGAAAGTTAACTTCGTTAACtttca                                                                                                                                            <  1:1151555/150‑4 (MQ=255)
        ccTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCgc                                                                                                                                            >  2:1439582/1‑149 (MQ=255)
        ccTCGCCGGTCCTGATGAAATCTACGTTTCACCTAGACAAATCAGCCGTTTCAACCTCCGCACTCTTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGa                                                                                                                                            <  1:330582/150‑1 (MQ=255)
         cTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTTAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACTTCGTTAACtttcag                                                                                                                                           <  1:604126/150‑5 (MQ=255)
         cTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACTTCGTTAACtttcag                                                                                                                                           <  1:188335/150‑5 (MQ=255)
         cTCGCCGGTACTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACTTCGTTAACtttcag                                                                                                                                           <  1:141913/150‑5 (MQ=255)
          tCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACTTCGTTAACtttcagc                                                                                                                                          <  2:43999/150‑6 (MQ=255)
          tCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACTTCGTTAACtttcagc                                                                                                                                          <  2:375502/150‑6 (MQ=255)
          tCGCCGGGCCTGAACACAACTACGATTCCCCAAGCCAAATCCGCCGTTTCAACATCCGCTCTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCAGCTGAAAGTTAACTTCGTTAACtttcagc                                                                                                                                          <  1:350655/150‑6 (MQ=255)
                  ccTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGaaa                                                                                                                                 <  1:1376115/150‑1 (MQ=255)
                    tgatgaCATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGaaaa                                                                                                                                >  1:1034688/1‑150 (MQ=255)
                    tgatgaCATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGAAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGaaaa                                                                                                                                <  2:821845/150‑1 (MQ=255)
                            tCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGca                                                                                                                        >  1:1355960/1‑150 (MQ=255)
                            tCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGca                                                                                                                        >  1:832702/1‑150 (MQ=255)
                             cTACGTTTCCCCTAGCCAAATCCGCCGATTCAACCTCCGCACTGGTGAAACCATCTCTGGTAAGAGTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTACGACAAACCTGAAAACGCCCGcaa                                                                                                                       >  1:1156677/1‑150 (MQ=255)
                              tACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGAAACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGCTGAACGCTATATTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAACCGCCCGcaac                                                                                                                      >  2:1553408/1‑150 (MQ=255)
                                   ttCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTTGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAAt                                                                                                                 <  1:675157/150‑1 (MQ=255)
                                    tCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCTCTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACTAACCTGAAAACGCCCGCAACAAAATc                                                                                                                <  2:864823/150‑1 (MQ=255)
                                    tCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATc                                                                                                                <  1:961202/150‑1 (MQ=255)
                                     ccccTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAACCGCACGCAACAAACTcc                                                                                                               >  1:267268/1‑150 (MQ=255)
                                     ccacTAGCCAAATCCGCCGATTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGAGAACGCTATTTTGCGCTGCTGAAAGTTACCGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATcc                                                                                                               >  2:1291272/4‑150 (MQ=255)
                                       ccTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTCACTTCGACAAACCTGAAAACGCCCGCAACAAAATCctc                                                                                                             >  2:71236/1‑150 (MQ=255)
                                       ccTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCctc                                                                                                             <  1:38827/150‑1 (MQ=255)
                                          aGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCttt                                                                                                          <  1:181190/150‑1 (MQ=255)
                                          aGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCttt                                                                                                          <  1:40646/150‑1 (MQ=255)
                                           gCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTg                                                                                                         >  1:645867/1‑150 (MQ=255)
                                            ccAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCACCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTgt                                                                                                        >  1:921380/1‑149 (MQ=255)
                                                aTCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAAc                                                                                                    <  2:818392/150‑1 (MQ=255)
                                                   cgccgTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCtt                                                                                                           <  1:1189501/140‑1 (MQ=255)
                                                   cgccgTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCtt                                                                                                           >  2:1189501/1‑140 (MQ=255)
                                                   cgccgTTTCAACCTCAGCACTGGGTATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTa                                                                                                 <  1:294455/150‑1 (MQ=255)
                                                       gTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAAccc                                                                                             >  1:1334435/1‑150 (MQ=255)
                                                          tCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGc                                                                                          >  2:1721830/1‑150 (MQ=255)
                                                           cAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCt                                                                                         >  1:390240/1‑150 (MQ=255)
                                                            aaCCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTg                                                                                        <  2:1416593/150‑1 (MQ=255)
                                                              ccTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCa                                                                                      >  2:649455/1‑150 (MQ=255)
                                                              ccTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCa                                                                                      >  1:945327/1‑150 (MQ=255)
                                                                 ccGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAACATCCTCTTTGAGAACTTAACCCCGCTGCACGc                                                                                   >  1:85568/1‑150 (MQ=255)
                                                                  cGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCa                                                                                  >  1:636478/1‑150 (MQ=255)
                                                                        ggTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTACCCCCGCTGCACGCAAActct                                                                            >  1:1274518/1‑150 (MQ=255)
                                                                        ggTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAActct                                                                            >  1:1144747/1‑150 (MQ=255)
                                                                        ggTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGAGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAActct                                                                            >  2:1193044/1‑150 (MQ=255)
                                                                         gTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACtctc                                                                           <  2:938474/150‑1 (MQ=255)
                                                                              aCCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTg                                                                      <  2:1320164/150‑1 (MQ=255)
                                                                                cATCTCTGGTATTATTCGCGCGGCGAAAGAAGGTTATCGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCg                                                                    <  2:916280/150‑1 (MQ=255)
                                                                                cATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCg                                                                    <  2:651962/150‑1 (MQ=255)
                                                                                cATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCg                                                                    <  2:446158/150‑1 (MQ=255)
                                                                                 aTCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGt                                                                   <  1:1482478/150‑1 (MQ=255)
                                                                                    tctGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATg                                                                <  2:1382699/150‑1 (MQ=255)
                                                                                     ctGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATgg                                                               >  2:1038960/1‑150 (MQ=255)
                                                                                         tAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGACCTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACg                                                           >  2:1711520/1‑150 (MQ=255)
                                                                                            gATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCTACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTCCACGCAAACTCTCGTCTGCGTATGGAACGTgg                                                        >  2:323982/1‑150 (MQ=255)
                                                                                              ttCGCCCGCAGAAAGAAGTTGAACGCTATGTTGCGCTTCTGAAAGATAACGAAGTTAACTTCGACAAACCTGAAATCGCCCGCAACTAAATCCTCTTTGAGAACTAAACCCCGCTGCACGCAACCTCTCTTATGCCTATGGAACGTGGTa                                                      <  2:121127/150‑1 (MQ=255)
                                                                                                cgcccgccGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAAc                                                    <  2:358296/150‑1 (MQ=255)
                                                                                                 gcccgccGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACg                                                   >  1:828289/1‑150 (MQ=255)
                                                                                                  cccgccGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACgg                                                  <  2:1102739/150‑1 (MQ=255)
                                                                                                     gccgAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTc                                               <  2:320314/150‑1 (MQ=255)
                                                                                                        gAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTAc                                            <  2:1034688/150‑1 (MQ=255)
                                                                                                           agaagGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGa                                         <  1:1158949/150‑1 (MQ=255)
                                                                                                             aagGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAg                                       <  2:765113/150‑1 (MQ=255)
                                                                                                               ggTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAATAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGAt                                     >  1:289600/1‑150 (MQ=255)
                                                                                                               ggTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGAt                                     <  2:253554/150‑1 (MQ=255)
                                                                                                                gTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGAtt                                    >  1:1270795/1‑150 (MQ=255)
                                                                                                                 tGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGAttt                                   <  1:1593018/150‑1 (MQ=255)
                                                                                                                 ggAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGAttt                                   <  2:1562964/149‑1 (MQ=255)
                                                                                                                    aCGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAAc                                <  1:891291/150‑1 (MQ=255)
                                                                                                                    aCGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAAc                                >  2:1387987/1‑150 (MQ=255)
                                                                                                                     cGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCATCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACt                               <  1:1242828/150‑1 (MQ=255)
                                                                                                                      gCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTg                              >  2:694345/1‑150 (MQ=255)
                                                                                                                       cTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTgg                             >  1:154269/1‑149 (MQ=255)
                                                                                                                        tATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCt                            <  2:769382/150‑1 (MQ=255)
                                                                                                                              gcgcTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTa                      <  1:1400715/150‑1 (MQ=255)
                                                                                                                               cgcTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTAc                     >  1:1324438/1‑150 (MQ=255)
                                                                                                                                 cagctgAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTg                   <  2:638417/148‑1 (MQ=255)
                                                                                                                                     tcgaaGTTAACTTCGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATc               >  2:1605657/4‑150 (MQ=255)
                                                                                                                                      gAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTAGCGTACTGGATCt              >  1:913854/1‑150 (MQ=255)
                                                                                                                                       gaaGTTAACTTCGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCt              <  1:812802/148‑1 (MQ=255)
                                                                                                                                       gaaGTTAACTTCGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATATAACTGCTCGCGTACTGGATCTg             >  1:1608753/2‑150 (MQ=255)
                                                                                                                                       aaaGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTg             >  1:324139/1‑150 (MQ=255)
                                                                                                                                        aaGTTAACGAAGTTAACATCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTgg            >  2:28134/1‑150 (MQ=255)
                                                                                                                                         aGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCACGCAACAAAATCATCTTTGAGAACTTAACCACGCTGCACGCAAACTCTAGTCTGCTTATGGAACGTGGTAAAGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATATGGc           >  1:1014993/1‑150 (MQ=255)
                                                                                                                                          gTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCa          >  1:421886/1‑150 (MQ=255)
                                                                                                                                           ttAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCAt         >  2:411225/1‑150 (MQ=255)
                                                                                                                                           tgaccGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAATATTTAACTGCTCGCGTACTGGATCTGGCAt         <  2:831548/146‑1 (MQ=255)
                                                                                                                                            tAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATc        >  2:922770/1‑150 (MQ=255)
                                                                                                                                             aaCGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCa       >  2:549805/1‑150 (MQ=255)
                                                                                                                                             aaCGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCa       >  2:1258464/1‑150 (MQ=255)
                                                                                                                                             aaCGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCa       >  1:611159/1‑150 (MQ=255)
                                                                                                                                              aCGAATTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCAc      <  2:348602/150‑1 (MQ=255)
                                                                                                                                              aCGAAGTTAACTTCGACAAACCTGAACACGCCCTCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTGGTCTGCGTATGGAACGTGGTAACGTTTCTACTGAAGATTTATCTGCTCGCGTACTGGATCTGGCATCAc      <  1:1030976/150‑1 (MQ=255)
                                                                                                                                               cGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCAcc     <  2:907948/150‑1 (MQ=255)
                                                                                                                                               cGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCAcc     >  2:422937/1‑150 (MQ=255)
                                                                                                                                               cGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCAcc     <  2:1435308/150‑1 (MQ=255)
                                                                                                                                                  aGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGGCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTAt  <  1:284990/150‑1 (MQ=255)
                                                                                                                                                  aGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACACCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTTAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTaa  >  2:1322525/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
AGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTAT  >  NC_000913/3966624‑3966919

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