Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 715,849 A→G 67.6% L146L (TTA→TTG)  ybfP → lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913715,8490AG67.6% 137.7 / 94.3 114L146L (TTA→TTG) ybfPlipoprotein
Reads supporting (aligned to +/- strand):  ref base A (22/15);  new base G (42/35);  total (64/50)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.78e-01

Only 100 of 124 total aligned reads displayed.

TGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCGGACCTGAACATGTCTT  >  NC_000913/715751‑715948
                                                                                                  |                                                                                                   
tgccctggCAAACCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTac                                                                                                    <  2:926807/94‑1 (MQ=255)
 gCCAGGGCAAGCCCGGCACTCGACCGATAGTTAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTgcg                                                                                                   <  1:848358/100‑1 (MQ=255)
  ccAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTgcga                                                                                                  >  2:1802798/1‑100 (MQ=255)
  ccAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTacga                                                                                                  >  1:234083/1‑100 (MQ=255)
   cAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTgcgac                                                                                                 >  1:880592/1‑100 (MQ=255)
   cAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTgcgac                                                                                                 >  2:655948/1‑100 (MQ=255)
   cAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGATCGATGATGCATGTCGGGATTATTGAGTTacgac                                                                                                 >  2:2645645/1‑100 (MQ=255)
    aGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTgcgacg                                                                                                >  2:2016079/1‑100 (MQ=255)
       gCAAGCCCGGCACTCGACAGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTAAGTTACGACGCTg                                                                                             >  2:1206129/1‑100 (MQ=255)
           gCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAg                                                                                         >  1:623828/1‑100 (MQ=255)
             ccGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTc                                                                                       >  2:2116493/1‑100 (MQ=255)
              cGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCa                                                                                      >  2:460017/1‑100 (MQ=255)
              cGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCa                                                                                      >  1:1351098/1‑100 (MQ=255)
               ggCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAg                                                                                     <  1:1642313/100‑1 (MQ=255)
                 cACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGaa                                                                                   <  2:1312029/100‑1 (MQ=255)
                   cTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATACATGTCGGGATTATTGATTTACGACGCTGCCAGTCAGaaaa                                                                                 >  1:199595/1‑100 (MQ=255)
                      gACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTg                                                                              >  1:1632976/1‑100 (MQ=255)
                      gACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTg                                                                              <  1:1375684/100‑1 (MQ=255)
                      gACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTg                                                                              <  1:1832314/100‑1 (MQ=255)
                       aCCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGc                                                                             <  1:2362733/100‑1 (MQ=255)
                        ccGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGcc                                                                            >  1:2035455/1‑100 (MQ=255)
                        ccGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGcc                                                                            >  2:274161/1‑100 (MQ=255)
                          gATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAg                                                                          <  2:595595/100‑1 (MQ=255)
                          gATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAg                                                                          >  1:522841/1‑100 (MQ=255)
                          gATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAg                                                                          <  1:338462/100‑1 (MQ=255)
                            tAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCa                                                                        >  1:1112319/1‑100 (MQ=255)
                              gTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAg                                                                      >  2:1304368/1‑100 (MQ=255)
                                cAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGaa                                                                    >  2:2290880/1‑100 (MQ=255)
                                cAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGaa                                                                    >  2:1869397/1‑100 (MQ=255)
                                 aaTGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAAt                                                                   >  1:1402055/1‑100 (MQ=255)
                                    gAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGa                                                                <  1:1157235/100‑1 (MQ=255)
                                    gAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGa                                                                >  1:1619405/1‑100 (MQ=255)
                                      ggATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGaaa                                                              <  1:2229066/100‑1 (MQ=255)
                                      ggATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGaaa                                                              >  1:1583306/1‑100 (MQ=255)
                                      ggATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGaaa                                                              <  1:254083/100‑1 (MQ=255)
                                        aTTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGc                                                            <  2:1637537/100‑1 (MQ=255)
                                            tGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATa                                                        <  1:660645/100‑1 (MQ=255)
                                             gTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATaa                                                       >  1:970579/1‑100 (MQ=255)
                                             gTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATaa                                                       >  2:1658682/1‑100 (MQ=255)
                                             gTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATaa                                                       >  1:2020063/1‑100 (MQ=255)
                                             gTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATaa                                                       >  2:346496/1‑100 (MQ=255)
                                              tCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGACAGTCAGAAAATTGCCAGCAAGAATTGAACGCGATaaa                                                      >  2:1691921/1‑100 (MQ=255)
                                               cAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAAt                                                     <  1:1524692/100‑1 (MQ=255)
                                               cAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAAt                                                     >  1:1535474/1‑100 (MQ=255)
                                               cAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAAt                                                     >  1:141654/1‑100 (MQ=255)
                                               cAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAAt                                                     >  1:1670733/1‑100 (MQ=255)
                                               cAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAAt                                                     <  2:654732/100‑1 (MQ=255)
                                               cAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAATCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAAt                                                     >  1:142173/1‑100 (MQ=255)
                                               cAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAAt                                                     >  1:2639963/1‑100 (MQ=255)
                                                aGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATa                                                    <  1:9863/100‑1 (MQ=255)
                                                 gCGAATGTGGATTTTTTGGCTCCATGATGCATGTCGGGATTATTGAGTGGCGACGCTGCAAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATa                                                    <  1:1208093/99‑1 (MQ=255)
                                                 gCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATac                                                   <  1:1129135/100‑1 (MQ=255)
                                                   gAAAGTGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATa                                                    <  1:1335365/97‑1 (MQ=255)
                                                   gAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATacac                                                 >  2:278119/1‑100 (MQ=255)
                                                   gAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATacac                                                 <  2:396988/100‑1 (MQ=255)
                                                     aaGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACAt                                               <  1:1805654/100‑1 (MQ=255)
                                                       gCGGATTTTTTGGCTCGATGATGCATGTCGGGTTTATTGAGTTACTACGCTGCCAGTCAGAAAATTGCAAGCAAGAATTGAAAGCGATAAATACACATTa                                             >  1:308919/1‑100 (MQ=255)
                                                       gCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTa                                             >  1:1529475/1‑100 (MQ=255)
                                                          ggatttttGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTaaaa                                          >  2:1373333/4‑100 (MQ=255)
                                                          gATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTaaaa                                          >  2:1383059/1‑100 (MQ=255)
                                                            ttttttGGCTCGATGATGCATGTCGGGATTATTGTGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATa                                        >  1:2318666/1‑100 (MQ=255)
                                                             tttttGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAg                                       <  2:903218/100‑1 (MQ=255)
                                                             tttttGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAg                                       <  2:2278359/100‑1 (MQ=255)
                                                              ttttGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGt                                      <  1:2645676/100‑1 (MQ=255)
                                                               tttGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTg                                     <  1:1302812/100‑1 (MQ=255)
                                                                ttGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCTAGAATTGAAAGCGATAAATACACATTAAAATAGTGa                                    <  1:1869414/100‑1 (MQ=255)
                                                                 tGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAg                                   >  1:2506479/1‑100 (MQ=255)
                                                                  ggCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGt                                  >  1:2365076/1‑100 (MQ=255)
                                                                    cTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTcc                                >  1:1740849/1‑100 (MQ=255)
                                                                    cTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTcc                                <  1:2127221/100‑1 (MQ=255)
                                                                          gatgCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAAt                          >  1:2030331/1‑100 (MQ=255)
                                                                           atgCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATg                         >  1:2282211/1‑100 (MQ=255)
                                                                           atgCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATg                         <  2:1147484/100‑1 (MQ=255)
                                                                            tgCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGa                        >  2:1517268/1‑100 (MQ=255)
                                                                              cATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGaaa                      <  1:549810/100‑1 (MQ=255)
                                                                                tGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAAcc                    >  1:2355158/1‑100 (MQ=255)
                                                                                tGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGaaaac                    >  2:1570600/1‑98 (MQ=255)
                                                                                 gTCGGGGTTATTAAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGAAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGACAGGGAATGAAAcc                    >  2:1008742/1‑99 (MQ=255)
                                                                                 gTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAAccc                   >  2:2536177/1‑100 (MQ=255)
                                                                                  tCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCg                  <  2:2006877/100‑1 (MQ=255)
                                                                                   cGGGATTGTTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCgg                 >  1:1392122/1‑100 (MQ=255)
                                                                                   cGGGATTATTGGGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATAGACATTAAAATAGTGAGTCCGAGAATGAAACCCgg                 >  2:1395487/1‑100 (MQ=255)
                                                                                   cGGGATTATTGAGTTGCGACGCTGCAAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCgg                 >  2:568901/1‑100 (MQ=255)
                                                                                    gggATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCgg                 >  2:1855490/1‑99 (MQ=255)
                                                                                    gggATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCGGa                <  2:1100728/100‑1 (MQ=255)
                                                                                    gggATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGACAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCg                  >  2:736194/1‑98 (MQ=255)
                                                                                    gggATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCGGa                <  2:20077/100‑1 (MQ=255)
                                                                                     ggATTATTGAGTTGCGACGCTGCCAGTCAGAAAGTTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCGGAc               >  1:1138424/1‑100 (MQ=255)
                                                                                     ggATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCGGa                >  1:579235/1‑99 (MQ=255)
                                                                                      gATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCGGAcc              <  1:2536207/100‑1 (MQ=255)
                                                                                      gATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCGGAcc              <  2:1351329/100‑1 (MQ=255)
                                                                                      gATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCGGAcc              <  2:2439443/100‑1 (MQ=255)
                                                                                        ttattGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCGGACCTg            >  1:2378581/1‑100 (MQ=255)
                                                                                           ttGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCGGACCTGAAc         <  2:409234/100‑1 (MQ=255)
                                                                                           ttGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGTAATGAAACCCGGACCTGAAc         <  2:1619389/100‑1 (MQ=255)
                                                                                           ttGAGTTACGACGCTACCACTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATGAAAACAGTGAGTCCAGGAATGAAACCCGGACCAGAAc         >  2:1764797/1‑100 (MQ=255)
                                                                                               gTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCGGACCTGAACATGt     >  1:2528751/1‑100 (MQ=255)
                                                                                               gTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCGGACCTGAACATGt     >  2:2210937/1‑100 (MQ=255)
                                                                                                ttACGACGCTACCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAATCTGGACATGAACATGTc    <  2:199593/100‑1 (MQ=255)
                                                                                                  acgacgCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCGGACCTGAACATGTCtt  <  2:955341/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
TGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCGGACCTGAACATGTCTT  >  NC_000913/715751‑715948

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑