Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,649,291 +C coding (1017/4962 nt) b2520 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,649,2871.C98.6% 733.6 / 1.3 214coding (1021/4962 nt)b2520hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base C (105/106);  total (106/108)
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

Only 100 of 257 total aligned reads displayed.

CAGCGCCGGAAGTTTTAAATCGAGTAGCGTTGTCTGACCGTCTTTACGCGCCAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCTTGCTGGGCCGCGCCGTT  >  NC_000913/2649195‑2649374
                                                                                             |                                                                                       
cAGCGCCGGAAGTTTTAAATCGAGTAGCGTTGTCTGACCGTCTTTACGCGCCAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATg                                                                                          >  2:7163908/1‑93 (MQ=255)
 aGCGCCGGAAGTTTTAAATCGAGTAGCGTTGTCTGACCGTCTTTACGCGCCAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGc                                                                                        <  2:6675911/93‑1 (MQ=255)
  gcgcCGGAAGTTTTAAATCGAGTAGCGTTGTCTGACCGTCTTTACGCGCCAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGcc                                                                                       >  2:2599340/1‑93 (MQ=255)
  gcgcCGGAAGTTTTAAATCGAGTAGCGTTGTCTGACCGTCTTTACGCGCCAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGcc                                                                                       <  1:1117106/93‑1 (MQ=255)
  gcgcCGGAAGTTTTAAATCGAGTAGCGTTGTCTGACCGTCTTTACGCGCCAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGcc                                                                                       <  1:6358275/93‑1 (MQ=255)
     ccGGAAGTTTTAAATCGAGTAGCGTTGTTTGACCGTCTTTACGCGCCAACAGTACTGCCGCGTTTTTATCATTTTCCAGCTGCACATGccccc                                                                                    >  1:4106915/1‑92 (MQ=255)
     ccGGAAGTTTTAAATCGAGTAGCGTTGTCTGACCGTCTTTACGCGCCAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGccccc                                                                                    <  2:6740438/93‑2 (MQ=255)
       ggAAGTTTTAAATCGAGTAGCGTTGTCTGACCGTCTTTACGCGCCAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGccccctg                                                                                  <  1:4726772/93‑4 (MQ=255)
       ggAAGTTTTAAATCGAGTAGCGTTGTCTGACCGTCTTTACGCGCCAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGccccctg                                                                                  <  2:7108211/93‑4 (MQ=255)
       ggAAGTTTTAAATCGAGTAGCGTTGTCTGACCGTCTTTACGCGCCAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGccccctg                                                                                  >  1:5208497/1‑90 (MQ=255)
       ggAAGTTTTAAATCGAGTAGCGTTGTCTGACCGTCTTTACGCGCCAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGccccctg                                                                                  <  2:759793/93‑4 (MQ=255)
        gAAGTTTTAAATCGAGTAGCGTTGTCTGACCGTCTTTACGCGCCAACAGTACTGCCGCGTTTTTATCATTTTCCAGCTGCACATGccccctga                                                                                 <  2:8772489/93‑5 (MQ=255)
        gAAGTTTTAAATCGAGTAGCGTTGTCTGACCGTCTTTACGCGCCAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGccccctga                                                                                 >  1:5985061/1‑89 (MQ=255)
        gAAGTTTTAAATCGAGTAGCGTTGTCTGACCGTCTTTACGCGCCAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGccccctga                                                                                 <  2:2727359/93‑5 (MQ=255)
        gAAGTTTTAAATCGAGTAGCGTTGTCTGACCGTCTTTACGCGCCAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGccccctga                                                                                 <  1:518076/93‑5 (MQ=255)
           gTTTTAAATCGAGTAGCGTTGTCTGACCGTCTTTACGCGCCAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCCTGAGCg                                                                               >  2:6475320/1‑93 (MQ=255)
            ttttAAATCGAGTAGCGTTGTCTGACCGTCTTTACGCGCCAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCCTGAGCGt                                                                              <  1:1439805/93‑1 (MQ=255)
            ttttAAATCGAGTAGCGTTGTCTGACCGTCTTTACGCGCCAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCCTGAGCGt                                                                              <  1:4989480/93‑1 (MQ=255)
              ttAAATCGAGTAGCGTTGTCTGACCGTCTTTACGCGCCAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCCTGAGCGTCa                                                                            >  1:698880/1‑93 (MQ=255)
              ttAAATCGAGTAGCGTTGTCTGACCGTCTTTACGCGCCAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCCTGAGCGTCa                                                                            >  2:2918536/1‑93 (MQ=255)
               tAAATCGAGTAGCGTTGTCTGACCGTCTTTACGCGCCAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCCTGAGCGTCAc                                                                           <  1:1809934/93‑1 (MQ=255)
                 aaTCGAGTAGCGTTGTCTGACCGTCTTTACGCGCCAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCCTGAGCGTCACTg                                                                         >  1:5906060/1‑93 (MQ=255)
                   tCGAGTAGCGTTGTCTGACCGTCTTTACGCGCCAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCCTGAGCGTCACTGGt                                                                       <  2:1577791/93‑1 (MQ=255)
                    cGAGTAGCGTTGTCTGACCGTCTTTACGCGCCAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCCTGAGCGTCACTGGtt                                                                      >  1:4304355/1‑93 (MQ=255)
                     gAGTAGCGTTGTCTGACCGTCTTTACGCGCCAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTg                                                                     <  2:2381488/93‑1 (MQ=255)
                     gAGTAGCGTTGTCTGACCGTCTTTACGCGCCAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTg                                                                     >  2:4265608/1‑93 (MQ=255)
                     gAGTAGCGTTGTCTGACCGTCTTTACGCGCCAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTg                                                                     >  1:7100229/1‑93 (MQ=255)
                      aGTAGCGTTGTCTGACCGTCTTTACGCGCCAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTGc                                                                    >  1:53318/1‑93 (MQ=255)
                       gTAGCGTTGTCTGACCGTCTTTACGCGCCAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTGcc                                                                   >  1:2918763/1‑93 (MQ=255)
                       gTAGCGTTGTCTGACCGTCTTTACGCGCCAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTGcc                                                                   <  2:444968/93‑1 (MQ=255)
                        tAGCGTTGTCTGACCGTCTTTACGCGCCAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTGCCt                                                                  >  2:3071458/1‑93 (MQ=255)
                          gCGTTGTCTGACCGTCTTTACGCGCCAACAGTTATGCCGCGTTTTTATCATGTTCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGa                                                                >  1:7816694/1‑93 (MQ=255)
                          gCGTTGTCTGACCGTCTTTACGCGCCAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGa                                                                <  2:2180644/93‑1 (MQ=255)
                            gTTGTCTGACCGTCTTTACGCGCCAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGt                                                              <  1:792149/93‑1 (MQ=255)
                            gTTGTCTGACCGTCTTTACGCGCCAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGt                                                              <  2:866688/93‑1 (MQ=255)
                               gTCTGACCGTCTTTACGCGCCAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGTCag                                                           >  2:7753212/1‑93 (MQ=255)
                                tCTGACCGTCTTTACGCGCCAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGTCaga                                                          >  2:8099377/1‑93 (MQ=255)
                                  tGACCGTCTTTACGCGCCAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGt                                                        >  1:1329873/1‑93 (MQ=255)
                                   gACCGTCTTTACGCGCCAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTc                                                       >  1:4824159/1‑93 (MQ=255)
                                   gACCGTCTTTACGCGCCAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTc                                                       <  2:5672711/93‑1 (MQ=255)
                                    aCCGTCTTTACGCGCCAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCt                                                      >  1:8685846/1‑93 (MQ=255)
                                    aCCGTCTTTACGCGCCAACAGTAACGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCCTGAGCGTCACTGGCTGCCTGAGCCAGAGTCt                                                      >  1:7627076/1‑93 (MQ=255)
                                      cGTCTTTACGCGCCAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGc                                                    >  1:6107530/1‑93 (MQ=255)
                                        tCTTTACGCGCCAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGccc                                                  <  1:1733080/93‑1 (MQ=255)
                                        tCTTTACGCGCCAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGccc                                                  <  1:2372110/93‑1 (MQ=255)
                                        tCTTTACGCGCCAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGccc                                                  >  1:2501053/1‑93 (MQ=255)
                                          tttACGCGCCAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCtt                                                <  1:4643326/93‑1 (MQ=255)
                                           ttACGCGCCAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCttt                                               <  2:913449/93‑1 (MQ=255)
                                           ttACGCGCCAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCttt                                               >  1:1949390/1‑93 (MQ=255)
                                            tACGCGCCAACAGTAATGCCGCGATTTTATCATTTTCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTtc                                              <  2:4749352/93‑1 (MQ=255)
                                             aCGCGCCAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTtct                                             <  2:5059147/93‑1 (MQ=255)
                                             aCGCGCCAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTtct                                             <  2:6505034/93‑1 (MQ=255)
                                              cgcgCCAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTtctc                                            >  2:5157287/1‑93 (MQ=255)
                                               gcgcCAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCa                                           <  1:3776358/93‑1 (MQ=255)
                                                  ccAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCAttt                                        >  2:6028192/1‑93 (MQ=255)
                                                  ccAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCAttt                                        <  1:5749984/93‑1 (MQ=255)
                                                   cAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTta                                       >  2:4437121/1‑93 (MQ=255)
                                                     aCAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTtaat                                     <  2:1284254/93‑1 (MQ=255)
                                                     aCAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTtaat                                     <  1:8764458/93‑1 (MQ=255)
                                                        gTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATaaa                                  <  1:3546727/93‑1 (MQ=255)
                                                        gTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATaaa                                  <  2:6575297/93‑1 (MQ=255)
                                                         tAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCCTTTAATAAag                                 >  2:2072692/1‑93 (MQ=255)
                                                          aaTGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCCTGAGCTTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAaga                                <  2:7885090/93‑1 (MQ=255)
                                                          aaTGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAaga                                <  2:1620706/93‑1 (MQ=255)
                                                             gCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGAc                             <  2:3733337/93‑1 (MQ=255)
                                                              ccGCGTTTTTATCATTTTCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACt                            <  1:8448570/93‑1 (MQ=255)
                                                                gcgTTTTTATCATTTTCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTc                          <  1:7999675/93‑1 (MQ=255)
                                                                gcgTGTGTATCATTTTCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTc                          <  1:4588668/93‑1 (MQ=255)
                                                                 cgTTTTTATCATTTTCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCa                         >  2:4691841/1‑93 (MQ=255)
                                                                  gTTTTTATCATTTTCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCaa                        >  1:2325928/1‑93 (MQ=255)
                                                                  gTTTTTATCATTTTCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCaa                        >  2:1866101/1‑93 (MQ=255)
                                                                  gTTTTTATCATTTTCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCaa                        <  1:339498/93‑1 (MQ=255)
                                                                  gTTTTTATCATTTTCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCaa                        <  2:7137892/93‑1 (MQ=255)
                                                                   tttttATCATTTTCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAAt                       <  2:402529/93‑1 (MQ=255)
                                                                   tttttATCATTTTCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAAt                       <  1:7143620/93‑1 (MQ=255)
                                                                    ttttATCATTTTCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAAtt                      >  1:2663324/1‑93 (MQ=255)
                                                                     tttATCATTTTCCAGCTGCACATGCCCCCTGATCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTc                     <  1:1462010/93‑1 (MQ=255)
                                                                      ttATCATTTTCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTcc                    >  2:365668/1‑93 (MQ=255)
                                                                       tATCATTTTCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCt                   <  1:978647/93‑1 (MQ=255)
                                                                       tATCATTTTCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCt                   <  2:4193848/93‑1 (MQ=255)
                                                                        aTCATTTTCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCtt                  <  2:6537340/93‑1 (MQ=255)
                                                                        aTCATTTTCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCtt                  >  2:3397951/1‑93 (MQ=255)
                                                                          cATTTTCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCTTGc                <  1:3058672/93‑1 (MQ=255)
                                                                            ttttCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCTTGCTg              >  1:8069326/1‑93 (MQ=255)
                                                                            ttttCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCTTGCTg              <  2:2536951/93‑1 (MQ=255)
                                                                             tttCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCTTGCTgg             >  2:5277255/1‑93 (MQ=255)
                                                                              ttCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCTTGCTggg            <  1:6469101/93‑1 (MQ=255)
                                                                              ttCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCTTGCTggg            >  2:8419506/1‑93 (MQ=255)
                                                                               tCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCTTGCTGGGc           <  2:4044858/93‑1 (MQ=255)
                                                                               tCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCTTGCTGGGc           <  2:7262452/93‑1 (MQ=255)
                                                                               tCCAGCTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCTTGCTGGGc           >  2:618870/1‑93 (MQ=255)
                                                                                ccAGCTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCTTGCTGGGcc          <  1:7664036/93‑1 (MQ=255)
                                                                                ccAGCTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCTTGCTGGGcc          <  1:6083652/93‑1 (MQ=255)
                                                                                   gCTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCTTGCTGGGCcgcg       <  2:5465048/93‑1 (MQ=255)
                                                                                    cTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCTTGCTGGGCCgcgc      <  1:2216014/93‑1 (MQ=255)
                                                                                    cTGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCTTGCTGGGCCgcgc      >  2:1070365/1‑93 (MQ=255)
                                                                                     tGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCTTGCTGGGCCGCGcc     >  2:7455389/1‑93 (MQ=255)
                                                                                     tGCACATGCCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCTTGCTGGGCCGCGcc     <  1:2047812/93‑1 (MQ=255)
                                                                                      gaacaTGCCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCTTGCTGGGCCGCGCCg    <  2:6352113/91‑1 (MQ=255)
                                                                                        acaTGCCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCTTGCTGGGCCGCGCCGtt  <  1:5426446/93‑1 (MQ=255)
                                                                                             |                                                                                       
CAGCGCCGGAAGTTTTAAATCGAGTAGCGTTGTCTGACCGTCTTTACGCGCCAACAGTAATGCCGCGTTTTTATCATTTTCCAGCTGCACATGCCCCTGAGCGTCACTGGTTGCCTGAGTCAGAGTCTGCCCTTTCTCATTTAATAAAGAGACTTCAATTCCTTGCTGGGCCGCGCCGTT  >  NC_000913/2649195‑2649374

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