Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,109,019 A→T 7.5% I104N (ATT→AAT)  PP_1881 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,109,0190AT 7.5% 242.7 / 8.2 80I104N (ATT→AAT) PP_1881hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (37/37);  new base T (3/3);  total (40/40)
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.30e-01

GAAACATATTGGTCGCCCTTGTGGTCGACCAGTTCGTAGCGAATCCACAGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGA‑GGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGACGGTCAGGTCGGTTTTCAGGCGGTCTTTGTGCA  >  NC_002947/2108875‑2109160
                                                                                                                                                 |                                                                                                                                             
gAAACATATTGGTCGCCCTTGTGGTCGACCAGTTCGTAGCGAATCCACAGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGA‑GGGGCGGGAGAAACTTTTCGACAGGGCCATTGGc                                                                                                                                           >  3:265072/1‑149 (MQ=255)
    aatatTGGTCGCACTTGTGGTCGACCAGTTCGTAGCGAATCCACAGTTTGCCATAGAGTGGAATGGACGCCTGGGGAAAAGACAGTGCTAACGGAAAGTCTACTGGCAAGA‑GCTGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCt                                                                                                                                        >  4:119528/2‑148 (MQ=255)
      tatTGGTCGCCCTTGTGGTCGACCAGTTCGTAGCGAATCCACAGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGA‑GGGGCGGGAGAAACTTTTCGACAGGGCCATTGGCCCTgtc                                                                                                                                     <  4:265072/149‑3 (MQ=255)
          ggTCGCCCTTGTGGTCGACCAGTTCGTAGCGAATCCACAGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGA‑GGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGgt                                                                                                                                  <  7:166550/148‑1 (MQ=255)
            tCGCCCTTGTGGTCGACCAGTTCGTAGCGAATCCACAGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGA‑GGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGtt                                                                                                                               <  3:119528/149‑1 (MQ=255)
                ccTTGTGGTCGACCAGTTCGTAGCGAATCCACAGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGA‑GGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGa                                                                                                                                      <  4:197399/138‑1 (MQ=255)
                ccTTGTGGTCGACCAGTTCGTAGCGAATCCACAGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGA‑GGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGa                                                                                                                                      >  3:197399/1‑138 (MQ=255)
                ccTTGTGGTCGACCAGTTCGTAGCGAATCCACAGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGA‑GGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCg                                                                                                                           <  5:59504/149‑1 (MQ=255)
                ccTTGTGGTCGACCAGTTCGTAGCGAATCCACAGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGA‑GGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCg                                                                                                                           >  3:275147/1‑149 (MQ=255)
                       gTCGACCAGTTCGTAGCGAATCCACAGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGA‑GGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCa                                                                                                                     <  6:27854/148‑1 (MQ=255)
                            ccAGTTCGTAGCGAATCCACAGTTTGCCATTGCGTGGAATGGACGCCAGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGA‑GGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAg                                                                                                               >  4:93717/1‑149 (MQ=255)
                             cAGTTCGTAGCGAATCCACAGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGA‑GGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAgg                                                                                                              >  7:2010/1‑149 (MQ=255)
                              aGTTCGTAGCGAATCCACAGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGA‑GGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGt                                                                                                             >  8:221770/1‑149 (MQ=255)
                              aGTTCGTAGCGAAGCCACAGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGA‑GGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGt                                                                                                             <  8:164054/149‑1 (MQ=255)
                                      gtgAATCCACAGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGA‑GGGGCGGGAGAAACTTTTCGACAGGGCCAATGGACCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCa                                                                                                     >  3:47804/3‑149 (MQ=255)
                                       cGAATCCACAGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGA‑GGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAg                                                                                                    <  8:2010/149‑1 (MQ=255)
                                        gAATCCACAGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGA‑GGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCa                                                                                                     >  3:258352/1‑147 (MQ=255)
                                        gAATCCACAGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGA‑GGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCa                                                                                                     <  4:258352/147‑1 (MQ=255)
                                        gAATCCACAGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGA‑GGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGt                                                                                                   >  7:233160/1‑149 (MQ=255)
                                         aaTCCACAGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGA‑GGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGtt                                                                                                  >  8:300817/1‑149 (MQ=255)
                                          aTCCTCAGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGA‑GGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGtt                                                                                                  <  2:36501/148‑1 (MQ=255)
                                              acaGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGA‑GGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGgt                                                                                                                                  >  2:144715/1‑112 (MQ=255)
                                              acaGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGA‑GGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGgt                                                                                                                                  <  1:144715/112‑1 (MQ=255)
                                              acaGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGA‑GGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAAc                                                                                             >  3:134349/1‑149 (MQ=255)
                                                 gTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGA‑GGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTc                                                                                          <  4:47804/149‑1 (MQ=255)
                                                     gCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGA‑GGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGag                                                                                       >  3:295939/1‑148 (MQ=255)
                                                      ccATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGGCGGGAGAAAATTTTCCACAGGCCCAATCGCCCTCAAGCTGTTTACGCCCCCAGTCAAGGTGACCTGCAGTTGTAACGTCGagt                                                                                      >  1:292043/1‑149 (MQ=255)
                                                           gCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGA‑GGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCa                                                                                                     <  3:269506/128‑1 (MQ=255)
                                                           gCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGA‑GGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCa                                                                                                     >  4:269506/1‑128 (MQ=255)
                                                             gTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGA‑GGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCAttt                                                                              <  8:233160/149‑1 (MQ=255)
                                                                        gCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGA‑GGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGcc                                                                   >  1:20024/1‑149 (MQ=255)
                                                                           tGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGA‑GGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGtt                                                                                                  <  6:158209/115‑1 (MQ=255)
                                                                           tGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGA‑GGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGtt                                                                                                  >  5:158209/1‑115 (MQ=255)
                                                                               gAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGA‑GGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGAt                                                                          <  5:59135/135‑1 (MQ=255)
                                                                               gAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGA‑GGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGAt                                                                          >  6:59135/1‑135 (MQ=255)
                                                                                aaaaGCCAGTGCTAACGGAAAGTCTCCTGGCACGT‑GGGGGGGGAGAAACTTTTCGACAGGGCCAATGGGCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGGAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGct                                                           <  2:292043/149‑2 (MQ=255)
                                                                                   aGCCAGTGCTAACGGAAAGTCTCCTGGCACGA‑GGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTg                                                                             <  6:127520/128‑1 (MQ=255)
                                                                                   aGCCAGTGCTAACGGAAAGTCTCCTGGCACGA‑GGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTg                                                                             >  5:127520/1‑128 (MQ=255)
                                                                                      cAGTGCTAACGGAAAGTCTCCTGGCACGA‑GGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAg                                                     <  2:20024/149‑1 (MQ=255)
                                                                                       aGTGCTAACGGAAAGTCTCCTGGCACGA‑GGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGa                                                    <  4:134349/149‑1 (MQ=255)
                                                                                       aGTGCTAACGGAAAGTCTCCTGGCACGA‑GGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGa                                                    <  7:300817/149‑1 (MQ=255)
                                                                                         tGCTAACGGAAAGTCTCCTGGCACGA‑GGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGAcc                                                  <  2:39238/149‑1 (MQ=255)
                                                                                            tAACGGAAAGTCTCCTGGCACGA‑GGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGGAGTTGTAACGTCGAGTAGTGATTTGGATCACTGCCCGGTGGCAGGTCAgcgcgtg                                               >  2:328925/1‑143 (MQ=255)
                                                                                            tAACGGAAAGTCTCCTGGCACGA‑GGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATg                                               >  1:335832/1‑149 (MQ=255)
                                                                                            tAACGGAAAGTCTCCTGGCACGA‑GGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATg                                               >  6:234275/1‑149 (MQ=255)
                                                                                                ggAAAGTCTCCTGGCACGA‑GGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGag                                                                                       >  3:229143/1‑105 (MQ=255)
                                                                                                ggAAAGTCTCCTGGCACGA‑GGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGag                                                                                       <  4:229143/105‑1 (MQ=255)
                                                                                                ggAAAGTCTCCTGGCACGA‑GGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCggg                                           >  6:38437/1‑149 (MQ=255)
                                                                                                ggAAAGTCTCCTGGCACGA‑GGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTCGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCgg                                            >  6:332961/1‑148 (MQ=255)
                                                                                                ggAAAGTCTCCTGGCACGA‑GGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTCGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCgg                                            <  5:332961/148‑1 (MQ=255)
                                                                                                 gAAAGTCTCCTGGCACGA‑GGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCgggg                                          <  5:38437/149‑1 (MQ=255)
                                                                                                        tcCTGGCACGA‑GGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGa                                   >  3:293710/1‑149 (MQ=255)
                                                                                                          cTGGCACGA‑GGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGACg                                 <  3:24779/149‑1 (MQ=255)
                                                                                                          cTGGCACGA‑GGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGACg                                 >  6:365629/1‑149 (MQ=255)
                                                                                                          cTGGCACGA‑GGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGACg                                 >  4:24779/1‑149 (MQ=255)
                                                                                                          cTGGCACGA‑GGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGACg                                 <  4:275147/149‑1 (MQ=255)
                                                                                                           tGGCACGA‑GGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGACgg                                >  6:359780/1‑149 (MQ=255)
                                                                                                                       gCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGc                                                            >  6:180086/1‑110 (MQ=255)
                                                                                                                       gCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGc                                                            <  5:180086/110‑1 (MQ=255)
                                                                                                                         gggAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGACGGTCAGGTCGGtttt                   <  5:359780/149‑1 (MQ=255)
                                                                                                                         gggAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGACGGTCAGGTCGGtttt                   >  5:328435/1‑149 (MQ=255)
                                                                                                                          ggAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGa                                   <  4:330033/132‑1 (MQ=255)
                                                                                                                          ggAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGa                                   >  3:330033/1‑132 (MQ=255)
                                                                                                                           gagaAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGACGGTCAGGTCGGtttt                   >  6:208418/1‑147 (MQ=255)
                                                                                                                           gagaAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGACGGTCAGGTCGGTTTTCa                 <  5:365629/149‑1 (MQ=255)
                                                                                                                              aaaCTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGACGGTCAGGTCGGTTTTCAGGc              >  8:30654/1‑149 (MQ=255)
                                                                                                                               aaCTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGACGGTCAGGTCGGTTTTCAGGCg             <  5:234275/149‑1 (MQ=255)
                                                                                                                               aaCTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGACGGTCAGGTCGGTTTTCAGGCg             <  1:328925/149‑1 (MQ=255)
                                                                                                                                  ttttCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGACGGTCAGGt                          <  5:243344/133‑1 (MQ=255)
                                                                                                                                  ttttCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGACGGTCAGGt                          >  6:243344/1‑133 (MQ=255)
                                                                                                                                  ttttCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGACGGTCAGGTCGGTTTTCAGGCGGTc          >  6:75833/1‑149 (MQ=255)
                                                                                                                                    ttCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTATCGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGACGGTCAGGTCGGTTTTCAGGCGGTCt         <  5:92644/148‑1 (MQ=255)
                                                                                                                                      cttcAGGGCCATTGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCggg                                           >  3:303629/4‑112 (MQ=255)
                                                                                                                                      cttcAGGGCCATTGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCggg                                           <  4:303629/109‑1 (MQ=255)
                                                                                                                                      cttcAGGGCCATTGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGACGGTCAGGTCGGTTTTCa                 >  3:165709/4‑138 (MQ=255)
                                                                                                                                      cttcAGGGCCATTGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGACGGTCAGGTCGGTTTTCa                 <  4:165709/135‑1 (MQ=255)
                                                                                                                                       gACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGACGGTCAGGTCGGTTTTCAGGCGGTCTTtgt     >  6:243938/1‑149 (MQ=255)
                                                                                                                                         cAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGTACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGTGATCTCGACGGTGAGGTCGGTTTTCAGGCGGTCTTTGTGc   <  3:93717/149‑1 (MQ=255)
                                                                                                                                          aGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGACGGTCAGGTCGGTTTTCAGGCGGTCTTTGTGCa  <  7:30654/149‑1 (MQ=255)
                                                                                                                                          aGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGCTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGACGGTCAGGTCGGTTTTCAGGCGGTCTTTGTGCa  <  7:221770/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                             
GAAACATATTGGTCGCCCTTGTGGTCGACCAGTTCGTAGCGAATCCACAGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGA‑GGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGACGGTCAGGTCGGTTTTCAGGCGGTCTTTGTGCA  >  NC_002947/2108875‑2109160

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: ‑