Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 969,798 T→C L396P (CTC→CCC)  ycaQ → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913969,7980TC99.5% 626.7 / ‑5.8 196L396P (CTC→CCC) ycaQhypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (85/110);  minor base G (1/0);  total (86/110)
Fisher's exact test for biased strand distribution p-value = 4.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 201 total aligned reads displayed.

GCTGCAAAAAGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCTCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCGATCCATCTG  >  NC_000913/969709‑969887
                                                                                         |                                                                                         
gCTGCAAAAAGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCtt                                                                                        >  1:4750141/1‑93 (MQ=255)
gCTGCAAAAAGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCtt                                                                                        <  1:6296313/93‑1 (MQ=255)
gCTGCAAAAAGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCtt                                                                                        <  2:2107883/93‑1 (MQ=255)
 cTGCAAAAAGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCttt                                                                                       <  2:2441209/93‑1 (MQ=255)
  tGCAAAAAGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTa                                                                                      >  2:5016524/1‑93 (MQ=255)
  tGCAAAAAGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTa                                                                                      >  2:414236/1‑93 (MQ=255)
      aaaaGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGa                                                                                  >  1:2305488/1‑93 (MQ=255)
       aaaGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGAt                                                                                 <  2:53553/93‑1 (MQ=255)
       aaaGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGAt                                                                                 <  2:4684592/93‑1 (MQ=255)
        aaGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGAtt                                                                                <  2:1185979/93‑1 (MQ=255)
         aGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTg                                                                               <  2:287892/93‑1 (MQ=255)
          gggCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGc                                                                              >  1:1486452/1‑93 (MQ=255)
          gggCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGc                                                                              >  1:5612138/1‑93 (MQ=255)
          gggCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCACGCAAGGCCCCTTTACTGATTGc                                                                              >  1:3108304/1‑93 (MQ=255)
           ggCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGcc                                                                             <  1:67910/93‑1 (MQ=255)
           ggCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGcc                                                                             <  1:2975519/93‑1 (MQ=255)
              tACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCa                                                                          >  2:1852098/1‑93 (MQ=255)
               aCGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCAc                                                                         <  1:3763782/93‑1 (MQ=255)
               aCGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCAc                                                                         <  2:823556/93‑1 (MQ=255)
               aCGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCAc                                                                         >  1:1039437/1‑93 (MQ=255)
                cGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCAcc                                                                        >  1:1258678/1‑93 (MQ=255)
                  tCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCgg                                                                      <  1:1119181/93‑1 (MQ=255)
                    aGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCt                                                                    >  1:3300818/1‑93 (MQ=255)
                      gCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTgg                                                                  <  2:1517077/93‑1 (MQ=255)
                      gCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTgg                                                                  <  1:3336917/93‑1 (MQ=255)
                        gATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGa                                                                >  2:1983254/1‑93 (MQ=255)
                        gATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGa                                                                <  1:4047334/93‑1 (MQ=255)
                        gATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGa                                                                >  1:2421387/1‑93 (MQ=255)
                         aTTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGaa                                                               >  2:1896605/1‑93 (MQ=255)
                         aTTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGaa                                                               >  1:1594807/1‑93 (MQ=255)
                          ttACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGaaa                                                              >  2:4033191/1‑93 (MQ=255)
                           tACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAAt                                                             >  1:3574844/1‑93 (MQ=255)
                           tACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAAt                                                             <  2:706719/93‑1 (MQ=255)
                            aCTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATa                                                            <  2:2613006/93‑1 (MQ=255)
                             cTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAg                                                           >  2:4252215/1‑93 (MQ=255)
                             cTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAg                                                           >  2:4418027/1‑93 (MQ=255)
                             cTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAg                                                           <  2:4602579/93‑1 (MQ=255)
                             cTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAg                                                           >  1:418995/1‑93 (MQ=255)
                              tGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGa                                                          >  1:4213225/1‑93 (MQ=255)
                               gATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGAc                                                         >  1:4627179/1‑93 (MQ=255)
                               gATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGAc                                                         >  1:947742/1‑93 (MQ=255)
                                aTTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGAcc                                                        >  1:5494708/1‑93 (MQ=255)
                                aTTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGAcc                                                        <  2:2044889/93‑1 (MQ=255)
                                 tttCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGAccc                                                       >  2:299167/1‑93 (MQ=255)
                                 tttCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGAccc                                                       <  2:1076722/93‑1 (MQ=255)
                                 tgtCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGAccc                                                       <  2:4900416/91‑1 (MQ=255)
                                  ttCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGAcccc                                                      <  2:3900070/93‑1 (MQ=255)
                                    cGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGt                                                    <  2:5908435/93‑1 (MQ=255)
                                    cGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGt                                                    <  1:4341843/93‑1 (MQ=255)
                                     gCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTc                                                   >  2:382805/1‑93 (MQ=255)
                                      cTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCg                                                  >  2:4483926/1‑93 (MQ=255)
                                       tAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGc                                                 <  2:4830710/93‑1 (MQ=255)
                                       tAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGc                                                 >  2:3380142/1‑93 (MQ=255)
                                        aaCTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCa                                                <  1:3433152/93‑1 (MQ=255)
                                        aaCTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCa                                                >  1:5436661/1‑93 (MQ=255)
                                        aaCTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCa                                                <  1:281777/93‑1 (MQ=255)
                                         aCTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCAt                                               <  2:5583263/93‑1 (MQ=255)
                                            ggCAGCAGGCAACGCGGGTGACATTAGGATTCTGCCCGCAAGCCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATaaa                                            <  1:4099269/93‑1 (MQ=255)
                                            ggCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATaaa                                            <  1:5796646/93‑1 (MQ=255)
                                             gcagcaGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTGTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAAt                                           <  1:73441/93‑1 (MQ=255)
                                              cagcagGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATg                                          <  1:824449/93‑1 (MQ=255)
                                              cagcagGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATg                                          >  2:27603/1‑93 (MQ=255)
                                                gcagGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGAt                                        <  1:1681268/93‑1 (MQ=255)
                                                gcagGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGAt                                        >  1:4886551/1‑93 (MQ=255)
                                                gcagGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGAt                                        <  1:1887909/93‑1 (MQ=255)
                                                 cagGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGAtt                                       >  2:2368555/1‑93 (MQ=255)
                                                  agGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTa                                      <  2:2576043/93‑1 (MQ=255)
                                                  agGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTa                                      <  1:2026683/93‑1 (MQ=255)
                                                   ggCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTAt                                     <  2:5208597/93‑1 (MQ=255)
                                                   ggCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTAt                                     >  2:1161429/1‑93 (MQ=255)
                                                   ggCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTAt                                     >  1:3484149/1‑93 (MQ=255)
                                                       aCGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATa                                 >  1:1277213/1‑93 (MQ=255)
                                                         gcgGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAg                               <  2:2146481/93‑1 (MQ=255)
                                                           gggTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGAt                             >  2:1801873/1‑93 (MQ=255)
                                                            ggTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGAtt                            >  2:5408725/1‑93 (MQ=255)
                                                            ggTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGAtt                            >  2:3249000/1‑93 (MQ=255)
                                                             gTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGAttt                           >  2:5501201/1‑93 (MQ=255)
                                                              tGACATTAGTATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTa                          <  1:1453055/93‑1 (MQ=255)
                                                              tGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTa                          <  2:5156821/93‑1 (MQ=255)
                                                               gACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTAt                         >  1:4976914/1‑93 (MQ=255)
                                                               gACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCAGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTAt                         <  1:5994744/93‑1 (MQ=255)
                                                                aCATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATg                        <  2:1474551/93‑1 (MQ=255)
                                                                 cATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGc                       <  1:2265137/93‑1 (MQ=255)
                                                                  aTTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGCCCCCGTCGCATAAATGATTATGATAAGATTTATGCa                      >  2:6233991/1‑93 (MQ=255)
                                                                   ttAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCAt                     <  2:3473464/93‑1 (MQ=255)
                                                                      ggATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCa                  <  1:5079912/93‑1 (MQ=255)
                                                                      ggATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCa                  >  1:2619213/1‑93 (MQ=255)
                                                                       gATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCAt                 >  2:1289607/1‑93 (MQ=255)
                                                                            tgcCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAAccc            <  2:610552/93‑1 (MQ=255)
                                                                             gcCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTTATAACACg           >  2:3426927/1‑93 (MQ=255)
                                                                             gcCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCg           <  1:200836/93‑1 (MQ=255)
                                                                              cccGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCGa          <  2:2581939/93‑1 (MQ=255)
                                                                              cccGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCGa          <  2:5778311/93‑1 (MQ=255)
                                                                               ccGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCGAt         >  2:1314366/1‑93 (MQ=255)
                                                                               ccGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCGAt         <  2:3524899/93‑1 (MQ=255)
                                                                                  cAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCGATCCa      <  2:6264992/93‑1 (MQ=255)
                                                                                   aaGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCGATCCAt     <  2:4186560/93‑1 (MQ=255)
                                                                                    aGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCGATCCATc    >  1:5984064/1‑93 (MQ=255)
                                                                                     ggCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCGATCCATCt   <  2:1488408/93‑1 (MQ=255)
                                                                                     ggCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCGATCCATCt   >  1:1783532/1‑93 (MQ=255)
                                                                                         |                                                                                         
GCTGCAAAAAGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCTCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCGATCCATCTG  >  NC_000913/969709‑969887

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