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.6% 971.5 / ‑5.7 270T306T (ACG→ACAb3613putative membrane protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (104/165);  total (105/165)
Fisher's exact test for biased strand distribution p-value = 3.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.42e-01

Only 100 of 284 total aligned reads displayed.

ACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACGCTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAGGTACTGAAGTTAAAGCGAT  >  NC_000913/3785292‑3785471
                                                                                          |                                                                                         
accgaccgAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACt                                                                                         >  2:4321063/1‑93 (MQ=255)
accgaccgAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACt                                                                                         >  2:2823785/1‑93 (MQ=255)
 ccgaccgaAAGCGAAAAACCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTg                                                                                        <  2:108610/93‑1 (MQ=255)
   gaccgaAAGCGAAAAATCGCCGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCa                                                                                      >  2:3188079/1‑93 (MQ=255)
     ccgaAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATc                                                                                    >  2:824771/1‑93 (MQ=255)
     ccgaAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATc                                                                                    <  1:2062445/93‑1 (MQ=255)
     ccgaAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGGGGTCAGGCATTCTGGCCGGTTCGCGGGCCGAGAGTGCATc                                                                                    >  2:3565711/1‑93 (MQ=255)
      cgaAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCg                                                                                   >  1:1533836/1‑93 (MQ=255)
      cgaAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCg                                                                                   >  2:2259/1‑93 (MQ=255)
      cgaAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCg                                                                                   >  2:3492934/1‑93 (MQ=255)
       gaAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGc                                                                                  <  1:2138115/93‑1 (MQ=255)
        aaaGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCt                                                                                 <  1:3495699/93‑1 (MQ=255)
        aaaGCGAAAAATCGCTGATGTCCCGTACCGGGGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGGTCGCGGGCCGACACTGCATCGCt                                                                                 >  2:3269492/1‑93 (MQ=255)
         aaGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTa                                                                                <  1:3370282/93‑1 (MQ=255)
          aGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTAt                                                                               >  2:2993132/1‑93 (MQ=255)
          aGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTAt                                                                               >  1:188251/1‑93 (MQ=255)
          aGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTAt                                                                               <  1:570861/93‑1 (MQ=255)
           gCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATg                                                                              >  1:4235325/1‑93 (MQ=255)
             gAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGc                                                                            <  1:1784200/93‑1 (MQ=255)
              aaaaaTCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCg                                                                           >  1:2357361/1‑93 (MQ=255)
                aaaTCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGaa                                                                         >  1:4557996/1‑93 (MQ=255)
                aaaTCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGaa                                                                         <  1:4460864/93‑1 (MQ=255)
                aaaTCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGaa                                                                         <  1:1018831/93‑1 (MQ=255)
                 aaTCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCGGGTTCGCGGGCCGACACTGCATCGCTATGGCGAAc                                                                        <  1:4115129/93‑1 (MQ=255)
                 aaTCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAAc                                                                        >  1:3395297/1‑93 (MQ=255)
                  aTCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACa                                                                       <  2:2675113/93‑1 (MQ=255)
                  aTCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACa                                                                       >  2:1446490/1‑93 (MQ=255)
                  aTCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACa                                                                       <  2:775912/93‑1 (MQ=255)
                  aTCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACa                                                                       <  1:1171683/93‑1 (MQ=255)
                   tCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAg                                                                      >  2:890598/1‑93 (MQ=255)
                   tCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAg                                                                      >  2:1612959/1‑93 (MQ=255)
                    cGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGc                                                                     <  2:2456591/93‑1 (MQ=255)
                    cGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGc                                                                     <  2:3860858/93‑1 (MQ=255)
                       tGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTAc                                                                  <  2:1697037/93‑1 (MQ=255)
                       tGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTAc                                                                  <  2:3196457/93‑1 (MQ=255)
                       tGATGCCCCGTACCGGTGGTCGGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTAc                                                                  <  2:321011/93‑1 (MQ=255)
                        gATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCGGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACa                                                                 <  1:901349/93‑1 (MQ=255)
                        gATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACa                                                                 >  2:1590368/1‑93 (MQ=255)
                         aTGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAg                                                                >  1:171519/1‑93 (MQ=255)
                         aTGGCCCGTACCGGTGGTCGGGGCGCCCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAg                                                                <  1:2025996/93‑1 (MQ=255)
                           gTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAggg                                                              >  2:1990439/1‑93 (MQ=255)
                           gTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAggg                                                              >  1:1810070/1‑93 (MQ=255)
                           gTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGCCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAggg                                                              <  2:559420/93‑1 (MQ=255)
                            tCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGt                                                             <  1:766906/93‑1 (MQ=255)
                            tCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGt                                                             <  2:2896441/93‑1 (MQ=255)
                              ccGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGa                                                           >  2:1354991/1‑93 (MQ=255)
                               cGTACGGGTGGCCGGGGCGCCCGGCGTGGCAAGGCATTCTGGCCGCTCCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGaa                                                          <  1:2820063/93‑1 (MQ=255)
                                 tACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAAtt                                                        <  1:4376591/93‑1 (MQ=255)
                                 tACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAAtt                                                        <  2:1310200/93‑1 (MQ=255)
                                 tACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAAtt                                                        <  1:961595/93‑1 (MQ=255)
                                  gccGGTGGTGTGGGCGCGCTGCTTGGTCAGGCATTCGGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTa                                                       <  2:4174895/92‑1 (MQ=255)
                                   ccGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTAc                                                      <  2:2749486/93‑1 (MQ=255)
                                   ccGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTAc                                                      <  1:752140/93‑1 (MQ=255)
                                   ccGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTAc                                                      <  1:2883574/93‑1 (MQ=255)
                                     ggtggtCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGc                                                    <  1:625866/93‑1 (MQ=255)
                                       gggtCGGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTg                                                  <  1:3167389/92‑1 (MQ=255)
                                        ggtCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGGCGACACTGTATAGCTATGGCGAACAGCTGCAGGGTGAATTACGCTgg                                                 >  2:4663018/1‑93 (MQ=255)
                                          tCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGaa                                               >  1:2382457/1‑93 (MQ=255)
                                           cTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGaaa                                              <  1:2626751/93‑1 (MQ=255)
                                             gggCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAgg                                            <  2:1504447/93‑1 (MQ=255)
                                              ggCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGt                                           <  2:1007317/93‑1 (MQ=255)
                                                cgctcCGCGTGTGAAGACCTTCTGGCCGCTTAGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTAt                                         <  2:1143357/93‑1 (MQ=255)
                                                 gggccgcGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATg                                        <  1:3988276/91‑1 (MQ=255)
                                                 gcgcCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGCATTACGCTGGAAAGGTATg                                        <  2:4524121/93‑1 (MQ=255)
                                                 gcgcCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATg                                        <  1:1764151/93‑1 (MQ=255)
                                                 gcgcCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATg                                        >  1:1861824/1‑93 (MQ=255)
                                                    ccgcGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGtt                                     >  1:3884535/1‑93 (MQ=255)
                                                      gcgTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTAt                                   <  2:2626516/93‑1 (MQ=255)
                                                        gTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCg                                 <  1:1686835/93‑1 (MQ=255)
                                                         tGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCgg                                <  2:643495/93‑1 (MQ=255)
                                                         tGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCgg                                <  2:502469/93‑1 (MQ=255)
                                                          gttCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGt                               <  2:3354709/91‑1 (MQ=255)
                                                          ggTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGACAGGTATGGTTATCgct                               >  2:232430/1‑91 (MQ=255)
                                                           gTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTg                              >  1:2965286/1‑93 (MQ=255)
                                                            tCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGc                             >  2:3705850/1‑93 (MQ=255)
                                                            tCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGc                             >  1:270583/1‑93 (MQ=255)
                                                             cAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCt                            >  2:3671911/1‑93 (MQ=255)
                                                             cAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCt                            <  1:663381/93‑1 (MQ=255)
                                                              aGGCATTCGGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCtt                           <  1:1282107/93‑1 (MQ=255)
                                                                  aTTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGa                       >  2:4005976/1‑93 (MQ=255)
                                                                   ttCTGGCCGGTTGCCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAATGTATGGTTATCGGTGCTTCTGaa                      <  2:579487/93‑1 (MQ=255)
                                                                   ttCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGGGAATTACGCTGGAAAGGTATGGTTATAGGTGCTTCTGaa                      >  2:1625951/1‑93 (MQ=255)
                                                                    tCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAg                     <  1:3819143/93‑1 (MQ=255)
                                                                    tCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAg                     >  1:3295889/1‑93 (MQ=255)
                                                                    tCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAg                     <  2:1881981/93‑1 (MQ=255)
                                                                     cTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAgg                    <  1:323919/93‑1 (MQ=255)
                                                                        gCCGGTCCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAGGTAc                 <  2:3801502/93‑1 (MQ=255)
                                                                         ccGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAGGTACt                <  1:2648593/93‑1 (MQ=255)
                                                                            gTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAGGTACTGaa             >  2:3414542/1‑93 (MQ=255)
                                                                             ttCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAGGTACTGAAg            <  1:3064554/93‑1 (MQ=255)
                                                                             ttCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAGGTACTGAAg            <  2:4184744/93‑1 (MQ=255)
                                                                             ttCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAGGTACTGAAg            <  1:424088/93‑1 (MQ=255)
                                                                             ttCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAGGTACTGAAg            <  2:2160644/93‑1 (MQ=255)
                                                                              tCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAGGTACTGAAGt           <  1:228356/93‑1 (MQ=255)
                                                                                gcgGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAGGTACTGAAGTTa         >  1:1819305/1‑93 (MQ=255)
                                                                                  gggCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAGGTACTGAAGTTaaa       <  1:3640959/93‑1 (MQ=255)
                                                                                   ggCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAGGTACTGAAGTTAAAg      <  2:2139159/93‑1 (MQ=255)
                                                                                    gCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAGGTACTGAAGTTAAAGc     <  2:3105340/93‑1 (MQ=255)
                                                                                    gCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAGGTACTGAAGTTAAAGc     <  2:732/93‑1 (MQ=255)
                                                                                      cGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAGGTACTGAAGTTAAAGCGa   <  1:3584563/93‑1 (MQ=255)
                                                                                          |                                                                                         
ACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACGCTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAGGTACTGAAGTTAAAGCGAT  >  NC_000913/3785292‑3785471

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