Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,338,629 T→C 7.4% noncoding (353/395 nt) PP_mr53 → CrcZ ncRNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,338,6290TC 7.4% 192.6 / 4.8 67noncoding (353/395 nt)PP_mr53CrcZ ncRNA
Reads supporting (aligned to +/- strand):  ref base T (34/28);  new base C (3/2);  total (37/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01

ACATCAGCGGCCAAAGCAATCCGTTTGCTCTTTACCCCCGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCTATTCTTGGCCGAACATTCCTTCACACAGTGCATCCCATG  >  NC_002947/5338483‑5338775
                                                                                                                                                  |                                                                                                                                                  
aCATCAGCGGCCAAAGCAATCCGTTTGCTCTTTACCCCCGTTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCtt                                                                                                                                                   <  8:125134/148‑1 (MQ=255)
  aTCAGCGGCCAAAGCAATCCGTTTGCTCTTTACCCCCGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAgg                                                                                                                                                <  7:148602/149‑1 (MQ=255)
     aGCGGCCAAAGCAATCCGTTTGCTCTTTACCCCCGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTc                                                                                                                                             <  1:111295/149‑1 (MQ=255)
      gCGGCCAAAGCAATCCGTTTGCTCTTTACCCCCGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTcc                                                                                                                                            >  7:64642/1‑149 (MQ=255)
         gCCAAAGCAATCCGTTTGCTCTTTACCCCCGGTTTGGGGTTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCCTAAGCTcccca                                                                                                                                         <  3:272046/149‑2 (MQ=255)
          ccAAAGCAATCCGTTTGCTCTTTACCCCCGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCCTAAGCTccccaa                                                                                                                                        <  3:48833/149‑3 (MQ=255)
            aaaGCAATCCGTTTGCTCTTTACCCCCGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGt                                                                                                                                      <  8:64642/149‑1 (MQ=255)
            aaaGCAATCCGTTTGCTCTTTACCCCCGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGt                                                                                                                                      >  2:192610/1‑149 (MQ=255)
             aaGCAATCCGTTTGCTCTTTACCCCCGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTc                                                                                                                                     <  3:124898/149‑1 (MQ=255)
                 aaTCCGTTTGCTCTTTACCCCCGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGttt                                                                                                                                 <  5:169075/149‑1 (MQ=255)
                      gTTTGCTCTTTACCCCCGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACCATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTAcc                                                                                                                            >  8:144341/1‑149 (MQ=255)
                       tttGCTCTTTACCCCCGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTAccc                                                                                                                           <  6:369712/149‑1 (MQ=255)
                                 aCCCCCGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTggcag                                                                                                                  <  6:70260/148‑1 (MQ=255)
                                   ccccGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTggcaggca                                                                                                               >  2:204980/1‑149 (MQ=255)
                                   ccccGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTggcaggca                                                                                                               >  1:8669/1‑149 (MQ=255)
                                    cccGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTc                                                                                                                                <  4:63855/131‑1 (MQ=255)
                                    cccGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTc                                                                                                                                >  3:63855/1‑131 (MQ=255)
                                     ccGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATc                                                                                                             >  8:292117/1‑149 (MQ=255)
                                            gggggTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTa                                                                                                       <  1:192610/148‑1 (MQ=255)
                                                    tcCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCATa                                                                                               >  1:346092/1‑148 (MQ=255)
                                                      cacaAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTaca                                                                                             <  8:131020/148‑1 (MQ=255)
                                                      cacaAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGAACGCAACGACTTCTTGGGGAGATTAGGCTCCACTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTacac                                                                                            >  1:86976/1‑149 (MQ=255)
                                                        caAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTa                                                                                                                                                  >  3:200130/1‑93 (MQ=255)
                                                        caAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTa                                                                                                                                                  <  4:200130/93‑1 (MQ=255)
                                                            cTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTcc                                                                                      >  6:343439/1‑149 (MQ=255)
                                                             ttCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTa                                                                                               <  3:56725/139‑1 (MQ=255)
                                                             ttCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTa                                                                                               >  4:56725/1‑139 (MQ=255)
                                                              tCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTAc                                                                                                      >  8:11952/1‑131 (MQ=255)
                                                               cGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTAcc                                                                                                                            <  3:32606/108‑1 (MQ=255)
                                                               cGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTAcc                                                                                                                            >  4:32606/1‑108 (MQ=255)
                                                                gagaTTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTa                                                                                               <  4:184877/136‑1 (MQ=255)
                                                                gagaTTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTa                                                                                               >  3:184877/1‑136 (MQ=255)
                                                                 agaTTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAg                                                                                 >  4:292032/1‑149 (MQ=255)
                                                                      tgtgGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTaaa                                                                             <  5:343439/148‑1 (MQ=255)
                                                                        tgGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGc                                                                                                  >  7:74138/1‑125 (MQ=255)
                                                                        tgGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGc                                                                                                  <  8:74138/125‑1 (MQ=255)
                                                                          gACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGc                                                                          >  8:242666/1‑147 (MQ=255)
                                                                               ggCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATc                                                                   >  4:367001/1‑149 (MQ=255)
                                                                               ggCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATc                                                                   >  6:239363/1‑149 (MQ=255)
                                                                                   cctatctaaaaCAAGAAGCCCAGCAAAAAACCAATACTAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCATTGTCGTTTCTACCCTTCTGGCATGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATcccc                                                                <  6:226833/141‑1 (MQ=255)
                                                                                          aaaaCAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGCGCTTAGGGTCCCCTTGTCGTTTCTACGGTTCTGGCAGGGATGCTGCTAGCCGCCTAGACCATTCCCCGAGTTAATGTTAGAATCGCCGCCcatga                                                        >  4:290528/1‑147 (MQ=255)
                                                                                          aaaaCAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTcccc                                                                                    <  4:143551/121‑1 (MQ=255)
                                                                                          aaaaCAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTcccc                                                                                    >  3:143551/1‑121 (MQ=255)
                                                                                           aaaCAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACACGCCTACACCATTCCCCGAGAAAATGCTAGAATcccc                                                                >  5:226833/1‑140 (MQ=255)
                                                                                                 gaagCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGc                                                  >  7:210767/1‑148 (MQ=255)
                                                                                                      ccAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATtct                                            >  4:78998/1‑149 (MQ=255)
                                                                                                         gCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCt                                         <  7:292117/149‑1 (MQ=255)
                                                                                                             aaaaaCAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTa                                                                        <  6:2646/114‑1 (MQ=255)
                                                                                                             aaaaaCAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTa                                                                        >  5:2646/1‑114 (MQ=255)
                                                                                                               aaaCAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGc                                                              <  3:68285/122‑1 (MQ=255)
                                                                                                               aaaCAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGc                                                              >  4:68285/1‑122 (MQ=255)
                                                                                                                        gagCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCTATTCTTGGCCGAACa                          >  7:183965/1‑149 (MQ=255)
                                                                                                                            aCGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCTATTCTTGGCCGAACATTcc                      >  4:96747/1‑149 (MQ=255)
                                                                                                                            aCGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCTATTCTTGGCCGAACATTcc                      >  4:102128/1‑149 (MQ=255)
                                                                                                                               cAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCTATTCTTGGCCGAACATTCCTTc                   <  1:204980/149‑1 (MQ=255)
                                                                                                                                aaCGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAg                                                                                 >  6:354469/1‑86 (MQ=255)
                                                                                                                                aaCGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAg                                                                                 <  5:354469/86‑1 (MQ=255)
                                                                                                                                 aCGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCTATTCTTGGCCGAACATTCCTTcac                 <  7:144341/149‑1 (MQ=255)
                                                                                                                                    aCTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCTATTCTTGGCCGAACATTCCTTcacaca              <  3:78998/149‑1 (MQ=255)
                                                                                                                                      ttcttGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCTATTCTTGGCCGAACATTCCTTCACACAGt            >  6:284575/1‑149 (MQ=255)
                                                                                                                                      ttcttGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCTATTCTTGGCCGAACATTCCTTCACACAGt            >  6:134890/1‑149 (MQ=255)
                                                                                                                                      gacaaggggAGCCTAAGCTCCCCTTGTCGTTTCTGCCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCTATTCTTGGCCGAACATTCCTTCACACAGt            >  6:374221/6‑149 (MQ=255)
                                                                                                                                      gacaaggggAGCCTAAGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCTATTCTTGGCCGAACATTCCTTCACACAGt            >  4:155448/6‑149 (MQ=255)
                                                                                                                                        caaggggAGCCTAAGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCTATTCTTGGCCGAACATTCCTTCACACAGTg           >  2:137111/4‑148 (MQ=255)
                                                                                                                                            gggAGCTTAGGCGCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCTATTCTTGGCCGAACATTCCTTCACACAGTGCATc       <  3:274585/148‑1 (MQ=255)
                                                                                                                                                gCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCTATTCTTGGCCGAACATTCCTTCACACAGTGCATCCCATg  <  3:292032/149‑1 (MQ=255)
                                                                                                                                                gCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCTATTCTTGGCCGAACATTCCTTCACACAGTGCATCCCATg  >  1:305275/1‑149 (MQ=255)
                                                                                                                                                gCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCTATTCTTGGCCGAACATTCCTTCACACAGTGCATCCCATg  <  1:151707/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
ACATCAGCGGCCAAAGCAATCCGTTTGCTCTTTACCCCCGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCTATTCTTGGCCGAACATTCCTTCACACAGTGCATCCCATG  >  NC_002947/5338483‑5338775

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