Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,757,042 T→A 7.6% I25F (ATC→TTC)  PP_3320 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,757,0420TA 7.6% 176.3 / 3.3 62I25F (ATC→TTC) PP_3320hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (28/29);  new base A (3/2);  total (31/31)
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

GTTGGATGCCAGCACGATATGCTGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGATTTGCAACCTGCGATAGATGTTATGTTAT  >  NC_002947/3756894‑3757190
                                                                                                                                                    |                                                                                                                                                    
gTTGGATGCCAGCACGATATGCTGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGAt                                                                                                                                                      <  5:233113/149‑1 (MQ=255)
    gATGCCAGCACGATATGCTGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCa                                                                                                                                                  <  5:36262/149‑1 (MQ=255)
    gATGCCAGCACGATATGCTGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCACCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTc                                                                                                                                                   <  1:267379/148‑1 (MQ=255)
      tGCCAGCACGATATGCTGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCAtg                                                                                                                                                <  3:292232/149‑1 (MQ=255)
      tGCCAGCACGATATGCTGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCAtg                                                                                                                                                <  7:199076/149‑1 (MQ=255)
             aCGATATGCTGACATTTCCCTCTTCTCCACATCGCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTgg                                                                                                                                         <  4:209005/149‑1 (MQ=255)
              cGATATGCTGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTggg                                                                                                                                        >  1:275856/1‑149 (MQ=255)
                   tGCTGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGt                                                                                                                                       >  7:172276/1‑145 (MQ=255)
                   tGCTGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGt                                                                                                                                       <  8:172276/145‑1 (MQ=255)
                   tGCTGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGAt                                                                                                                                    <  5:144692/148‑1 (MQ=255)
                    gCTGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATgg                                                                                                                                  >  7:322981/1‑149 (MQ=255)
                      tGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCg                                                                                                                                >  7:277551/1‑149 (MQ=255)
                        aCATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTgg                                                                                                                                         <  2:340726/138‑1 (MQ=255)
                        aCATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTgg                                                                                                                                         >  1:340726/1‑138 (MQ=255)
                            ttCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGAt                                                                                                                           >  2:298957/1‑148 (MQ=255)
                                      tcCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGc                                                                                                                                 <  5:365278/132‑1 (MQ=255)
                                      tcCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGc                                                                                                                                 >  6:365278/1‑132 (MQ=255)
                                        cacaTCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGctgctg                                                                                                                                          <  2:324140/121‑1 (MQ=255)
                                        cacaTCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGAGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAggggg                                                                                                              >  8:283591/1‑149 (MQ=255)
                                              ccGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAg                                                                                                        <  3:153569/149‑1 (MQ=255)
                                                        tAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGc                                                                                              >  8:36846/1‑149 (MQ=255)
                                                        tAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGc                                                                                              >  2:164698/1‑149 (MQ=255)
                                                            gTAACAAGCGCCCAGGGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTCGGTGATGGCGGGGATCAGCCCCATGGGGTTACAGGGTATGTGGCGAAg                                                                                          >  4:246428/1‑149 (MQ=255)
                                                            gTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAg                                                                                          <  8:192057/149‑1 (MQ=255)
                                                               aCAAGCGCCCAGAGTTGCTGAGCACCATCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGcac                                                                                       <  8:322981/149‑1 (MQ=255)
                                                               aCAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGGGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGcac                                                                                       <  3:243179/149‑1 (MQ=255)
                                                                 aaGCGCCCAGAGTTGCTGAGCACCAGAAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGaaca                                                                                      >  4:46490/1‑148 (MQ=255)
                                                                   gcgcCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGACGATGTCATGCTGCTGGGTGATGGCGGGGAGCAGCCCCCGGGGGTTCCAGGGTATGTGGAGAAGCACAGCg                                                                                   >  2:100596/1‑149 (MQ=255)
                                                                        cAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTc                                                                                >  6:142075/1‑147 (MQ=255)
                                                                                        cAGCAATGTACTGAAGCAGCGCATTTTCGTTCATCCAGGCAAGAATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCACAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTccc                                                              >  6:88813/1‑149 (MQ=255)
                                                                                          gCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGc                                                            <  4:122220/149‑1 (MQ=255)
                                                                                              tGTACTGAAGGAGCTCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCtgt                                                        <  4:160913/149‑1 (MQ=255)
                                                                                                 aCTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGcc                                                     >  3:109851/1‑149 (MQ=255)
                                                                                                  cTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCt                                                    <  2:275856/149‑1 (MQ=255)
                                                                                                   tGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGAt                                                                                                                                                      <  6:30332/50‑1 (MQ=255)
                                                                                                   tGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGAt                                                                                                                                                      >  5:30332/1‑50 (MQ=255)
                                                                                                   tGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCCTGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTg                                                   >  5:147212/1‑149 (MQ=255)
                                                                                                   tGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTg                                                   >  2:77814/1‑149 (MQ=255)
                                                                                                   tGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTAGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCCGTCAAGGCGAAGCGGCTTCCCGCCTTTGCCTg                                                   >  1:150676/1‑149 (MQ=255)
                                                                                                     aaGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCt                                                 <  5:142075/149‑1 (MQ=255)
                                                                                                          gcgcCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAgg                                            <  8:277551/149‑1 (MQ=255)
                                                                                                          gcgcCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAgg                                            <  2:5671/149‑1 (MQ=255)
                                                                                                            gcCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGt                                          <  7:36846/149‑1 (MQ=255)
                                                                                                               ttttCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGAGAACCGGCTTCCCGCCTGTGCCTGAACAAGGGt                                          >  7:234526/1‑146 (MQ=255)
                                                                                                                     ttCATCGAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCTGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACAt                                 <  3:246428/149‑1 (MQ=255)
                                                                                                                     ttCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACAt                                 >  4:246140/1‑149 (MQ=255)
                                                                                                                      tCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACAt                                 >  4:246767/1‑148 (MQ=255)
                                                                                                                        aTCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGAt                              >  2:294822/1‑149 (MQ=255)
                                                                                                                        aTCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGAt                              >  6:213896/1‑149 (MQ=255)
                                                                                                                         tCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGAtt                             <  1:164698/149‑1 (MQ=255)
                                                                                                                         tCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGAtt                             <  7:283591/149‑1 (MQ=255)
                                                                                                                         tCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGAtt                             <  4:109851/149‑1 (MQ=255)
                                                                                                                         tCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATCTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGAtt                             <  6:152972/149‑1 (MQ=255)
                                                                                                                           tcacccAGCAGCATGACATCGGCAAAGTAATGATGATGGGTGATGGCGGGGATCATCCCCAAGGGGTTCCAGGGTATGTGGCGAAGCACATCGTTCAAGGCGAAGCGGCTTACAGCATGTTCCTGCTCAAGGGTAAGGAACATGAttt                            >  7:56964/5‑148 (MQ=255)
                                                                                                                             ggCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGATTTGCa                         <  1:298957/149‑1 (MQ=255)
                                                                                                                              cccAGCAGCATGACATCGGCGAAGTCCTGATGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGATTTGc                          >  7:275881/2‑147 (MQ=38)
                                                                                                                              cccAGCAGCATGACATCGGCGAAGTCCTGATGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGATTTGc                          <  8:275881/146‑1 (MQ=38)
                                                                                                                                  gcagcaTGACATCGGCGAAGTCCTGATGCT‑GGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCt                                                    <  2:76482/116‑1 (MQ=37)
                                                                                                                                  gcagcaTGACATCGGCGAAGTCCTGATGCT‑GGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCt                                                    >  1:76482/1‑116 (MQ=37)
                                                                                                                                            ttCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGATTTGCAACCTGCGATAGatgt          >  3:87655/1‑149 (MQ=255)
                                                                                                                                            ttCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGATTTGCAACCTGCGATAGatgt          >  4:189153/1‑149 (MQ=255)
                                                                                                                                            ttCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGATTTGCAACCTGCGATAGat            >  7:267245/1‑147 (MQ=255)
                                                                                                                                            ttCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGATTTGCAACCTGCGATAGat            <  8:267245/147‑1 (MQ=255)
                                                                                                                                            ttCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCATGCTCAAGGGTAAGGAACATGATTTGCAACCTGCGATAGatgt          >  5:150679/1‑149 (MQ=255)
                                                                                                                                                ccGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGTGTATGTGTCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGATTTGCAACCTGCGATAGatgttatg      <  5:88813/149‑1 (MQ=255)
                                                                                                                                                    tGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGATTTGCAACCTGCGATAGATgttatgttat  <  3:246767/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GTTGGATGCCAGCACGATATGCTGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGATTTGCAACCTGCGATAGATGTTATGTTAT  >  NC_002947/3756894‑3757190

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