Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,304,479 C→G 5.6% P216A (CCG→GCG)  oppC → oligopeptide ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,304,4790CG 5.6% 245.7 / 6.4 72P216A (CCG→GCG) oppColigopeptide ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (33/35);  new base G (2/2);  total (35/37)
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.68e-02

GATTTTCGTGGCGATAGGCATGGTTTCCTGGCTGGATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTGGGGTGCATTGCTGAGTGATGGCGCGAACTCGAT  >  NC_000913/1304331‑1304627
                                                                                                                                                    |                                                                                                                                                    
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:159136/149‑1 (MQ=255)
                                         cTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTg                                                                                                             <  1:289392/149‑1 (MQ=255)
                                          tCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTgg                                                                                                            >  3:469905/1‑149 (MQ=255)
                                          tCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTgg                                                                                                            >  2:78302/1‑149 (MQ=255)
                                          tCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTgg                                                                                                            <  4:530531/149‑1 (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                                                                                        <  2:455991/149‑1 (MQ=255)
                                                              cccTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCttt                                                                                        <  1:389934/149‑1 (MQ=255)
                                                              cccTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCttt                                                                                        <  3:130305/149‑1 (MQ=255)
                                                                cTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGa                                                                                      >  4:394904/1‑149 (MQ=255)
                                                                 tGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGaa                                                                                     >  3:138254/1‑149 (MQ=255)
                                                                 tGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGaa                                                                                     >  3:200505/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                                                   >  4:19084/1‑149 (MQ=255)
                                                                                                   acaAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAag                                                   >  3:238728/1‑149 (MQ=255)
                                                                                                    caAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAaga                                                  <  2:230364/149‑1 (MQ=255)
                                                                                                    caAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAaga                                                  <  1:293050/149‑1 (MQ=255)
                                                                                                    caAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAaga                                                  <  3:313968/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                                          <  2:251732/149‑1 (MQ=255)
                                                                                                            cGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTaa                                          <  4:564881/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:176854/149‑1 (MQ=255)
                                                                                                                                 tGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTGGGGTGCATTGCTGa                     >  1:340182/1‑149 (MQ=255)
                                                                                                                                 tGTTACTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTGGGGTGCATTGCTGa                     <  4:161757/149‑1 (MQ=255)
                                                                                                                                   ttATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTGGGGTGCATTGCTGAGt                   >  2:486024/1‑149 (MQ=255)
                                                                                                                                   ttATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTGGGGTGCATTGCTGAGt                   >  2:523035/1‑149 (MQ=255)
                                                                                                                                    tATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTgg                                                                  <  2:224908/101‑1 (MQ=255)
                                                                                                                                    tATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTgg                                                                  >  1:224908/1‑101 (MQ=255)
                                                                                                                                     aTTCGGCACAATGTGGCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTa                                                                           >  3:526677/1‑91 (MQ=255)
                                                                                                                                     aTTCGGCACAATGTGGCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTa                                                                           <  4:526677/91‑1 (MQ=255)
                                                                                                                                     aTTCGGCACAATGTGGCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTGGGGt                               >  1:267630/1‑135 (MQ=255)
                                                                                                                                     aTTCGGCACAATGTGGCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTGGGGt                               <  2:267630/135‑1 (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)
                                                                                                                                                   gCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGAACCTCTTTTAATCTTTCCTTAGCTTCCAGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTGGGGTGCATTGCTGAGTGATGGCGCGAACTCGa   <  4:193709/149‑1 (MQ=255)
                                                                                                                                                    ccGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTGGGGTGCATTGCTGAGTGATGGCGCGAACTCGAt  <  1:488697/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GATTTTCGTGGCGATAGGCATGGTTTCCTGGCTGGATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTGGGGTGCATTGCTGAGTGATGGCGCGAACTCGAT  >  NC_000913/1304331‑1304627

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: