Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,127,469 G→A G254G (GGC→GGT b2981 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,127,4690GA99.4% 1073.9 / NA 311G254G (GGC→GGTb2981hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (169/140);  minor base C (0/2);  total (169/142)
Fisher's exact test for biased strand distribution p-value = 2.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 315 total aligned reads displayed.

TTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGCAT  >  NC_000913/3127380‑3127558
                                                                                         |                                                                                         
ttCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCAccc                                                                                        >  2:7241703/1‑93 (MQ=255)
ttCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCAccc                                                                                        >  2:2888864/1‑93 (MQ=255)
ttCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCAccc                                                                                        >  1:803520/1‑93 (MQ=255)
 tCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCa                                                                                       <  2:4159779/93‑1 (MQ=255)
 tCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCa                                                                                       <  1:7249279/93‑1 (MQ=255)
  cGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCAt                                                                                      <  2:3509726/93‑1 (MQ=255)
  cGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCAt                                                                                      <  2:2129444/93‑1 (MQ=255)
  cGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCAt                                                                                      <  1:5267115/93‑1 (MQ=255)
   ggCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATg                                                                                     <  2:5399692/93‑1 (MQ=255)
    gCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATgg                                                                                    >  2:6904871/1‑93 (MQ=255)
      aaGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTc                                                                                  <  1:7763468/93‑1 (MQ=255)
      aaGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTc                                                                                  >  1:5941779/1‑93 (MQ=255)
      aaGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTc                                                                                  >  1:5132995/1‑93 (MQ=255)
       aGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTcg                                                                                 >  2:7865828/1‑93 (MQ=255)
       aGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTcg                                                                                 <  2:7625408/93‑1 (MQ=255)
        gCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTcgc                                                                                >  1:305649/1‑93 (MQ=255)
        gCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTcgc                                                                                >  2:5345798/1‑93 (MQ=255)
        gCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCAATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTcgc                                                                                >  1:3650472/1‑93 (MQ=255)
         ccGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTcgcg                                                                               >  1:6917081/1‑93 (MQ=255)
          cGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGa                                                                              >  2:3086476/1‑93 (MQ=255)
          cGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGa                                                                              >  2:1102054/1‑93 (MQ=255)
           gACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGaa                                                                             >  1:3791037/1‑93 (MQ=255)
           gACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGaa                                                                             >  1:7344467/1‑93 (MQ=255)
            aCAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAAt                                                                            >  1:7151272/1‑93 (MQ=255)
            aCAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAAt                                                                            >  2:5606693/1‑93 (MQ=255)
            aCAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAAt                                                                            >  1:4713707/1‑93 (MQ=255)
             cAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAAtt                                                                           >  1:2029566/1‑93 (MQ=255)
             cAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAAtt                                                                           <  1:5389556/93‑1 (MQ=255)
              aaTAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTc                                                                          >  1:4393146/1‑93 (MQ=255)
               aTAAACGCGCAGGGGGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCa                                                                         >  2:3818794/1‑93 (MQ=255)
                tAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCAt                                                                        <  2:699765/93‑1 (MQ=255)
                tAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCAt                                                                        <  1:701933/93‑1 (MQ=255)
                 aaaCGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATg                                                                       >  1:2500596/1‑93 (MQ=255)
                   aCGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGtt                                                                     <  2:5252058/93‑1 (MQ=255)
                   aCGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATCGTCGCGAATTCATGtt                                                                     <  2:5220203/93‑1 (MQ=255)
                          gggTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTaaa                                                              <  1:288422/93‑1 (MQ=255)
                          gggTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTaaa                                                              >  1:430475/1‑93 (MQ=255)
                          gggTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTaaa                                                              >  2:4177162/1‑93 (MQ=255)
                           ggTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAAc                                                             <  2:1262816/93‑1 (MQ=255)
                             tGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACgg                                                           <  1:1404862/93‑1 (MQ=255)
                             tGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACgg                                                           >  1:4518778/1‑93 (MQ=255)
                              gCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGc                                                          <  1:7081299/93‑1 (MQ=255)
                                 gAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGcaca                                                       <  1:2343274/93‑1 (MQ=255)
                                 gAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGcaca                                                       >  1:1285134/1‑93 (MQ=255)
                                 gAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGcaca                                                       >  1:3222964/1‑93 (MQ=255)
                                 gAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGcaca                                                       >  2:2590708/1‑93 (MQ=255)
                                 gAATCAATCTCAGCGATATAGATCAGATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGcaca                                                       <  2:5471974/93‑1 (MQ=255)
                                  aatcaatcTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACaa                                                      >  1:5386003/1‑93 (MQ=255)
                                  aatcaatcTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACaa                                                      <  2:3042876/93‑1 (MQ=255)
                                       atcTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATccc                                                 <  2:42838/93‑1 (MQ=255)
                                         ctcAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCa                                               <  1:3963864/93‑1 (MQ=255)
                                         ctcAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCa                                               <  1:6287935/93‑1 (MQ=255)
                                          tcAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCat                                              >  1:6025214/1‑92 (MQ=255)
                                          tcAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAg                                              >  1:5702783/1‑93 (MQ=255)
                                          tcAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAg                                              >  2:7414411/1‑93 (MQ=255)
                                            aGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGtt                                            <  2:6471783/93‑1 (MQ=255)
                                               gATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTtgtg                                         >  2:7500548/1‑93 (MQ=255)
                                                atatAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGATCGGTAAACGGCACAAATCCCCAGTTGtgc                                        <  1:3478333/93‑1 (MQ=255)
                                                 tataGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGtgct                                       >  1:3947386/1‑93 (MQ=255)
                                                 tataGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGtgct                                       <  2:1492257/93‑1 (MQ=255)
                                                    aGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTgcc                                    >  1:1078585/1‑93 (MQ=255)
                                                     gATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTgcca                                   >  1:4182131/1‑93 (MQ=255)
                                                     gATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTgcca                                   >  1:4666951/1‑93 (MQ=255)
                                                     gATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTgcca                                   <  2:7322974/93‑1 (MQ=255)
                                                     gATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTgcca                                   >  1:2637360/1‑93 (MQ=255)
                                                      aTCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTgccag                                  >  1:7513008/1‑93 (MQ=255)
                                                      aTCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTgccag                                  <  1:7667247/93‑1 (MQ=255)
                                                       tCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTgccagc                                 >  1:5516928/1‑93 (MQ=255)
                                                        cATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTgccagcc                                >  2:4679157/1‑93 (MQ=255)
                                                        cATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTgccagcc                                >  1:1666623/1‑93 (MQ=255)
                                                         atatCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTgccagcca                               <  2:5743012/93‑1 (MQ=255)
                                                         atatCCACCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTgccagcca                               <  1:6676919/93‑1 (MQ=255)
                                                          tatCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGccagccag                              >  1:1931534/1‑93 (MQ=255)
                                                          tatCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGGTGGCCTGcccgccag                              >  2:2012163/1‑93 (MQ=255)
                                                           atCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCaga                             >  2:1644783/1‑93 (MQ=255)
                                                           atCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCaga                             <  2:3564160/93‑1 (MQ=255)
                                                           atCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCaga                             <  2:4827395/93‑1 (MQ=255)
                                                            tCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCagag                            >  2:1755593/1‑93 (MQ=255)
                                                            tCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCagag                            <  1:6491362/93‑1 (MQ=255)
                                                            tCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCagag                            >  1:1177781/1‑93 (MQ=255)
                                                              gTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGtt                          <  2:2327145/93‑1 (MQ=255)
                                                               tCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTg                         >  1:2147461/1‑93 (MQ=255)
                                                                    cACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAAt                    >  2:2451657/1‑93 (MQ=255)
                                                                     aCCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAAtc                   <  2:5899206/93‑1 (MQ=255)
                                                                     aCCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAAtc                   <  2:6721938/93‑1 (MQ=255)
                                                                        aGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCa                <  1:3582423/93‑1 (MQ=255)
                                                                        aGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCa                <  1:6122836/93‑1 (MQ=255)
                                                                          atatTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCAAGTTGAGCTGCCAGCCAGAGTTGAAAACATCACg              >  1:1730317/1‑93 (MQ=255)
                                                                           tatTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGc             >  2:3411144/1‑93 (MQ=255)
                                                                            atTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCa            >  2:5587190/1‑93 (MQ=255)
                                                                            atTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCa            >  1:505968/1‑93 (MQ=255)
                                                                             tttAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCGGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAg           <  2:3222962/93‑1 (MQ=255)
                                                                             tttAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAg           >  2:1796689/1‑93 (MQ=255)
                                                                             tttAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAg           <  2:2683878/93‑1 (MQ=255)
                                                                              ttAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAgg          <  2:190597/93‑1 (MQ=255)
                                                                              ttAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAgg          <  1:154213/93‑1 (MQ=255)
                                                                               tAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGa         <  1:300157/93‑1 (MQ=255)
                                                                                aaGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGAt        <  2:6139284/93‑1 (MQ=255)
                                                                                  gTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCt      >  2:1956520/1‑93 (MQ=255)
                                                                                    tGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGc    >  2:5503075/1‑93 (MQ=255)
                                                                                         |                                                                                         
TTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGCAT  >  NC_000913/3127380‑3127558

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