Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,944,114 (G)6→7 noncoding (411/2904 nt) rrlC → 23S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,944,1081.G99.2% 469.7 / ‑2.8 123noncoding (405/2904 nt)rrlC23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (63/59);  total (63/60)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.08e-01

Only 100 of 163 total aligned reads displayed.

CTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGG‑GGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACC  >  NC_000913/3944008‑3944159
                                                                                                     |                                                    
cTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAatat                                                       >  1:7256946/1‑101 (MQ=35)
cTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAatat                                                       <  1:7256510/101‑1 (MQ=38)
cTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAatat                                                       <  1:7256331/101‑1 (MQ=38)
cTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAatat                                                       <  2:7256508/101‑1 (MQ=38)
 tGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑g                                                     >  2:7256949/1‑101 (MQ=35)
 tGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑g                                                     >  1:7256950/1‑101 (MQ=35)
  ggAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACTTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑gg                                                    >  2:7256955/1‑101 (MQ=25)
  ggAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑gg                                                    <  1:7256357/101‑1 (MQ=38)
  ggAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑gg                                                    >  1:7256953/1‑101 (MQ=35)
  ggAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑gg                                                    >  1:7256954/1‑101 (MQ=35)
  ggAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑gg                                                    <  1:7256405/101‑1 (MQ=38)
  ggAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑g                                                     >  1:7256957/1‑100 (MQ=37)
   gAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑ggg                                                   >  1:7256958/1‑101 (MQ=35)
    aaaGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑ggg‑g                                                 <  1:7256586/101‑1 (MQ=38)
     aaGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑ggg‑gg                                                <  2:7256795/101‑1 (MQ=35)
      aGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑ggg‑ggg                                               >  2:7256959/1‑101 (MQ=35)
      aGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑ggg‑ggg                                               <  2:7256631/101‑1 (MQ=35)
       ggCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑ggg‑gggg                                              <  1:7256316/101‑2 (MQ=25)
       ggCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑ggg‑gggg                                              <  2:7256308/101‑2 (MQ=25)
         cgcgcgATACAGGGGGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑ggg‑ggggac                                            >  1:7256962/1‑98 (MQ=255)
          gcgcgATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑ggg‑gggg                                              >  2:7256963/1‑97 (MQ=39)
           cgcgATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑ggg‑ggggacca                                          <  2:7256177/101‑6 (MQ=255)
           cgcgATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑ggg‑ggggacca                                          >  2:7256964/1‑96 (MQ=255)
            gcgATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑ggg‑gg                                                <  1:7256464/94‑1 (MQ=37)
            gcgATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGTGGGACCAt                                          <  2:7255973/101‑1 (MQ=255)
             cgATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCAtc                                         >  2:7256969/1‑101 (MQ=255)
             cgATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCAtc                                         >  1:7256966/1‑101 (MQ=255)
             cgATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCAtc                                         >  1:7256968/1‑101 (MQ=255)
               aTACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCAtcct                                       >  1:7256971/1‑101 (MQ=255)
                tACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCAtcctc                                      <  1:7256560/101‑1 (MQ=255)
                 aCAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCAtcctcc                                     >  2:7256975/1‑101 (MQ=255)
                 aCAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCAtcctcc                                     <  1:7256408/101‑1 (MQ=255)
                 aCAGGGTGACAGCCCCGTACACAAAAATGCACATATTGAGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCAtcctcc                                     >  1:7256974/1‑101 (MQ=255)
                  cAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCa                                    >  2:7256976/1‑101 (MQ=255)
                  cAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCa                                    <  2:7256570/101‑1 (MQ=255)
                  cAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCa                                    >  1:7256977/1‑101 (MQ=255)
                   aGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCaa                                   >  2:7256979/1‑101 (MQ=255)
                   aGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCaa                                   >  1:7256981/1‑101 (MQ=255)
                   aGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑ggg‑gggga                                             <  2:7256124/90‑3 (MQ=38)
                   aGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑ggg‑gggg                                              >  1:7256982/1‑88 (MQ=255)
                     ggTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCAAgg                                 <  2:7256442/101‑1 (MQ=255)
                      gTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCAAGGc                                <  1:7256112/101‑1 (MQ=255)
                      gTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCAAGGc                                <  1:7256439/101‑1 (MQ=255)
                       tGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCAAGGCt                               <  2:7256585/101‑1 (MQ=255)
                        gACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCAAGGCTa                              >  2:7256985/1‑101 (MQ=255)
                         aCAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCAAGGCTaa                             <  2:7256288/101‑1 (MQ=255)
                         aCAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCAAGGCTaa                             <  1:7256424/101‑1 (MQ=255)
                         aCAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCAAGGCTaa                             >  1:7256986/1‑101 (MQ=255)
                          cAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCAAGGCTaaa                            <  2:7256441/101‑1 (MQ=255)
                           aGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCAAGGCTAAAt                           >  2:7256992/1‑101 (MQ=255)
                           aGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCAAGGCTAAAt                           >  2:7256991/1‑101 (MQ=255)
                           aGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCAAGGCTAAAt                           >  2:7256989/1‑101 (MQ=255)
                             ccccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCAAGGCTAAATAc                         >  1:7256995/1‑101 (MQ=255)
                             ccccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCAAGGCTAAATAc                         <  1:7256577/101‑1 (MQ=255)
                             ccccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCAAGGCTAAATAc                         >  1:7256994/1‑101 (MQ=255)
                              cccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCAAGGCTAAATACt                        <  2:7256674/101‑1 (MQ=255)
                                 gTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCAAGGCTAAATACTCct                     <  1:7255892/101‑1 (MQ=255)
                                 gTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCAAGGCTAAATACTCct                     <  2:7256353/101‑1 (MQ=255)
                                 ctACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCAAGGCTAAATACTCct                     >  1:7256998/2‑101 (MQ=255)
                                  tACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCAAGGCTAAATACTCctg                    <  1:7256817/101‑1 (MQ=255)
                                  tACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCAAGGCTAAATACTCctg                    <  1:7256435/101‑1 (MQ=255)
                                  tACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCAAGGCTAAATACTCctg                    <  2:7256543/101‑1 (MQ=255)
                                  tACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCAAGGCTAAATACTCctg                    >  2:7256999/1‑101 (MQ=255)
                                    cacaAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCAAGGCTAAATACTc                       >  1:7257002/1‑96 (MQ=255)
                                    cacaAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCAAGGCTAAATACTCctgac                  <  1:7256254/101‑1 (MQ=255)
                                     acaAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCAAGGCTAAATACTCctgact                 <  1:7256494/101‑1 (MQ=255)
                                     acaAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCAAGGCTAAATACTCctgact                 >  1:7257006/1‑101 (MQ=255)
                                     acaAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCAAGGCTAAATACTCctgact                 >  1:7257005/1‑101 (MQ=255)
                                     acaAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCAAGGCTAAATACTCctgact                 >  1:7257003/1‑101 (MQ=255)
                                      caAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCAAGGCTAAATACTCctgactg                <  2:7256555/101‑1 (MQ=255)
                                       aaaaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCAAGGCTAAATACTCctgactga               <  1:7256613/101‑1 (MQ=255)
                                       aaaaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCAAGGCTAAATACTCctgactga               <  1:7256364/101‑1 (MQ=255)
                                       aaaaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCAAGGCTAAATACTCctgactga               >  2:7257008/1‑101 (MQ=255)
                                        aaaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCAAGGCTAAATACTCCtgactgac              >  1:7257009/1‑101 (MQ=255)
                                         aaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCAAGGCTAAATACTCCTGACTGAcc             >  1:7257011/1‑101 (MQ=255)
                                          aaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCg            >  1:7257012/1‑101 (MQ=255)
                                          aaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCg            >  2:7257013/1‑101 (MQ=255)
                                           aTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGa           <  2:7256484/101‑1 (MQ=255)
                                            tGCACATATTGTGAGCTCGATGAGTGGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGAt          <  2:7256694/101‑1 (MQ=255)
                                            tGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGAt          >  1:7257016/1‑101 (MQ=255)
                                             gCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTTTCTGAATATGGGG‑GGGACCATCCTCCAAGGCTAAATACTCCTGACTGAcc             <  2:7256661/97‑1 (MQ=255)
                                             gCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATa         >  2:7257017/1‑101 (MQ=255)
                                             gCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATa         <  1:7256222/101‑1 (MQ=255)
                                              cacaTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAg        >  1:7257021/1‑101 (MQ=255)
                                              cacaTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAg        >  2:7257019/1‑101 (MQ=255)
                                               acaTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGt       <  2:7256341/101‑1 (MQ=255)
                                               acaTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGt       <  1:7256500/101‑1 (MQ=255)
                                               acaTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGt       >  2:7257022/1‑101 (MQ=255)
                                                caTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTg      >  1:7257023/1‑101 (MQ=255)
                                                caTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTg      <  1:7256345/101‑1 (MQ=255)
                                                caTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTg      <  2:7256593/101‑1 (MQ=255)
                                                 atatTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGa     >  1:7257028/1‑101 (MQ=16)
                                                 atatTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGa     >  1:7257025/1‑101 (MQ=16)
                                                 atatTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGa     <  1:7256512/101‑1 (MQ=16)
                                                 atatTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGa     <  2:7256507/101‑1 (MQ=16)
                                                 atatTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑gg                                                    >  2:7257029/1‑54 (MQ=38)
                                                   atTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAAc   >  1:7257030/1‑101 (MQ=9)
                                                   atTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAAc   <  2:7256782/101‑1 (MQ=9)
                                                    ttGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAAcc  <  1:7256384/101‑1 (MQ=2)
                                                    ttGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGG‑GGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAAcc  <  2:7256429/101‑1 (MQ=2)
                                                                                                     |                                                    
CTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGG‑GGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACC  >  NC_000913/3944008‑3944159

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