Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,809,082 T→C V42A (GTT→GCT)  b1729 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,809,0820TC88.0% 528.7 / 51.6 191V42A (GTT→GCT) b1729hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (9/14);  new base C (80/88);  total (89/102)
Fisher's exact test for biased strand distribution p-value = 5.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

Only 100 of 199 total aligned reads displayed.

TCGTTGTACTGCTGTTTGCGCTGGCTCAGACCCGCCATAAACAGTGGAGTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGTTTTTGGCCTTGCCCTGCATACCATTTATGGTTCTGACAGCCAGGTACTTAAAGATTCTGTACAGTGGTTTAACATCGTTGGTAACGGCT  >  NC_000913/1808992‑1809171
                                                                                          |                                                                                         
tCGTTGTACTTCTGTTTGCGCTGGCTCAGCCCCGCCATAAACAGTGGAGTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGctt                                                                                         <  2:773810/93‑1 (MQ=255)
  gTTGTACTGCTGTTTGCGCTGGCTCAGACCCGCCATAAACAGTGGAGTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGctttt                                                                                       >  1:849651/1‑93 (MQ=255)
   ttGTTCTGCTGTTCGCGCTTCCTCAACCCCGCCATAAACAGTGGAGTCTGGCGAAAAAAGTTCTGGTGTGTCTGGTGATGGGTGTTGCTTTTg                                                                                      <  1:720963/93‑1 (MQ=255)
   ttGTACTGCTGTTTGCGCTGGCTCAGACCCGCCATAAACAGTGGAGTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGTTTTTg                                                                                      >  1:1081790/1‑93 (MQ=255)
   ttGTACTGCTGTTTGCGCTGGCTCAGACCCGCCATAAACAGTGGAGTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGTTTTTg                                                                                      >  1:93200/1‑93 (MQ=255)
   ttGTACTGCTGTTTGCGCTGGCTCAGACCCGCCATAAACAGTGGAGTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTg                                                                                      >  1:310215/1‑93 (MQ=255)
   ttGTACTGCTGTTTGCGCTGGCTCAGACCCGCCATAAACAGTGGAGTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTg                                                                                      >  2:5582041/1‑93 (MQ=255)
   ttGTACTGCTGTTTGCGCTGGCTCAGACCCGCCATAAACAGTGGAGTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTg                                                                                      >  2:5962542/1‑93 (MQ=255)
   ttGTACTGCTGTTTGCGCTGCCTCAGCCCCGCCATAAACAGTGGAGTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTg                                                                                      <  2:5750958/93‑1 (MQ=255)
   ttGTACTGCTGTTTGCGCGGGCTCAGCCCCGCCATAAACAGTGGAGTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTg                                                                                      <  2:5125450/93‑1 (MQ=255)
    tGTACTGCTGTTTGCGCTGGCTCAGACCCGCCATAAACAGTGGAGTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTgg                                                                                     >  1:6126469/1‑93 (MQ=255)
      tcgtgctgTTGGCGGGGGCTCACCCCCGCCATAAACAGTGGAGTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGcc                                                                                   <  2:107839/90‑1 (MQ=255)
      tcctgctgTTGGCGGTCGCCCAGCCCCGCCATAAACAGTGGAGTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGcc                                                                                   <  2:2439328/91‑1 (MQ=255)
       tctgctgTTTGCGCTGGCTCAGCCCCGCCATAAACAGTGGAGTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGCCt                                                                                  <  1:5283180/92‑1 (MQ=255)
       tctgctgTTGGCGCTGCCTCAGCCCCGCCATAAACAGTGGAGTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGCCt                                                                                  <  2:1296985/92‑1 (MQ=255)
       aCTGTTGTTGGCGCGGACTCAGCCCCGCCATAAACAGTGGAGTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGTTTTTGGCCt                                                                                  <  1:320826/93‑1 (MQ=255)
       aCTGCTGTTTGCGCTGTCTCAGACCCGCCATAAACAGTGGAGTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGCCt                                                                                  <  1:6160676/93‑1 (MQ=255)
       aCTGCTGTTTGCGCTGGCTCAGACCCGCCATAAACAGTGGAGTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTATgtcaa                                                                                  >  1:3715339/1‑89 (MQ=255)
       aCTGCTGTTTGCGCTGGCTCAGACCCGCCATAAACAGTGGAGCCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGCCt                                                                                  >  2:4670281/1‑93 (MQ=255)
        ctggtgTTGGCGGTGCCTGACCCCCGCCATAAACAGTGGAGTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGTTTTTGGCCtt                                                                                 <  2:6181297/93‑1 (MQ=255)
        ctgctgTTTGCGCTGGCTCAGACCCGCCATAAACAGTGGAGTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGCCtt                                                                                 >  1:1022915/1‑93 (MQ=255)
          tctgTTTGCGCTGGCTCAGCCCCGCCATAAACAGTGGAGTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGCCTTGc                                                                               <  1:5587894/92‑1 (MQ=255)
          gctgTTTGCGCTGTCTCACCCCCGCCATAAACAGTGGAGTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGCCTTGc                                                                               <  1:2035611/93‑1 (MQ=255)
          gctgTTTGCGCTGGCTCAGACCCGCCATAAACAGTGGAGTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGTTTTTGGCCTTGc                                                                               >  1:2854221/1‑93 (MQ=255)
           ctgTTTGCGCTGCCTCAGCCCCGCCATAAACAGTGGAGTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGCCTTGcc                                                                              <  2:1187689/93‑1 (MQ=255)
               ttGCGCGGCCTCGGCCCCGCCATAAACAGTGGAGTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGCCTTGCCCTGc                                                                          <  2:6102766/93‑1 (MQ=255)
                  cgcTGGCTCAGACCCGCCATAAACAGTGGAGTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGCCTTGCCCTGCATa                                                                       >  1:3046653/1‑93 (MQ=255)
                  cgcTGGCTCAGACCCGCCATAAACAGTGGAGTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGCCTTGCCCTGCATa                                                                       >  1:1760186/1‑93 (MQ=255)
                   gcTGCCTCAGACCCGCCATAAACAGTGGAGTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGCCTTGCCCTGCATAc                                                                      <  1:421286/93‑1 (MQ=255)
                    cTGGCTCAGACCCGCCATAAACAGTGGAGTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGCCTTGCCCTGCATAcc                                                                     >  1:5032731/1‑93 (MQ=255)
                     tGGCTCAGACCCGCCATAAACAGTGGAGTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGCCTTGCCCTGCATACCa                                                                    <  1:6410010/93‑1 (MQ=255)
                      ggCTCAGCCCCGCCATAAACAGTGGAGTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGCCTTGCCCTGCATACCAt                                                                   <  2:3152863/93‑1 (MQ=255)
                       ccTCAGACCCGCCATAAACAGTGGAGTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGCCTTGCCCTGCATACCAtt                                                                  <  2:3483660/92‑1 (MQ=255)
                         tCAGACCCGCCATAAACAGTGGAGTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGCCTTGCCCTGCATACCATTTa                                                                <  1:2515730/93‑1 (MQ=255)
                          cacccccGCCATAAACAGTGGAGTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGTTTTTGGCCTTGCCCTGCATACCATTTAt                                                               <  1:4574522/89‑1 (MQ=255)
                           cgccccGCCATAAACAGTGGAGTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGCCTTGCCCTGCATACCATTTATg                                                              <  2:3258217/90‑1 (MQ=255)
                           aGACCCGCCATAAACAGTGGAGTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGCCTTGCCCTGCATACCATTTATg                                                              >  2:219676/1‑93 (MQ=255)
                           aGACCCGCCATAAACAGTGGAGTCTGGCGAAAAAAGTGCTGGTGGGGGTGGTGATGGGTGTGGGTTTTGGTCTTGCCCTGCATACCATTTATg                                                              >  2:4857125/1‑93 (MQ=255)
                            gACCCGCCATAAACAGTGGAGTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGTTTTTGGCCTTGCCCTGCATACCATTTATgg                                                             >  2:2114638/1‑93 (MQ=255)
                            gACCCGCCATAAACAGTGGAGTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGCCTTGCCCTGCATACCATTTATgg                                                             >  2:1426696/1‑93 (MQ=255)
                             aCCCGCCATAAACAGTGGAGTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGCCTTGCCCTGCATACCATTTATGGt                                                            >  2:3012945/1‑93 (MQ=255)
                               ccGCCATAAACAGTGGAGTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGTTTTTGGCCTTGCCCTGCATACCATTTATGGTTc                                                          <  1:6000973/93‑1 (MQ=255)
                               ccGCCATAAACAGTGGAGTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGTTTTTGGCCTTGCCCTGCATACCATTCATGGTTc                                                          <  1:6312676/93‑1 (MQ=255)
                               ccGCCATAAACAGTGGAGTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGCCTTGCCCTGCATACCATTTATGGttt                                                          >  1:1907452/1‑92 (MQ=255)
                               ccGCCATAAACAGTGGAGTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGCCTTGCCCTGCATACCATTTATGGTTc                                                          <  2:784717/93‑1 (MQ=255)
                               ccGCCATAAACAGTGGAGTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGCCTTGCCCTGCATACCATTTATGGTTc                                                          >  2:424463/1‑93 (MQ=255)
                               ccGCCATAAACAGTGGAGTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGCCTTGCCCTGCATACCATTTATGGTTc                                                          <  2:706773/93‑1 (MQ=255)
                                 gCCATAAACAGTGGAGTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGCCTTGCCCTGCATACCATTTATGGTTCTg                                                        <  2:3700787/93‑1 (MQ=255)
                                 gCCATAAACAGTGGAGTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGCCTTGCCCTGCATACCATTTATGGTTCTg                                                        >  2:1004794/1‑93 (MQ=255)
                                    aTAAACAGTGGAGTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGCCTTGCCCTGCATACCATTTATGGTTCTGACa                                                     <  1:1694933/93‑1 (MQ=255)
                                      aaaCAGTGGAGTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGTTTTTGGCCTTGCCCTGCATACCATTTATGGTTCTGACAGc                                                   <  2:4349451/93‑1 (MQ=255)
                                      aaaCAGTGGAGTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGCCTTGCCCTGCATACCATTTATGGTTCTGACAGc                                                   <  1:3516098/93‑1 (MQ=255)
                                       aaCAGTGGAGTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGATTTTTGGCCTTGCCCTGCATACCATTTATGGTTCTGACAGcc                                                  <  2:3567555/93‑1 (MQ=255)
                                        aCAGTGGAGTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGCCTTGCCCTGCATACCATTTATTGTTCTGACAGCCa                                                 >  1:1266919/1‑93 (MQ=255)
                                        aCAGTGGAGTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGCCTTGCCCTGCATACCATTTATGGTTCTGACAGCCa                                                 <  2:747464/93‑1 (MQ=255)
                                         cAGTGGAGTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGTTTTTGGCCTTGCCCTGCATACCATTTATGGTTCTGACAGCCAg                                                <  1:648266/93‑1 (MQ=255)
                                          aGTGGAGTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGCCTTGCCCTGCATACCATTTATGGTTCTGACAGCCAgg                                               <  2:1707443/93‑1 (MQ=255)
                                          aGTGGAGTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGCCTTGCCCTGCATACCATTTATGGTTCTGACAGCCAgg                                               <  1:5957051/93‑1 (MQ=255)
                                           gTGGAGTCTGGCGAAAAAAGTGCTGGTGGGGCTGGTGATGGGTGTGGCTTTTGGCCTTGCCCTGCATACCATTTATGGTTCTGAGAGCCAGGt                                              >  1:384235/1‑93 (MQ=255)
                                            tGGAGTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGCCTTGCCCTGCATACCATTTATGGTTCTGACAGCCAGGTa                                             <  1:217186/93‑1 (MQ=255)
                                            tGGAGTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGCCTTGCCCTGCATACCATTTATGGTTCTGACAGCCAGGTa                                             >  2:658281/1‑93 (MQ=255)
                                             ggAGTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGCCTTGCCCTGCATACCATTTATGGTTCTGACAGCCAGGTAc                                            >  1:3059060/1‑93 (MQ=255)
                                              gAGTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGCCTTGCCCTGCATACCATTTATGGTTCTGACAGCCAGGTACt                                           >  2:1918793/1‑93 (MQ=255)
                                                gTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGCCTTGCCCTGCATACCATTTATGGTTCTGACAGCCAGGTACTTa                                         <  2:5384414/93‑1 (MQ=255)
                                                 tCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGCCTTGCCCTGCATACCATTTATGGTTCTGACAGCCAGGTACTTaa                                        <  1:2968917/93‑1 (MQ=255)
                                                  cTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGCCTTGCCCTGCATACCATTTATGGTTCTGACAGCCAGGTACTTaaa                                       <  2:4030051/93‑1 (MQ=255)
                                                  cTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGCCTTGCCCTGCATACCATTTATGGTTCTGACAGCCAGGTACTTaaa                                       <  1:1830933/93‑1 (MQ=255)
                                                  cTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGCCTTGCCCTGCATACCATTTATGGTTCTGACAGCCAGGTACTTaaa                                       <  2:1546042/93‑1 (MQ=255)
                                                  cTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGCCTTGCCCTGCATACCATTTATGGTTCTGACAGCCAGGTACTTaaa                                       <  2:2750234/93‑1 (MQ=255)
                                                     gCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGCCTTGCCCTGCATACCATTTATGGTTCTGACAGCCAGGTACTTAAAGAt                                    <  2:2286279/93‑1 (MQ=255)
                                                     gCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGCCTTGCCCTGCATACCATTTATGGTTCTGACAGCCAGGTACTTAAAGAt                                    <  1:6210429/93‑1 (MQ=255)
                                                       gAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGCCTTGCCCTGCATACCATTTATGGTTCTGACAGCCAGGTACTTAAAGATTc                                  >  1:4877282/1‑93 (MQ=255)
                                                       gAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGCCTTGCCCTGCATACCATTTATGGTTCTGACAGCCAGGTACTTAAAGATTc                                  <  1:3885186/93‑1 (MQ=255)
                                                        aaaaaaGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGCCTTGCCCTGCATACCATTTATGGTTCTGACAGCCAGGTACTTAAAGATTCt                                 >  1:2634664/1‑93 (MQ=255)
                                                         aaaaaGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGCCTTGCCCTGCATACCATTTATGGTTCTGACAGCCAGGTACTTAAAGATTCTg                                >  1:5880545/1‑93 (MQ=255)
                                                         aaaaaGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGCCTTGCCCTGCATACCATTTATGGTTCTGACAGCCAGGTACTTAAAGATTCTg                                >  1:6080321/1‑93 (MQ=255)
                                                           aaaGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGCCTTGCCCTGCATACCATTTATGGTTCTGACAGCCAGGTACTTAAAGATTCTGTa                              <  2:518499/93‑1 (MQ=255)
                                                           aaaGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGCCTTGCCCTGCATACCATTTATGGTTCTGACAGCCAGGTACTTAAAGATTCTGTa                              <  2:2598490/93‑1 (MQ=255)
                                                            aaGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGCCTTGCCCTGCATACCATTTATGGTTCTGACAGCCAGGTACTTAAAGATTCTGTAc                             >  1:1449427/1‑93 (MQ=255)
                                                            aaGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGCCTTGCCCTGCATACCATTTATGGTTCTGACAGCCAGGTACTTAAAGATTCTGTAc                             >  1:546977/1‑93 (MQ=255)
                                                             aGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGCCTTGCCCTGCATACCATTTATGGTTCTGACAGCCAGGTACTTAAAGATTCTGTACa                            <  1:2870016/93‑1 (MQ=255)
                                                             aGTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGCCTTGCCCTGCATACCATTTATGGTTCTGACAGCCAGGTACTTAAAGATTCAGTACa                            <  1:1990787/93‑1 (MQ=255)
                                                              gTGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGCCTTGCCCTGCATACCATTTATGGTTCTGACAGCCAGGTACTTAAAGATTCTGTACAg                           >  2:4111441/1‑93 (MQ=255)
                                                               tGCTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGCCTTGCCCTGCATACCATTTATGGTTCTGACAGCCAGGTACTTAAAGATTCTGTACAGt                          <  2:623420/93‑1 (MQ=255)
                                                                gCTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGCCTTGCCCTGCATACCATTTATGGTTCTGACAGCCAGGTACTTAAAGATTCTGTACAGTg                         >  2:1220478/1‑93 (MQ=255)
                                                                 cTGGTGGGTCTGGTGATGGGTGTGGCTTTTGGCCTTGCCCTGCATACCATTTATGGTTCTGACAGCCAGGTACTTAAAGATTCTGTACAGTgg                        >  2:627444/1‑93 (MQ=255)
                                                                   ggtggGTCTGGTGATGGGTGTGGCTTTTGGCCTTGCCCTGCATACCATTTATGGTTCTGACAGCCAGGTACTTAAAGATTCTGTACAGTGGtt                      >  1:5091701/1‑93 (MQ=255)
                                                                      gggTCTGGTGATGGGTGTGGCTTTTGGCCTTGCCCTGCATACCATTTATGGTTCTGACAGCCAGGTACTTAAAGATTCTGTACAGTGGTTTaa                   >  2:439809/1‑93 (MQ=255)
                                                                       ggTCTGGTGATGGGTGTGGCTTTTGGCCTTGCCCTGCATACCATTTATGGTTCTGACAGCCAGGTACTTAAAGATTCTGTACAGTGGTTTAAc                  <  1:4025867/93‑1 (MQ=255)
                                                                          cTGGTGATGGGTGTGGCTTTTGGCCTTGCCCTGCATACCATTTATGGTTCTGACAGCCAGGTACTTAAAGATTCTGTACAGTGGTTTAACATc               <  1:3290083/93‑1 (MQ=255)
                                                                          cTGGTGATGGGTGTGGCTTTTGGCCTTGCCCTGCATACCATTTATGGTTCTGACAGCCAGGTACTTAAAGATTCTGTACAGTGGTTTAACATc               >  1:4142188/1‑93 (MQ=255)
                                                                             gTGATGGGTGTGGCTTTTGGCCTTGCCCTGCATACCATTTATGGTTCTGATAGCCAGGTACTTAAAGATTCTGTACAGTGGTTTAACATCGtt            >  1:964153/1‑93 (MQ=255)
                                                                              tGATGGGTGTGGCTTTTGGCCTTGCCCTGCATACCATTTATGGTTCTGACAGCCAGGTACTTAAAGATTCTGTACAGTGGTTTAACATCGTTg           >  2:5282383/1‑93 (MQ=255)
                                                                               gATGGGTGTGGCTTTTGGCCTTGCCCTGCATACCATTTATGGTTCTGACAGCCAGGTACTTAAAGATTCTGTACAGTGGTTTAACATCGTTgg          >  2:1352255/1‑93 (MQ=255)
                                                                               gATGGGTGTGGCTTTTGGCCTTGCCCTGCATACCATTTATGGTTCTGACAGCCAGGTACTTAAAGATTCTGTACAGTGGTTTAACATCGTTgg          <  1:1406701/93‑1 (MQ=255)
                                                                                 tGGGTGTGGCTTTTGGCCTTGCCCTGCATACCATTTATGGTTCTGACAGCCAGGTACTTAAAGATTCTGTACAGTGGTTTAACATCGTTGGTa        >  2:4976360/1‑93 (MQ=255)
                                                                                 tGGGTGTGGCTTTTGGCCTTGCCCTGCATACCATTTATGGTTCTGACAGCCAGGTACTTAAAGATTCTGTACAGTGGTTTAACATCGTTGGTa        <  2:1802612/93‑1 (MQ=255)
                                                                                  gggTGTGGCTTTTGGCCTTGCCCTGCATACCATTTATGGTTCTGACAGCCAGGTACTTAAAGATTCTGTACAGTGGTTTAACATCGTTGGTaa       <  1:1655470/93‑1 (MQ=255)
                                                                                      gtgGCTTTTGGCCTTGCCCTGCATACCATTTATGGTTCTGACAGCCAGGTACTTAAAGATTCTGTACAGTGGTTTAACATCGTTGGTAACGGc   <  2:4328259/93‑1 (MQ=255)
                                                                                       tgGCTTTTGGCCTTGCCCTGCATACCATTTATGGTTCTGACAGCCAGGTACTTAAAGATTCTGTACAGTGGTTTAACATCGTTGGTAACGGCt  >  1:5741523/1‑93 (MQ=255)
                                                                                          |                                                                                         
TCGTTGTACTGCTGTTTGCGCTGGCTCAGACCCGCCATAAACAGTGGAGTCTGGCGAAAAAAGTGCTGGTGGGTCTGGTGATGGGTGTGGTTTTTGGCCTTGCCCTGCATACCATTTATGGTTCTGACAGCCAGGTACTTAAAGATTCTGTACAGTGGTTTAACATCGTTGGTAACGGCT  >  NC_000913/1808992‑1809171

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