Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,647,298 T→C T1004A (ACG→GCG)  b2520 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,647,2980TC99.3% 949.8 / ‑3.6 299T1004A (ACG→GCG) b2520hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base C (158/139);  total (160/139)
Fisher's exact test for biased strand distribution p-value = 5.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

Only 100 of 314 total aligned reads displayed.

CACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCG  >  NC_000913/2647208‑2647387
                                                                                          |                                                                                         
cACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGccc                                                                                         >  2:650152/1‑93 (MQ=255)
cACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGccc                                                                                         >  2:6968769/1‑93 (MQ=255)
cACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGccc                                                                                         >  2:7855644/1‑93 (MQ=255)
 aCTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCa                                                                                        >  2:5870443/1‑93 (MQ=255)
  cTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAg                                                                                       >  2:9740158/1‑93 (MQ=255)
  cTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAg                                                                                       >  1:8395677/1‑93 (MQ=255)
   tCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAgg                                                                                      >  1:3994297/1‑93 (MQ=255)
     gCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTATCGAAGTCATCTGCCGCCCAGGCt                                                                                    >  2:6927904/1‑93 (MQ=255)
     gCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCt                                                                                    <  1:2537560/93‑1 (MQ=255)
     gCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCt                                                                                    >  2:1916019/1‑93 (MQ=255)
     gCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCt                                                                                    <  1:2169820/93‑1 (MQ=255)
       cATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTg                                                                                  <  2:8491376/93‑1 (MQ=255)
       cATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTg                                                                                  <  2:628496/93‑1 (MQ=255)
       cATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTg                                                                                  >  2:2944231/1‑93 (MQ=255)
       cATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTg                                                                                  >  2:3792951/1‑93 (MQ=255)
       cATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTg                                                                                  <  2:4064519/93‑1 (MQ=255)
        aTAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTgc                                                                                 >  1:6981484/1‑93 (MQ=255)
        aTAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTgc                                                                                 <  2:5286794/93‑1 (MQ=255)
        aTAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTgc                                                                                 >  2:8925912/1‑93 (MQ=255)
           aaGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGcc                                                                              <  1:8845425/93‑1 (MQ=255)
           aaGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGcc                                                                              >  1:4687413/1‑93 (MQ=255)
               gcgGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGa                                                                          <  2:1323923/93‑1 (MQ=255)
               gcgGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGa                                                                          <  1:8546069/93‑1 (MQ=255)
                cgGCATGTTCAGTTCAGCATGCACCGGTGTGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGAc                                                                         <  2:7424705/93‑1 (MQ=255)
                 ggCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGAcg                                                                        <  1:5728465/93‑1 (MQ=255)
                  gCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGAcgc                                                                       >  2:9223184/1‑93 (MQ=255)
                     tGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCa                                                                    >  2:7954523/1‑93 (MQ=255)
                     tGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCa                                                                    >  1:4276742/1‑93 (MQ=255)
                     tGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCa                                                                    >  1:9545862/1‑93 (MQ=255)
                      gttcagttcaGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCaa                                                                   >  2:4741259/1‑93 (MQ=255)
                      gttcagttcaGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCaa                                                                   >  2:2229562/1‑93 (MQ=255)
                       ttcagttcagCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAAt                                                                  <  2:4853223/93‑1 (MQ=255)
                       ttcagttcagCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAAt                                                                  <  2:3622325/93‑1 (MQ=255)
                       ttcagttcagCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAAt                                                                  <  2:9472051/93‑1 (MQ=255)
                       ttcagttcagCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAAt                                                                  <  1:3804693/93‑1 (MQ=255)
                        tcagttcagCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAAtt                                                                 <  1:8099697/93‑1 (MQ=255)
                            ttcagCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCAcc                                                             >  1:5564254/1‑93 (MQ=255)
                             tcagCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCg                                                            >  1:4961165/1‑93 (MQ=255)
                              cagCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGt                                                           <  2:9830350/93‑1 (MQ=255)
                                gCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTa                                                         >  1:4709089/1‑93 (MQ=255)
                                 cAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTaa                                                        <  1:2769881/93‑1 (MQ=255)
                                   aTCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTaaaa                                                      >  1:8722509/1‑93 (MQ=255)
                                   aTCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTaaaa                                                      >  1:6626282/1‑93 (MQ=255)
                                    tCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCCATTCACCGTTAAAAt                                                     >  2:5487391/1‑93 (MQ=255)
                                    tCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAAt                                                     <  1:9174875/93‑1 (MQ=255)
                                    tCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAAt                                                     >  1:1782546/1‑93 (MQ=255)
                                    tCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAAt                                                     >  2:2659442/1‑93 (MQ=255)
                                       ccGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCTCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGc                                                  <  1:6285602/93‑1 (MQ=255)
                                       ccGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGc                                                  >  1:6057987/1‑93 (MQ=255)
                                         ggTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCGATGACGCGCAATTCACCGTTAAAATCGCCa                                                <  1:301249/93‑1 (MQ=255)
                                         ggTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCa                                                <  1:7353226/93‑1 (MQ=255)
                                         ggTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCa                                                <  2:6165854/93‑1 (MQ=255)
                                               gCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGc                                          <  1:3902235/93‑1 (MQ=255)
                                                 aaCTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAg                                        >  1:1621968/1‑93 (MQ=255)
                                                  aCTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAgt                                       >  2:8735560/1‑93 (MQ=255)
                                                   cTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAgtg                                      >  2:567563/1‑93 (MQ=255)
                                                    tATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAgtgt                                     >  2:4896840/1‑93 (MQ=255)
                                                    tATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAgtgt                                     >  2:1515024/1‑93 (MQ=255)
                                                     aTCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTa                                    >  2:5025956/1‑93 (MQ=255)
                                                      tCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTaa                                   <  2:3977679/93‑1 (MQ=255)
                                                       cACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAAc                                  >  2:3109307/1‑93 (MQ=255)
                                                        cctttactttCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAAcc                                 <  2:2067857/92‑1 (MQ=255)
                                                          tttactttCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGa                               >  2:10018700/1‑93 (MQ=255)
                                                           ttactttCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAg                              >  2:6645509/1‑93 (MQ=255)
                                                              ctttCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCt                           <  1:10014241/93‑1 (MQ=255)
                                                               tttCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCtt                          >  1:5870463/1‑93 (MQ=255)
                                                               tttCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCtt                          >  1:3253365/1‑93 (MQ=255)
                                                                    ttGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGccc                     >  2:2481262/1‑93 (MQ=255)
                                                                    ttGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGccc                     >  1:9424284/1‑93 (MQ=255)
                                                                    ttGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGccc                     >  2:5139371/1‑93 (MQ=255)
                                                                    ttGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGccc                     >  1:7757346/1‑93 (MQ=255)
                                                                    ttGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGccc                     >  2:8255263/1‑93 (MQ=255)
                                                                    ttGCGACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGccc                     <  1:9890539/93‑1 (MQ=255)
                                                                     tGCTTCCGAAGTCATCTGACGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTACCCGAGCCTTCGCCCt                    <  2:4761108/93‑1 (MQ=255)
                                                                      gCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTg                   <  2:768058/93‑1 (MQ=255)
                                                                      gCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTg                   >  1:4383997/1‑93 (MQ=255)
                                                                       cTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGt                  >  1:4883695/1‑93 (MQ=255)
                                                                        tACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGtt                 <  2:4560591/93‑1 (MQ=255)
                                                                        tACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGtt                 >  1:3510605/1‑93 (MQ=255)
                                                                         aCCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATTGCCAATCGGCAGTGTAACCGCGCCTTCGCCCTGTTc                >  1:7175338/1‑93 (MQ=255)
                                                                           cGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGt              <  1:1974442/93‑1 (MQ=255)
                                                                            gAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGtt             >  2:8064226/1‑93 (MQ=255)
                                                                            gAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGtt             >  1:8570560/1‑93 (MQ=255)
                                                                            gAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGtt             >  2:1184922/1‑93 (MQ=255)
                                                                             aaGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTg            >  2:4182629/1‑93 (MQ=255)
                                                                             aaGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTg            <  2:1450296/93‑1 (MQ=255)
                                                                              aGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTgg           >  1:4153955/1‑92 (MQ=255)
                                                                               gTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAg          >  1:6891830/1‑93 (MQ=255)
                                                                                tCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGc         <  1:1133128/93‑1 (MQ=255)
                                                                                tCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGc         >  1:6389992/1‑93 (MQ=255)
                                                                                 cATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCg        <  2:9299381/93‑1 (MQ=255)
                                                                                  aTCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGt       >  1:529009/1‑93 (MQ=255)
                                                                                  aTCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGt       <  1:9722820/93‑1 (MQ=255)
                                                                                  aTCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGt       <  1:6741516/93‑1 (MQ=255)
                                                                                  aTCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGt       >  2:7211949/1‑93 (MQ=255)
                                                                                     tGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTAc    <  2:1953571/93‑1 (MQ=255)
                                                                                     tGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTAc    <  2:3224022/93‑1 (MQ=255)
                                                                                     tGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTAc    <  2:3272993/93‑1 (MQ=255)
                                                                                     tGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGACCGTTAc    >  1:3195854/1‑93 (MQ=255)
                                                                                       ccGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCg  <  1:1232743/93‑1 (MQ=255)
                                                                                          |                                                                                         
CACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCG  >  NC_000913/2647208‑2647387

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