Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,785,382 G→A T306T (ACG→ACA b3613 → putative membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,785,3820GA99.5% 1317.4 / ‑4.6 365T306T (ACG→ACAb3613putative membrane protein
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (171/191);  total (172/193)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

Only 100 of 379 total aligned reads displayed.

ACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACGCTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAGGTACTGAAGTTAAAGCGAT  >  NC_000913/3785292‑3785471
                                                                                          |                                                                                         
accgaccgAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACt                                                                                         <  1:5799023/93‑1 (MQ=255)
  cgaccgaAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGc                                                                                       >  1:1153398/1‑93 (MQ=255)
  cgaccgaAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGc                                                                                       >  2:6468857/1‑93 (MQ=255)
  cgaccgaAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGc                                                                                       <  1:3427345/93‑1 (MQ=255)
    accgaAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCAt                                                                                     >  1:6980409/1‑93 (MQ=255)
     ccgaAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATc                                                                                    <  1:2886536/93‑1 (MQ=255)
        aaaGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCt                                                                                 >  2:6658110/1‑93 (MQ=255)
         aaGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCtt                                                                                >  2:7245029/1‑92 (MQ=255)
         aaGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTa                                                                                >  2:9018526/1‑93 (MQ=255)
         aaGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTa                                                                                <  1:3625692/93‑1 (MQ=255)
         aaGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTa                                                                                >  2:5579562/1‑93 (MQ=255)
         aaGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTa                                                                                >  1:9376790/1‑93 (MQ=255)
          aGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTAt                                                                               <  2:3958617/93‑1 (MQ=255)
           gCGAAAAATCGCTGATGTCCCGTACCGGTGGGCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGCCCCTGTCTCGCTCTg                                                                              >  1:6372018/1‑93 (MQ=255)
            cGAAAAATCGCTGATGTCCCGTCCCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATgg                                                                             <  2:6354701/93‑1 (MQ=255)
            cGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATgg                                                                             >  1:6748570/1‑93 (MQ=255)
            cGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATgg                                                                             <  1:10029962/93‑1 (MQ=255)
               aaaaTCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGa                                                                          >  1:3987017/1‑93 (MQ=255)
                aaaTCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGaa                                                                         >  2:9549576/1‑93 (MQ=255)
                 actCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGAGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAAc                                                                        <  2:8201467/91‑1 (MQ=255)
                 aaTCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAAc                                                                        <  1:4690175/93‑1 (MQ=255)
                 aaTCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAAc                                                                        >  2:927205/1‑93 (MQ=255)
                  aTCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACa                                                                       <  2:3718827/93‑1 (MQ=255)
                  aTCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGCCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACa                                                                       <  1:7255134/93‑1 (MQ=255)
                   tCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACac                                                                      >  1:2573470/1‑92 (MQ=255)
                      cTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTa                                                                   >  1:5030118/1‑93 (MQ=255)
                      cTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGAGGGCCGACACTGCATCGCTATGGCGAACAGCTa                                                                   >  1:9731857/1‑93 (MQ=255)
                       tGATGTCCCGTGCCGGTGGTGTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTAc                                                                  <  1:182801/93‑1 (MQ=255)
                       tGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTAc                                                                  <  2:5487983/93‑1 (MQ=255)
                       tGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTAc                                                                  >  2:4827330/1‑93 (MQ=255)
                       tGATGTCCCGTACCGGTGGGCTGGGCGCGCCGCGTGATCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAACTAc                                                                  >  2:7326016/1‑93 (MQ=255)
                         aTGTCCCGTCCCGGGGTCCTGGGCGCGCCGCGTGTTCGGGCATTCTGGCGGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAg                                                                <  2:2499968/93‑1 (MQ=255)
                         aTGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGGCGCAACTGCCTCGATGTGGGGAACAGACACAg                                                                >  1:7147564/1‑93 (MQ=255)
                          tGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAgg                                                               >  1:8613404/1‑93 (MQ=255)
                          tGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAgg                                                               <  1:1936855/93‑1 (MQ=255)
                          tGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAgg                                                               >  2:2754325/1‑93 (MQ=255)
                          tGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAgg                                                               <  1:4084160/93‑1 (MQ=255)
                              ccGTACCGGTGGTCTGGGCGGGCGGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGa                                                           <  1:2337159/93‑1 (MQ=255)
                               cGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGaa                                                          >  2:9479183/1‑93 (MQ=255)
                                gTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAAt                                                         >  1:611158/1‑93 (MQ=255)
                                  gcgggtggtCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTa                                                       <  1:8147587/90‑1 (MQ=255)
                                   ccGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTAc                                                      <  1:2164743/93‑1 (MQ=255)
                                   ccGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGCGAATTAc                                                      >  1:8471261/1‑93 (MQ=255)
                                   ccGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCGGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTAc                                                      <  2:411009/93‑1 (MQ=255)
                                    cGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACg                                                     >  2:5887045/1‑93 (MQ=255)
                                     ggtggtCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGc                                                    >  1:7535070/1‑93 (MQ=255)
                                     ggtggtCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGc                                                    >  2:909461/1‑93 (MQ=255)
                                     ggtggtCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGc                                                    >  2:2782/1‑93 (MQ=255)
                                     ggtggtCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGCGACTTACGc                                                    >  1:4106389/1‑93 (MQ=255)
                                     gggggtCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGc                                                    <  1:8869509/90‑1 (MQ=255)
                                       tggtCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTg                                                  >  1:9966117/1‑93 (MQ=255)
                                       gggtCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTg                                                  <  2:7423122/92‑1 (MQ=255)
                                          tCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGaa                                               <  2:6118435/93‑1 (MQ=255)
                                           cTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGaaa                                              <  2:5104088/93‑1 (MQ=255)
                                           cTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGaaa                                              <  1:1482399/93‑1 (MQ=255)
                                            tGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAg                                             <  2:9407262/93‑1 (MQ=255)
                                             gggCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAgg                                            <  2:4205517/93‑1 (MQ=255)
                                               gcgcgcCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTa                                          >  2:1051047/1‑93 (MQ=255)
                                               gcgcgcCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTa                                          >  1:4902350/1‑93 (MQ=255)
                                                cgcgcCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTAt                                         <  2:3627603/93‑1 (MQ=255)
                                                  cgccgcGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATgg                                       >  2:2950317/1‑93 (MQ=255)
                                                   gccgcGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGt                                      >  1:5268494/1‑93 (MQ=255)
                                                     cgcgTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGGGAATTACGCTGGAAAGGTATGGTTa                                    >  2:2147368/1‑93 (MQ=255)
                                                      gcgTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTAt                                   <  1:4377527/93‑1 (MQ=255)
                                                      gcgTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTAt                                   <  1:2753939/93‑1 (MQ=255)
                                                      gcgTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTAt                                   >  1:830759/1‑93 (MQ=255)
                                                        gTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCg                                 <  1:1728422/93‑1 (MQ=255)
                                                         tagTCAGGCATTGTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCgg                                <  2:8689000/91‑1 (MQ=255)
                                                         tGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCgg                                >  2:3558771/1‑93 (MQ=255)
                                                         tGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCgg                                >  2:619310/1‑93 (MQ=255)
                                                          ggTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGt                               >  1:3892180/1‑93 (MQ=255)
                                                          ggTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGt                               >  1:1206557/1‑93 (MQ=255)
                                                           gTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTg                              <  2:3288151/93‑1 (MQ=255)
                                                              aGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCtt                           >  2:7078032/1‑93 (MQ=255)
                                                               ggCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTc                          <  2:2224208/93‑1 (MQ=255)
                                                               ggCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTc                          >  2:5197109/1‑93 (MQ=255)
                                                                   ttCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGaa                      <  1:785413/93‑1 (MQ=255)
                                                                   ttCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGaa                      <  2:7781350/93‑1 (MQ=255)
                                                                    tCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAg                     <  1:871690/93‑1 (MQ=255)
                                                                    tCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAg                     >  1:631358/1‑93 (MQ=255)
                                                                    tCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAg                     <  1:8943114/93‑1 (MQ=255)
                                                                      tGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAGGt                   >  2:1871448/1‑93 (MQ=255)
                                                                       ggCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAGGTa                  <  2:8374654/93‑1 (MQ=255)
                                                                           ggTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAGGTACTGa              >  1:7755683/1‑93 (MQ=255)
                                                                           ggTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGAATGGTCATCGGTGCTTCTGAAGGTACTGa              >  2:6932168/1‑93 (MQ=255)
                                                                             ttCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAGGTACTGAAg            <  2:5331901/93‑1 (MQ=255)
                                                                               cgcgGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAGGTACTGAAGtt          <  2:7570824/93‑1 (MQ=255)
                                                                                gcgGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAGGTACTGAAGTTa         <  2:7980298/93‑1 (MQ=255)
                                                                                gcgGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAGGTACTGAAGTTa         <  1:5351069/93‑1 (MQ=255)
                                                                                gcgGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAGGTACTGAAGTTa         <  1:3836153/93‑1 (MQ=255)
                                                                                 cgGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAGGTACTGAAGTTaa        <  1:1864980/93‑1 (MQ=255)
                                                                                 cgGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAGGTACTGAAGTTaa        >  2:4754945/1‑93 (MQ=255)
                                                                                 cgGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAGGTACTGAAGTTaa        >  1:5550175/1‑93 (MQ=255)
                                                                                  gggCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAGGTACTGAAGTTaaa       <  2:4348190/93‑1 (MQ=255)
                                                                                   ggCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAGGTACTGAAGTTAAAg      <  2:4535965/93‑1 (MQ=255)
                                                                                   ggCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAGGTACTGAAGTTAAAg      >  2:879491/1‑93 (MQ=255)
                                                                                    gCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAGGTACTGAAGTTAAAGc     <  1:7725366/93‑1 (MQ=255)
                                                                                     ccGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAGGTACTGAAGTTAAAGCg    <  1:2801211/93‑1 (MQ=255)
                                                                                      cGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAGGTACTGAAGTTAAAGCGa   >  1:3103495/1‑93 (MQ=255)
                                                                                       gACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAGGTACTGAAGTTAAAGCGAt  >  2:2631682/1‑93 (MQ=255)
                                                                                          |                                                                                         
ACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACGCTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAGGTACTGAAGTTAAAGCGAT  >  NC_000913/3785292‑3785471

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