Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,304,474 T→A 5.4% I214N (ATT→AAT)  oppC → oligopeptide ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,304,4740TA 5.4% 253.9 / 4.9 74I214N (ATT→AAT) oppColigopeptide ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (34/36);  new base A (2/2);  total (36/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.85e-01

CTGCTGATTTTCGTGGCGATAGGCATGGTTTCCTGGCTGGATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTGGGGTGCATTGCTGAGTGATGGCG  >  NC_000913/1304326‑1304617
                                                                                                                                                    |                                                                                                                                               
ctgctgATTTTCGTGGCGATAGGCATGGTTTCCTGGCTGGATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACAt                                                                                                                                                 >  2:176854/1‑149 (MQ=255)
   ctgATTTTCGTGGCGATAGGCATGGTTTCCTGGCTGGATATGGCTCGCATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATtgt                                                                                                                                              >  3:619255/1‑149 (MQ=255)
     gATTTTCGTGGCGATAGGCATGGTTTCCTGGCTGGATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGc                                                                                                                                            <  2:93676/149‑1 (MQ=255)
     gATTTTCGTGGCGATAGGCATGGTTTCCTGGCTGGATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGc                                                                                                                                            <  4:42708/149‑1 (MQ=255)
       ttttCGTGGCGATAGGCATGGTTTCCTGGCTGGATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCg                                                                                                                                          >  3:564881/1‑149 (MQ=255)
         ttCGTGGCGATAGGCATGGTTTCCTGGCTGGATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGaa                                                                                                                                        <  3:509804/149‑1 (MQ=255)
                  aTAGGCATGGTTTCCTGGCTGGATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCgg                                                                                                                               >  4:351626/1‑149 (MQ=255)
                      gCATGGTTTCCTGGCTGGATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGAGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGCATGTACTCGGtgtg                                                                                                                           >  1:512351/1‑149 (MQ=255)
                      gCATGGTTTCCTGGCTGGATATGGCTCGTATTGTGCGCGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGtgtg                                                                                                                           >  1:466671/1‑149 (MQ=255)
                       cATGGTTTCCTGGCTGGATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTTTGCCGAATGTACTCGGTgtgg                                                                                                                          >  2:305197/1‑149 (MQ=255)
                       cATGGTTTCCTGGCTGGATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTgtgg                                                                                                                          <  2:193994/149‑1 (MQ=255)
                        aTGGTTTCCTGGCTGGATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTgtggt                                                                                                                         >  3:193709/1‑149 (MQ=255)
                           gTTTCCTGGCTGGATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTgtggtggt                                                                                                                      <  1:167594/149‑1 (MQ=255)
                             ttCCTGGCTGGATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGtggtggtgg                                                                                                                    <  4:10146/149‑1 (MQ=255)
                              tCCTGGCTGGATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTggtggtggt                                                                                                                   <  2:148848/149‑1 (MQ=255)
                                              cTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTg                                                                                                   <  1:289392/149‑1 (MQ=255)
                                              cTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTg                                                                                                   <  1:159136/149‑1 (MQ=255)
                                               tCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTgg                                                                                                  <  4:530531/149‑1 (MQ=255)
                                               tCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTgg                                                                                                  >  3:469905/1‑149 (MQ=255)
                                               tCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTgg                                                                                                  >  2:78302/1‑149 (MQ=255)
                                                     tgtgCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCa                                                                                            <  1:247552/149‑1 (MQ=255)
                                                          gTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATg                                                                                       >  1:251732/1‑149 (MQ=255)
                                                               cAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGTCCAGCATGATCct                                                                                  >  1:145778/1‑149 (MQ=255)
                                                                aaaCCCTGAGTCTGAAGCACAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCctc                                                                                 >  1:382935/1‑149 (MQ=255)
                                                                  aCCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCtt                                                                               >  4:447990/1‑149 (MQ=255)
                                                                   cccTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCttt                                                                              <  1:389934/149‑1 (MQ=255)
                                                                   cccTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCttt                                                                              <  3:130305/149‑1 (MQ=255)
                                                                   cccTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCttt                                                                              <  2:455991/149‑1 (MQ=255)
                                                                     cTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGa                                                                            >  4:394904/1‑149 (MQ=255)
                                                                      tGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGaa                                                                           >  3:200505/1‑149 (MQ=255)
                                                                      tGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGaa                                                                           >  3:138254/1‑149 (MQ=255)
                                                                          tCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCAGGTTCCCAGCATGATCCTATTTGAATCtt                                                                       <  4:112783/149‑1 (MQ=255)
                                                                                  gcAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGc                                                               <  3:529973/149‑1 (MQ=255)
                                                                                         gTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTgg                                                        >  2:488697/1‑149 (MQ=255)
                                                                                             aTTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGtt                                                    <  1:158514/149‑1 (MQ=255)
                                                                                              ttGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTg                                                   >  4:372293/1‑149 (MQ=255)
                                                                                               tGAGGTGACACAAGTTGGCGGTGTATCGACATCGGGCACTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTTGTGGTGGTCTACGCAACACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCGTGGGGTTgg                                                  <  2:509583/149‑1 (MQ=255)
                                                                                                   gcggcACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTAc                                              >  2:353655/1‑149 (MQ=255)
                                                                                                     ggcACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGc                                            <  2:512351/149‑1 (MQ=255)
                                                                                                      gcACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCa                                           >  4:313968/1‑149 (MQ=255)
                                                                                                       cacaAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCaa                                          >  3:225675/1‑149 (MQ=255)
                                                                                                        acaAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAag                                         >  3:238728/1‑149 (MQ=255)
                                                                                                        acaAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAag                                         >  4:19084/1‑149 (MQ=255)
                                                                                                         caAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAaga                                        <  2:230364/149‑1 (MQ=255)
                                                                                                         caAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAaga                                        <  3:313968/149‑1 (MQ=255)
                                                                                                         caAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAaga                                        <  1:293050/149‑1 (MQ=255)
                                                                                                            gTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGcc                                     >  2:241348/1‑149 (MQ=255)
                                                                                                              tGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGt                                   <  3:139091/149‑1 (MQ=255)
                                                                                                                gCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTa                                 <  4:106924/149‑1 (MQ=255)
                                                                                                                 cGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTaa                                <  4:564881/149‑1 (MQ=255)
                                                                                                                 cGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTaa                                <  2:251732/149‑1 (MQ=255)
                                                                                                                  ggTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAag                               <  4:318675/149‑1 (MQ=255)
                                                                                                                   gtgtATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAagc                              >  1:136085/1‑149 (MQ=255)
                                                                                                                       aTCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCt                          <  1:116010/149‑1 (MQ=255)
                                                                                                                       aTCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCt                          <  1:526128/149‑1 (MQ=255)
                                                                                                                         cGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTgg                        <  4:458682/149‑1 (MQ=255)
                                                                                                                            cATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTGGGGt                     <  4:619255/149‑1 (MQ=255)
                                                                                                                             aTCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTGGGGTg                    >  4:273823/1‑149 (MQ=255)
                                                                                                                                   cATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCttt                                                                              <  4:99312/85‑1 (MQ=255)
                                                                                                                                   cATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCttt                                                                              >  3:99312/1‑85 (MQ=255)
                                                                                                                                      tGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTGGGGTGCATTGCTGa           >  1:340182/1‑149 (MQ=255)
                                                                                                                                      tGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTGGGGTGCATTGCTGa           <  1:176854/149‑1 (MQ=255)
                                                                                                                                      tGTTACTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTGGGGTGCATTGCTGa           <  4:161757/149‑1 (MQ=255)
                                                                                                                                        ttATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTGGGGTGCATTGCTGAGt         >  2:523035/1‑149 (MQ=255)
                                                                                                                                        ttATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTGGGGTGCATTGCTGAGt         >  2:486024/1‑149 (MQ=255)
                                                                                                                                         tATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTgg                                                        <  2:224908/101‑1 (MQ=255)
                                                                                                                                         tATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTgg                                                        >  1:224908/1‑101 (MQ=255)
                                                                                                                                          aTTCGGCACAATGTGGCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTa                                                                 <  4:526677/91‑1 (MQ=255)
                                                                                                                                          aTTCGGCACAATGTGGCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTa                                                                 >  3:526677/1‑91 (MQ=255)
                                                                                                                                          aTTCGGCACAATGTGGCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTGGGGt                     <  2:267630/135‑1 (MQ=255)
                                                                                                                                          aTTCGGCACAATGTGGCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTGGGGt                     >  1:267630/1‑135 (MQ=255)
                                                                                                                                          aTTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAg                                                                >  4:577219/1‑92 (MQ=255)
                                                                                                                                          aTTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAg                                                                <  3:577219/92‑1 (MQ=255)
                                                                                                                                              gCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTGGGGTGCATTGCTGAGTGATGgc   >  3:211401/1‑149 (MQ=255)
                                                                                                                                               ccACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTGGGGTGCATTGCTGAGTGATGgcg  >  3:227091/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
CTGCTGATTTTCGTGGCGATAGGCATGGTTTCCTGGCTGGATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTGGGGTGCATTGCTGAGTGATGGCG  >  NC_000913/1304326‑1304617

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: