Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,905,846 C→T 6.4% M59I (ATG→ATA)  PP_4314 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,905,8460CT 6.4% 183.9 / 2.0 62M59I (ATG→ATA) PP_4314hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (28/30);  new base T (2/2);  total (30/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

GCCGGTAGTACAGGTGCGTACCACAGCGGCCACAGAAAGCACGCTGGGCCCAGTCGGACGAATCATAGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCACAGGCAGGTGCCGTGTTTTTCCAGG‑GCCATGGGTGGCTCCTTGAACGGAAAGATGGGCCTGCAATTATAGAAGCCACCCAGGG  >  NC_002947/4905698‑4905994
                                                                                                                                                    |                                                                                                                                                     
gCCGGTAGTACAGGTGCGTACCACAGCGGCCACAGAAAGCACGCTGGGCCCAGTCGGACGAATCATAGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGcac                                                                                                                                                       <  4:46224/149‑1 (MQ=255)
  cGGTAGTACAGGTGCGTACCACAGCGGCCACAGAAAGCACGCTGGGCCCAGTCGGACGAATCATAGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACAt                                                                                                                                                     <  6:636/149‑1 (MQ=255)
   ggTAGTACAGGTGCGTACCACAGCGGCCACAGAAAGCACGCTGGGCCCAGTCGGACGAATCATAGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACAtg                                                                                                                                                    >  1:283436/1‑149 (MQ=255)
     tagtaCAGGTGCGTACCACAGCGGCCACAGAAAGCACGCTGGGCCCAGTCGGACGAATCATAGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACAtgct                                                                                                                                                  <  5:181317/149‑1 (MQ=255)
       gtaCAGGTGCGTACCACAGCGGCCACAGAAAGCACGCTGGGCCCAGTCGGACGAATCATAGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACAtgctgc                                                                                                                                                <  3:111898/149‑1 (MQ=255)
          cAGGTGCGTACCACAGCGGCCACAGAAAGCACGCTGGGCCCAGTCGGACGAATCATAGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCAGCATGTGCCGCaa                                                                                                                                              <  6:325729/148‑1 (MQ=255)
            ggTGCGTACCACAGCGGCCACAGAAAGCACGCTGGGCCCAGTCGGACGAATCATAGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGa                                                                                                                                           >  7:257656/1‑149 (MQ=255)
             gTGCGTACCACAGCGGCCACAGAAAGCACGCTGGGCCCAGTCGGACGAATCATAGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGAc                                                                                                                                          <  8:257656/149‑1 (MQ=255)
              tGCGTACCACAGCGGCCACAGAAAGCACGCTGGGCCCAGTCGGACGAATCATAGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCAGCATGTGCCGCaagt                                                                                                                                            >  7:128850/1‑144 (MQ=255)
              tGCGTACCACAGCGGCCACAGAAAGCACGCTGGGCCCAGTCGGACGAATCATAGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCAGCATGTGCCGCaagt                                                                                                                                            <  8:128850/146‑3 (MQ=255)
                    ccACAGCGGCCACAGAAAGCACGCTGGGCCCAGTCGGACGAATCATAGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCt                                                                                                                                   >  4:122942/1‑149 (MQ=255)
                        aGCGGCCACAGAAAGCACGCTGGGCCCAGTCGGACGAATCATAGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGata                                                                                                                               <  1:112056/149‑1 (MQ=255)
                         gCGGCCACAGAAAGCACGCTGGGCCCAGTCGGACGAATCATAGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGTACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGatat                                                                                                                              <  1:255946/149‑1 (MQ=255)
                           ggCCACAGAAAGCACGCTGGGCCCAGTCGGACGAATCATAGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATgg                                                                                                                            <  3:179288/149‑1 (MQ=255)
                           ggCCACAGAAAGCACGCTGGGCCCAGTCGGACGAATCATAGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATgg                                                                                                                            <  8:63593/149‑1 (MQ=255)
                            gCCACAGAAAGCACGCTGGGCCCAGTCGGACGAATCATAGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATggg                                                                                                                           <  1:196086/149‑1 (MQ=255)
                                caGAAAGCACGCTGGGCCCAGTCGGACGAATCATAGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACAtgctg                                                                                                                                                 <  2:7956/123‑1 (MQ=255)
                                 aGAAAGCACGCTGGGCCCAGTCGGACGAATCATAGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTg                                                                                                                      >  5:118743/1‑149 (MQ=255)
                                       cACGCTGGGCCCAGTCGTATGAATCATAGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACg                                                                                                                <  7:288395/149‑1 (MQ=255)
                                        aCGCTGGGCCCAGTCGGACGAATCATAGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACgg                                                                                                               >  3:267859/1‑149 (MQ=255)
                                               gCCCAGTCGGACGAATCATAGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGt                                                                                                        >  5:132855/1‑149 (MQ=255)
                                                 ccAGTCGGACGAATCATAGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCa                                                                                                      <  1:154150/149‑1 (MQ=255)
                                                       ggACGAATCATAGACGCTGGGTGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCcgcg                                                                                                >  2:119723/1‑149 (MQ=255)
                                                       ggACGAATCATAGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCcgcg                                                                                                >  7:66690/1‑149 (MQ=255)
                                                       ggACGAATCATAGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCcgc                                                                                                 <  3:122942/148‑1 (MQ=255)
                                                       ggACGAATCATAGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCcgc                                                                                                 <  6:291088/148‑1 (MQ=255)
                                                        gACGAATCATAGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCgt                                                                                               <  6:118743/149‑1 (MQ=255)
                                                           gAATCATAGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCgt                                                                                               >  5:174298/1‑146 (MQ=255)
                                                           gAATCATAGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCgt                                                                                               <  6:174298/146‑1 (MQ=255)
                                                           gAATCATAGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTc                                                                                            >  3:119755/1‑149 (MQ=255)
                                                                 tAGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGca                                                                                      >  7:268320/1‑149 (MQ=255)
                                                                  aGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGca                                                                                      >  1:303773/1‑148 (MQ=255)
                                                                   gACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGca                                                                                      >  7:83109/1‑147 (MQ=255)
                                                                      gCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCAcaggc                                                                                 <  4:267859/149‑1 (MQ=255)
                                                                      gCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCAcagg                                                                                  >  4:165254/1‑148 (MQ=255)
                                                                      gCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCAcagg                                                                                  >  1:153684/1‑148 (MQ=255)
                                                                         gggCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCAcaggcagg                                                                              >  4:70349/1‑149 (MQ=255)
                                                                             gcACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCACAGGCAGGTGcc                                                                          >  7:74388/1‑149 (MQ=255)
                                                                                   ccTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCACAGGCAGGTGCCGTGttt                                                                    <  3:78345/149‑1 (MQ=255)
                                                                                    cTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCACAGGCAGGTGCCGTGtttt                                                                   <  3:165254/149‑1 (MQ=255)
                                                                                        gATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATTCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCACAGGCAGGTGCCGTGTTTTTCCa                                                               <  2:249197/149‑1 (MQ=255)
                                                                                          tCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCACAGGCAGGTGCCGTGTTTTTCCAgg                                                             >  5:327720/1‑149 (MQ=255)
                                                                                           cAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCACAGGCAGGTGCCGTGTTTTTCCAgg‑g                                                           >  2:164315/1‑149 (MQ=255)
                                                                                            aGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCACAGGCAGGTGCCGTGTTTTTCCAGG‑Gc                                                          <  1:119723/149‑1 (MQ=255)
                                                                                             gCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCACAGGCAGGTGCCGTGTTTTTCCAGG‑Gc                                                          >  6:234304/1‑148 (MQ=255)
                                                                                                 aGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCACAGGCAGGTGCCGTGTTTTTCCAGG‑GCCATgg                                                     <  8:268320/149‑1 (MQ=255)
                                                                                                   gCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCACAGGCAGGTGCCGTGTTTTTCCAGG‑GCCATGGGt                                                   <  2:153684/149‑1 (MQ=255)
                                                                                                         tgcAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCACAGGCAGGTGCCGTGTTTTTCCAGG‑GCCATGGGTGGCTcc                                             >  2:202508/1‑149 (MQ=255)
                                                                                                         tgcAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCACAGGCAGGTGCCGTGTTTTTCCAGG‑GCCATGGGTGGCTcc                                             >  4:135013/1‑149 (MQ=255)
                                                                                                            aGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCACAGGCAGGTGCCGTGTTTTTCCAGG‑GCCATGGGTGGCTCCtt                                           <  6:327720/148‑1 (MQ=255)
                                                                                                            aGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCACAGGCAGGTGCCGTGTTTTTCCAGG‑GCCATGGGTGGCTCCTTg                                          <  8:175216/149‑1 (MQ=255)
                                                                                                                   caccaGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCACAGGCAGGTGCCGTGTTTTTCCAGG‑GCCATGGGTGGCTCCTTGAACGGa                                    >  3:21233/1‑148 (MQ=255)
                                                                                                                   caccaGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCACAGGCAGGTGCCGTGTTTTTCCAGG‑GCCATGGGTGGCTCCTTGAACGGa                                    >  4:270251/1‑148 (MQ=255)
                                                                                                                   caccaGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCACAGGCAGGTGCCGTGTTTTTCCAGG‑GCCATGGGTGGCTCCTTGAAAGGa                                    >  2:252471/1‑148 (MQ=255)
                                                                                                                      cagcagCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCACAGGCAGGTGCCGTGTTTTTCCAGG‑GCCATGGGTGGCTCCTTGAACGGa                                    >  1:314115/1‑145 (MQ=255)
                                                                                                                        gcagcGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCACAGGCAGGTGCCGTGTTTTTCCAGG‑GCCATGGGTGGCTCCTTGAACGGAAAGATg                              <  8:74388/149‑1 (MQ=255)
                                                                                                                          agcGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCACAGGCAGGTGCCGTGTTTTTCCAGG‑GCCATGGGTGGCTCCTTGAACGGAAAGATggg                            <  4:21233/149‑1 (MQ=255)
                                                                                                                                ccgccgGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCACAGGCAGGTGCCGTGTTTTTCAAGG‑GCCATGGGTGGCTCCTTGAACGGAAAGATGGGCCTGc                       >  4:5221/1‑148 (MQ=255)
                                                                                                                                   ccgGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGAATTGTCGACGGTGACCGTCAGCCGCTTGTCGCCGCACAGGCAGGTGCCTTGTTTTTCCAGC‑GCCATGGGTGGCTCCTTGAACGGAAAGATGGGCCTGCAAtt                   <  2:283436/149‑1 (MQ=255)
                                                                                                                                      ttgtcattGCAGCATGTGCCGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCACAGGCAGGTGCCGTGTTTTTCCAGG‑GCCATGGGTGGCTCCTTGAACGGAAAGATGGGCCTGCAATta                  >  5:265507/7‑147 (MQ=255)
                                                                                                                                      gTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCACAGGCAGGTGCCGTGTTTTTCCAGG‑GCCATGGGTGGCTCCTTGAACGGAAAGATGGGCCTGCAATtata                <  3:270251/149‑1 (MQ=255)
                                                                                                                                        ccACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCACAGGCAGGTGCCGTGTTTTTCCAGG‑GCCATGGGTGGCTCCTTGAAc                                       >  8:158935/1‑124 (MQ=255)
                                                                                                                                           cTTGCGGCACATTCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCACAGGCAGGTGCCGTGTTTTTCCAGG‑GCCATGGGTGGCTCCTTGAACGGAAAGATGGGCCTGCAATTATAGAAGc           <  6:307329/149‑1 (MQ=255)
                                                                                                                                                    caTGCTGGAATGTCAAGCGCTGTTATTGGTTTTGTCGTCGGTGACCGTCAGCCGGGTGTCGCCGCACAGGCAGGTGCCGTGTTTTTCCAGGTATCATGGGTGGCTCCTTGAACGGAAAGATGGGCCTGCAATTATAGAAGCCACCCAgg   <  1:252471/149‑1 (MQ=255)
                                                                                                                                                    caTGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCACAGGCAGGTGCCGTGTTTTTCCAGG‑GCCATGGGTGGCTCCTTGAACGGAAAGATGGGCCTGCAATTATAGAAGCCACCAAggg  <  4:119755/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCCGGTAGTACAGGTGCGTACCACAGCGGCCACAGAAAGCACGCTGGGCCCAGTCGGACGAATCATAGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCACAGGCAGGTGCCGTGTTTTTCCAGG‑GCCATGGGTGGCTCCTTGAACGGAAAGATGGGCCTGCAATTATAGAAGCCACCCAGGG  >  NC_002947/4905698‑4905994

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