Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,213,182 A→T 5.0% intergenic (+23/‑52) rrfE → / → yjaA 5S ribosomal RNA of rrnE operon/stress‑induced protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,213,1820AT 5.0% 258.3 / 2.2 90intergenic (+23/‑52)rrfE/yjaA5S ribosomal RNA of rrnE operon/stress‑induced protein
Reads supporting (aligned to +/- strand):  ref base A (42/43);  new base T (3/2);  total (45/45)
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

AGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGTAA  >  NC_000913/4213034‑4213328
                                                                                                                                                    |                                                                                                                                                   
aGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCAt                                                                                                                                                    >  2:949986/1‑150 (MQ=255)
 gAATTTGCCTGGCGGGAGTAGCGCGGTGGTCCCTCCTGACCCCATGACGAACTCAGAAGTGAAACGCCGTATCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATa                                                                                                                                                   <  1:963275/150‑1 (MQ=255)
     ttGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑c                                                                                                                                              <  1:1082046/150‑1 (MQ=255)
       gCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCg                                                                                                                                            >  2:1057945/1‑150 (MQ=255)
        ccTGGCGGCTGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGTCATCAAAATAGAAAAACCCCGGTCCATAAGG‑CCgg                                                                                                                                           <  2:907592/150‑1 (MQ=255)
        ccTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCgg                                                                                                                                           <  1:445925/150‑1 (MQ=255)
         cTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCggg                                                                                                                                          >  2:573341/1‑150 (MQ=255)
         cTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGGCCATAAGG‑CCggg                                                                                                                                          >  1:1124529/1‑150 (MQ=255)
          tGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCgggg                                                                                                                                         <  2:273994/150‑1 (MQ=255)
          tGGCGCCCGTAGCGCGGGGGGACACCCTGACCCCATGACTCACTCAGAACTGAAACGCCGTCGCGCCGACGGGAGTGTGGGGGCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGGCCATAAGG‑CCgggg                                                                                                                                         <  2:67618/150‑1 (MQ=14)
            gcggcAGTAGCGCGGGGGGCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGGAGTGTGGGGTCTCCCCATGCGATAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGG‑CCTTATGGACCGGGGtt                                                                                                                                       <  1:638920/150‑1 (MQ=14)
              ggcAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGtttt                                                                                                                                     >  1:1005933/1‑150 (MQ=255)
              ggcAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGtttt                                                                                                                                     >  2:259752/1‑150 (MQ=255)
                 aGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTg                                                                                                                                  <  1:1000145/150‑1 (MQ=255)
                  gTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGc                                                                                                                                 <  1:100455/150‑1 (MQ=255)
                  gTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCAGAAGG‑CCGGGGTTTTTTgg                                                                                                                                 >  1:903914/1‑149 (MQ=255)
                  gTAGCGCGGTGGTCCCACATGACCCCTTGCCGACATCAGAAGTGAAACGCCCTAGCGCCGATGGGAGTGGGGGGTCTCCCCATGCGAGAGTAGGTAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGc                                                                                                                                 <  2:221616/150‑1 (MQ=18)
                      cgcgGTGGTCCCACCTGACCCCATGTCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCatat                                                                                                                             >  1:63720/1‑150 (MQ=255)
                      cgcgGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCatat                                                                                                                             >  1:688950/1‑150 (MQ=255)
                            ggtCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAatta                                                                                                                       >  2:923476/1‑150 (MQ=255)
                            ggtCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCAGTAGCGCCGTTGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAatta                                                                                                                       <  2:509916/150‑1 (MQ=255)
                               cccACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTAttt                                                                                                                    >  1:1162927/1‑150 (MQ=255)
                               cccACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTAttt                                                                                                                    <  2:120420/150‑1 (MQ=255)
                                  aCCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCa                                                                                                                 <  1:1186256/150‑1 (MQ=255)
                                       aCCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCatgatg                                                                                                            <  1:1218431/150‑1 (MQ=255)
                                        ccccATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCAtgatga                                                                                                           >  1:1123228/1‑150 (MQ=255)
                                        ccccATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCAtgatga                                                                                                           <  1:392229/150‑1 (MQ=255)
                                         cccATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGaa                                                                                                          <  1:352424/150‑1 (MQ=255)
                                           cATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGTTGGTAGTGTGTGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAgg                                                                                                        <  2:1197313/150‑1 (MQ=255)
                                            aTGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAggg                                                                                                       >  2:1179878/1‑150 (MQ=255)
                                            aTGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGGCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAggg                                                                                                       <  2:303324/150‑1 (MQ=255)
                                                    cTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCAt                                                                                               >  1:156459/1‑150 (MQ=255)
                                                               aaCGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATGATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGt                                                                                    >  2:648900/1‑150 (MQ=255)
                                                                 cGCCTTAGCGCCGATTGTAGTGTGGGGTCTCCCCATGCGAGAGTATGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGAACATAAGT‑CAGTGGTTTGTTGCATATCAATTATTTGCATTATGAAGGGAATCTCATGTCAGTTCTGtat                                                                                  >  2:1238065/1‑150 (MQ=25)
                                                                   ccGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTatat                                                                                <  2:1036577/150‑1 (MQ=255)
                                                                    cGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCGATAAGG‑CCGGGGTTTTTTGCATATCATTTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATc                                                                               >  1:690451/1‑150 (MQ=255)
                                                                    cGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATc                                                                               <  2:51444/150‑1 (MQ=255)
                                                                      tAGCTCCGTTGGTTGTGGTGTGTCTCCCGATGCGAGAGTAGGGAACTGCCTTGAATGATATTAGAATAAGCGCGGTCCATAATG‑CCGGGGTTTTTTGCATTTCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCa                                                                             <  1:248411/150‑1 (MQ=255)
                                                                      tAGCGCAGTTGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCCGGCATCAATTAAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCa                                                                             <  2:664758/150‑1 (MQ=38)
                                                                         cgcCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAAt                                                                          >  2:1044386/1‑150 (MQ=255)
                                                                          gcCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAAtt                                                                         >  2:1100008/1‑150 (MQ=255)
                                                                            cGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCATTTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATtcg                                                                       <  1:1067616/150‑1 (MQ=255)
                                                                            cGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATtcg                                                                       <  2:1005933/150‑1 (MQ=255)
                                                                                ggTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTcgtcgt                                                                   >  1:990511/1‑150 (MQ=255)
                                                                                    gtgtGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTGCATAAGG‑CCGGGGTTTTTTGCATAACAATTATTTGCATGATGAAGGGAATCTCCTGTCAGTTCTGTAGATCCAAATTCGTCGTAAta                                                               >  2:112081/1‑149 (MQ=255)
                                                                                    gtgtGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATc                                                               >  2:1100292/1‑150 (MQ=255)
                                                                                      gtgGGCTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGCGGTTTTTTGCATATCAATTATTTGCATGATGAAGCGAATCTCATCTCAGTTCTGTATATCCAAATTCGTCGTAATCaa                                                             >  2:691853/1‑150 (MQ=255)
                                                                                       tgGGGGCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCaaa                                                            <  2:1015202/150‑1 (MQ=255)
                                                                                         gggTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTATAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAAtt                                                          <  1:253380/150‑1 (MQ=255)
                                                                                              tcCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGt                                                     >  1:485283/1‑150 (MQ=255)
                                                                                               ccccATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGtt                                                    <  2:158775/150‑1 (MQ=255)
                                                                                                cccATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTCATCAAATTACTGTTc                                                   >  1:1317663/1‑150 (MQ=255)
                                                                                                     gCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGAt                                              >  1:1165412/1‑150 (MQ=255)
                                                                                                     gCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGAt                                              <  2:310186/150‑1 (MQ=255)
                                                                                                     gCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGAt                                              <  2:774084/150‑1 (MQ=255)
                                                                                                      cGAGAGTATGGTACTTCCAGGCATCAAATTAGAAAACCCCCGGTCCATCAGG‑CCGGGGTTTTTTTCATTTAAAATATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATCTCCAAAATCGGCGTAATCAAATTACTGTTCGCGATc                                             <  2:110754/150‑1 (MQ=255)
                                                                                                          agTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGa                                         <  1:844123/150‑1 (MQ=255)
                                                                                                             cgggAACTGCAAGCCATCAAATTAGAAAAACCCCGGGCCTTAAGG‑CAGGGGTTTTTTGCATATCAATTATTTGCATGCTGAAGTGAATCTACTGTCAGTTCTGTATATCCAAATTCGGCGTAATCAAATTCCTGTTCGCGATCTTGAAAg                                      <  2:176866/149‑1 (MQ=255)
                                                                                                              gggAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCTAATTACTGTTCGCGATCTTGAAAGc                                     <  1:729336/150‑1 (MQ=255)
                                                                                                              gggAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGc                                     >  1:767744/1‑150 (MQ=255)
                                                                                                              gggAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCAT‑AGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAACGGAATCTCATGTCAGTTCAGTATATCCAAATTCGTAGTAATCAAATTACTGTTCGCGATCTTGAAAGCa                                    >  2:1252882/1‑150 (MQ=255)
                                                                                                               ggAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCa                                    <  2:688950/150‑1 (MQ=255)
                                                                                                                   cTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACg                                <  1:1013019/150‑1 (MQ=255)
                                                                                                                   cTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACg                                <  2:1162927/150‑1 (MQ=255)
                                                                                                                     gCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTg                              >  1:575740/1‑150 (MQ=255)
                                                                                                                      ccAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTAGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGa                             >  2:690728/1‑150 (MQ=255)
                                                                                                                       cAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAACCGTGaa                            >  1:261685/1‑150 (MQ=255)
                                                                                                                       cAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGaa                            >  2:377283/1‑150 (MQ=255)
                                                                                                                       cAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGaa                            <  1:259752/150‑1 (MQ=255)
                                                                                                                        aGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGTAAACGTGAAg                           <  1:366661/150‑1 (MQ=255)
                                                                                                                        aGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTAGTCGTAATCAATTTACTGTTCGCGATCTTGAAAGCAAACGTGAAg                           <  2:623815/150‑1 (MQ=255)
                                                                                                                          gCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAgtg                         <  1:923412/150‑1 (MQ=255)
                                                                                                                             tCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCa                      <  1:1035772/150‑1 (MQ=255)
                                                                                                                              cAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAg                     >  1:1050561/1‑150 (MQ=255)
                                                                                                                               aaaTTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAgg                    >  1:586628/1‑150 (MQ=255)
                                                                                                                                 aTTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGgag                  >  2:1113252/1‑150 (MQ=255)
                                                                                                                                   gcaaaaaaCCCCGG‑CCTTATGGACCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGAc                >  1:270431/4‑150 (MQ=255)
                                                                                                                                     gAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGc              <  2:39770/150‑1 (MQ=255)
                                                                                                                                     aaaaaaCCCCGG‑CCTTATGGACCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGCCGc              >  1:335199/2‑150 (MQ=255)
                                                                                                                                     aaaaaaCCCCGG‑CCTTATGGACCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGc              >  1:697864/2‑150 (MQ=255)
                                                                                                                                       aaaaCCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTg            >  2:1258108/1‑150 (MQ=255)
                                                                                                                                       aaaaCCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAAAGTGAAGTGTCAGGAGACGCTg            >  2:1200954/1‑150 (MQ=255)
                                                                                                                                         aaCCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTATTCAAATTACTGTTCGCGATCTTGAACGCAAACGTGAAGTGTCAGGAGACGCTGcc          >  1:816827/1‑150 (MQ=255)
                                                                                                                                         aaCCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGcc          <  1:397933/150‑1 (MQ=255)
                                                                                                                                           ccccGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGATTCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCtt        >  2:1254700/1‑150 (MQ=255)
                                                                                                                                            cccGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTc       >  2:812311/1‑150 (MQ=255)
                                                                                                                                            cccGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTc       >  2:517909/1‑150 (MQ=255)
                                                                                                                                              cGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAg     <  1:492590/150‑1 (MQ=255)
                                                                                                                                               ggTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGCCGCTGCCTTCAGt    >  2:566917/1‑150 (MQ=255)
                                                                                                                                                 tCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGTaa  <  2:468033/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGTAA  >  NC_000913/4213034‑4213328

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑