Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,196,268 Δ2 bp intergenic (‑169/+63) thiC ← / ← rsd phosphomethylpyrimidine synthase/regulator of sigma D

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,196,2680C.94.9% 235.4 / 6.0 60intergenic (‑169/+64)thiC/rsdphosphomethylpyrimidine synthase/regulator of sigma D
Reads supporting (aligned to +/- strand):  ref base C (0/3);  new base . (36/20);  total (37/23)
Fisher's exact test for biased strand distribution p-value = 5.45e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009134,196,2690C.96.5% 234.9 / 3.0 60intergenic (‑170/+63)thiC/rsdphosphomethylpyrimidine synthase/regulator of sigma D
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base . (36/20);  minor base T (1/1);  total (37/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGC  >  NC_000913/4196139‑4196403
                                                                                                                                 ||                                                                                                                                      
acgGGTGAAGCAGGAGACGGTCGCCGCAGGGGCGGTGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATGAGCCTGATCAGGTTCCGCGGATCCCGAATAAACTGTCTCAGCCAGTTAAGGCACT‑‑GACAAGAAATTAGc                                                                                                                          <  2:340520/143‑1 (MQ=255)
  ggATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCCGGCTCCGCGGACCCCGCTTACCCGGTCTCTGCCAGTCAGGGACATGACAAGAAACTAGccc                                                                                                                        >  1:241233/1‑143 (MQ=255)
    aggAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACT‑‑GACAAGAAATTAGCCCCg                                                                                                                      <  1:302936/141‑1 (MQ=255)
        aGCAAGAGACGATCGCCGCAGGGGCGATGCATTAAATGTCGGATGCGAGTTCCCTTCGCAGGTATTAGCATGATCAGGGACCGCGGAGGCCGAAAAAACGTACTCAGCCAGTTAAGGCACT‑‑GACAAGAAATTAGCCCCGCaaa                                                                                                                  >  2:394837/1‑143 (MQ=255)
             gagaCGATCGCTGCAGGGGCGATGCGATGGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTTCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACT‑‑GACAAGAAATTAGCCCCGCAAATgggg                                                                                                             <  2:131424/143‑1 (MQ=255)
              agaCGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACT‑‑GACAAGAAATTAGCCCCGCAAATGGGGc                                                                                                            <  2:30787/143‑1 (MQ=255)
                  gATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCGTTAGGGCGCGGTACCCTGCTAAGGCGCCGCTGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACT‑‑GACAAGAAATTAGCCCCGCAAATGGGGCATTg                                                                                                        <  2:494325/143‑1 (MQ=255)
                  gATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACT‑‑GACAAGAAATTAGCCCCGCAAATGGGGCATTg                                                                                                        >  2:1047137/1‑143 (MQ=255)
                   aTCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACT‑‑GACAAGAAATTAGCCCCGCACATGGGGCATTGa                                                                                                       >  2:816596/1‑143 (MQ=255)
                          cAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTCCGCATGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACT‑‑GACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTaaa                                                                                                <  1:756950/143‑1 (MQ=255)
                             gggCGATGCGATAGCAGATTAGCTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACT‑‑GACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTa                                                                                             <  2:467171/143‑1 (MQ=255)
                                cgatgagatAGCAGATTAACGATTGTTCCCTTCGCAGGTATTATCCTGATCAGTTTCCGCGGATCCCGAATAACCGGTCTCAGCCAGTTAAGGCATT‑‑GACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTAcgc                                                                                          <  2:425695/143‑1 (MQ=255)
                                 gatgcgatAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACT‑‑GACAAGAAACTAGACCCGCAAGTGGGGCATGGAATGTAAATTCcgcg                                                                                         >  2:123509/1‑143 (MQ=255)
                                   tgcgatAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACT‑‑GACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGtt                                                                                       <  2:1020994/143‑1 (MQ=255)
                                         aGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACT‑‑GACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAgc                                                                                 >  2:420084/1‑143 (MQ=255)
                                                   ctctTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACT‑‑GACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCa                                                                       >  1:664757/1‑143 (MQ=255)
                                                    tctTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACT‑‑GACAAGAAATTAGCCCCGCAAATGAGGCATTGAATGTAAAGTACGCGTTAACAGCGCAAAACTCaa                                                                      >  2:835520/1‑143 (MQ=255)
                                                      ttGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACT‑‑GACAAGAAATTAGCCCCGCCAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGc                                                                    >  2:710556/1‑143 (MQ=255)
                                                         ttCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACT‑‑GACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAgg                                                                 >  2:715709/1‑143 (MQ=255)
                                                           cccTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACT‑‑GACAAGAAATTAGAGCGGCAAATGGGGCATGGAATGTAAATTACGCGTTAGCAGCGCAGAACACAAGCAGGAt                                                               >  1:933968/1‑143 (MQ=255)
                                                                 gCAGGTATTAGCCTGATCAGGTTAAGCGGATCCCGAATAAACGGTCTCAGCAAGTTAAGGCACT‑‑GACAAGAAAGTAGCGCGGCAAATGGGAAGTAGAATGTAAATTACGCGTAAACAGCGCAGAACTCAAGCAGAATGGTTGa                                                         >  2:1072423/1‑143 (MQ=255)
                                                                   aGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACT‑‑GACAAGAAATTAGCCCCGCAAATGGGGCATCGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGAcg                                                       >  2:659317/1‑143 (MQ=255)
                                                                     gTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACT‑‑GACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGCcgcg                                                     >  1:393242/1‑143 (MQ=255)
                                                                      gatgaGCCTGATCAGGTGCCGCGGATCCCGAATAAACGGGCTCAGCCAGTTAAGGAA‑‑CAGTTAAGAACTTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACgcgg                                                    <  2:668876/139‑1 (MQ=255)
                                                                       aTTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACT‑‑GACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACgcggg                                                   >  1:833094/1‑143 (MQ=255)
                                                                           gCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACT‑‑GACAAGAAATTAGCCCCTCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGggcggc                                               <  1:102078/143‑1 (MQ=255)
                                                                           gCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACT‑‑GACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGggcggc                                               >  2:514009/1‑143 (MQ=255)
                                                                           gCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACT‑‑GACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGggcggc                                               <  1:1003875/143‑1 (MQ=255)
                                                                            ccTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACT‑‑GACAAGAAATTAGCCCCGCAAATGGGGCAGTGAATGTAAAGTACGCGCTAACAGCGCAGAACGCAAGCAGGATGGTTGACgagggcgggg                                              >  1:332454/1‑141 (MQ=255)
                                                                             cTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACT‑‑GACCAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGCTGCTTGACGCGGGCGGCGt                                             >  1:453463/1‑143 (MQ=255)
                                                                             cTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACT‑‑GACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGt                                             >  1:577469/1‑143 (MQ=255)
                                                                                 tCAGGGTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACT‑‑GACAAGAAATTCGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAg                                         <  2:549963/143‑1 (MQ=255)
                                                                                  cAGGTTCCGCGGATCCGGAGTAGGCGGTCTCGGCGGGTTAAGG‑‑GTCTGGCGGGCGGTTAGGCCCGCGGTTGGGGCATTGAATGTAAATTACGCGTTGGCGGCGCCGCTCTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGc                                        <  2:21353/143‑1 (MQ=255)
                                                                                    ggTTCCGCGGTTGCCGAATAAACGGTCTCAGCCAGTTAAGGCACT‑‑GACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCAc                                      <  2:576180/143‑1 (MQ=255)
                                                                                     gTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACT‑‑GACAAGAAATTAGCCCCGCAAATGGGGCATTGGATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCAcc                                     >  1:141917/1‑143 (MQ=255)
                                                                                       tCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACT‑‑GACTAGAAATTAGCCCCGCAACTGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACcag                                   >  2:287432/1‑143 (MQ=255)
                                                                                         cgcgGATCCCGAATAAACGGTCTGAGCCAGTTAAGGCACT‑‑GACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACcagca                                 <  1:394325/143‑1 (MQ=255)
                                                                                         cgcgGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACT‑‑GACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACcagca                                 >  1:502592/1‑143 (MQ=255)
                                                                                              aTCCCGAATAAACGGTCTCAGCCAGTTAAGGCACT‑‑GACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATc                            <  1:622784/143‑1 (MQ=255)
                                                                                                 ccGAATAAACGGTCTCAGCCAGTTAAGGCACT‑‑GACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAGCTCAAGCAGGATGTTTGACGCGCGCGGCGTCACGCACCAGCAGCATCGGc                         >  2:105646/1‑143 (MQ=255)
                                                                                                  cGAATAAACGGTCTCAGCCAGTTAAGGCACT‑‑GACCAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGAGCAGCGTCAAGCACCAGCAGAATCAGCt                        >  2:412761/1‑143 (MQ=255)
                                                                                                  cGAATAAACGGTCTCAGCCAGTTAAGGCACT‑‑GACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCt                        >  1:1063152/1‑143 (MQ=255)
                                                                                                  cGAATAAACGGTCTCAGCCAGTTAAGGCACT‑‑GACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCt                        >  2:79645/1‑143 (MQ=255)
                                                                                                   gAATAAACGGTCTCAGCCAGTTAAGGCACT‑‑GACAAGAAAGTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCtt                       <  1:394837/143‑1 (MQ=255)
                                                                                                   gAATAAACGGTCTCAGCCAGATAAGGCACT‑‑GACAAGAAATTAGCCCCGCAAATGGGGAATTGAATGTAAATTACGCGTTAACAGCGCAGCGCTCAAGCAGCATGTTTGACGCGGGCGGCGTAAAGCACCAGCAGAATCAGCta                       >  1:535347/1‑142 (MQ=255)
                                                                                                        aaCGGTCTCAGCCAGTTAAGGCACT‑‑GACGAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACACCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCCCCTGCAGAATCAGcgcacctt                  >  2:475449/1‑136 (MQ=255)
                                                                                                         aCGGTCTCAGCCAGTTAAGGCACT‑‑GACAAGAAATTAGCCCCGCAAATGGTGCATTGAATGTAAATTACGGGATAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTc                 >  1:911166/1‑143 (MQ=255)
                                                                                                         aCGGTCTCAGCCAGTTAAGGCACT‑‑GACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTc                 >  1:865654/1‑143 (MQ=255)
                                                                                                         aCGGTCTCAGCCAGTTAAGGCACT‑‑GACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTc                 >  1:176174/1‑143 (MQ=255)
                                                                                                          cGGTCTCAGCCAGTTAAGGCACT‑‑GACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTcc                <  2:815613/143‑1 (MQ=255)
                                                                                                              ctcAGCCAGTTAAGGCCCGGACAAGAGAGTTGCCCGGCAAGTGGGGCATTGAATGTAAATTGCGCGTTGACAGCGGGGAACTCAAGCGGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGca            <  2:186845/143‑1 (MQ=255)
                                                                                                                cAGCCAGTTAAGGCACT‑‑GACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCaa                                                                      <  2:958404/83‑1 (MQ=255)
                                                                                                                cAGCCAGTTAAGGCACT‑‑GACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCaa                                                                      >  1:958404/1‑83 (MQ=255)
                                                                                                                 aGCCAGTTAAGGCACT‑‑GACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGcacca         >  1:825426/1‑141 (MQ=255)
                                                                                                                 aGCCAGTTAAGGCACT‑‑GACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACaa         >  1:404433/1‑143 (MQ=255)
                                                                                                                 aGCCAGTTAAGGCACT‑‑GACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACaa         >  2:810993/1‑143 (MQ=255)
                                                                                                                 aGCCAGTTAAGGCACT‑‑GACAAGAAATTAGCACCGCAAATGGGGCATTGCATGTAAATTACGCGATAACAGCGCAGAACTCAAGCACGATGTTTGACGCGGGCCGCGGACAGCCCCCACAGAACCACCGTAACTTCCAGCACaa         >  2:214500/1‑143 (MQ=255)
                                                                                                                      gTTAAGGCACT‑‑GACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAgcgc    >  1:591445/1‑143 (MQ=255)
                                                                                                                      gTTAAGGCACT‑‑GACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAgcgc    <  1:593315/143‑1 (MQ=255)
                                                                                                                        tAAGGCACT‑‑GACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAgcgcgc  <  2:833094/143‑1 (MQ=255)
                                                                                                                                 ||                                                                                                                                      
ACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGC  >  NC_000913/4196139‑4196403

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