Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,357,665 G→A 7.7% G254G (GGG→GGA yfaD → transposase_31 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,357,6650GA 7.7% 187.3 / 3.9 77G254G (GGG→GGAyfaDtransposase_31 family protein
Reads supporting (aligned to +/- strand):  ref base G (28/42);  new base A (2/4);  total (31/46)
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

ATGACAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAG  >  NC_000913/2357517‑2357811
                                                                                                                                                     |                                                                                                                                                  
aTGACAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAggg                                                                                                                                                    >  1:827273/1‑149 (MQ=255)
   aCAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCAt                                                                                                                                                 <  2:693116/149‑1 (MQ=255)
    cagccagcTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATc                                                                                                                                                <  1:1343629/149‑1 (MQ=255)
    cagacagcTACAAACCCTGTTTATGTATCAGCGGCAAGCCGGCGATACCGCCCGGGTGACCCGTTTACTTGATGAGCCTGCCGAACCGTCCCCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATc                                                                                                                                                <  2:1509437/149‑1 (MQ=12)
       ccagcTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCaaa                                                                                                                                             >  2:659020/1‑149 (MQ=255)
       ccagcTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCaaa                                                                                                                                             >  2:829262/1‑149 (MQ=255)
       ccagcTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGACATCaaa                                                                                                                                             >  1:1405051/1‑149 (MQ=255)
          gcTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATtg                                                                                                                                          <  2:1520529/149‑1 (MQ=255)
          gcTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATtg                                                                                                                                          <  1:659020/149‑1 (MQ=255)
            tACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATAGCTCCCGTTTCGCCCGTTTTATTGAGGAGATTGCCGAACGTTCAACACTACAAAAGGAGAGATTAATTCGTATTGATGATCGGCTACGGCAGGATGGGCTTCAATTtggc                                                                                                                                        >  2:724432/1‑149 (MQ=255)
             aCAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATtggct                                                                                                                                       <  1:402702/149‑1 (MQ=255)
             aCAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTGTATTGAGGAGATTGCCGAACGGTCACCACTACAAAATGAGCGATTAATGACTATTGCTGAACGGCTTCGGCAGGAAGGGCATGAAATtggcg                                                                                                                                       >  2:1488233/1‑148 (MQ=38)
              cAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATtggctg                                                                                                                                      >  2:345451/1‑149 (MQ=255)
               aaaCACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCCCCACTACAAAAGGAGAGATTAATGACTATTGCTTAACGGCTACGGCAGGAAGGGCATCAAATtggctgg                                                                                                                                     >  1:688219/1‑149 (MQ=255)
               aaaCACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATAAAATtggctgg                                                                                                                                     >  2:646962/1‑149 (MQ=255)
                 acacTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTGCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCa                                                                                                                                   <  2:1319866/149‑1 (MQ=255)
                 acacTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTGCACCCGTTTTATTGAGGAGATTGCAGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCa                                                                                                                                   >  2:563038/1‑149 (MQ=255)
                  cacTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAg                                                                                                                                  <  1:43463/149‑1 (MQ=255)
                    cTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGATTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGa                                                                                                                                <  2:827273/149‑1 (MQ=255)
                    cTGTTTAATTATCAGCTGCAATGCGGCGATACCTCCCGTTTCACCGCGGGTTTTGAGGAGATTGCCGAACGTTCACCAATACAATTGGAGAGATTAATCACTATTGCAGTACGGCTACGGCAGGAATTGCATCAAATTGGCTGGCAGGa                                                                                                                                <  2:483473/149‑1 (MQ=25)
                      gTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGAATAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAg                                                                                                                              <  1:1562643/149‑1 (MQ=255)
                       tttAATTATCTGCTGCAATGCTTCGATACCTCCCGGTTCACCCGTTTTATTGAGGAGATTGACGAACGTTCACCACTACAAAAGGAGAGATTAATGCCTATTACTGCACGGCTACGGCAGGAAGGACCTAAAATTGGCTGGCAGCaatg                                                                                                                             >  2:1382618/1‑147 (MQ=255)
                       tttAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGACAGATTAATGACTATAGCTGACCTGATACGGCAGGATGGGAATCAAATTGGCTGGCAGGAAgt                                                                                                                             >  2:312185/1‑148 (MQ=38)
                        ttAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGt                                                                                                                            <  2:16572/149‑1 (MQ=255)
                          aaTTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTaa                                                                                                                          >  1:209044/1‑148 (MQ=255)
                            ttATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGGTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGc                                                                                                                        <  2:1072286/149‑1 (MQ=255)
                              aTCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCAt                                                                                                                      >  2:636142/1‑149 (MQ=255)
                               tCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGACCGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATg                                                                                                                     <  1:275017/149‑1 (MQ=255)
                                  gctgcAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGATGGTATGCATGAAc                                                                                                                  >  2:1330571/1‑149 (MQ=255)
                                    tgcAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCAACACTACAAAATGAGAGATTAATGACTATTGATGAACGGCTACGGCACGAAGGGCATCAAATTTGATGTCAGTAAGGTATGCATGAACaa                                                                                                                >  2:113934/1‑149 (MQ=255)
                                         tGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCAt                                                                                                           >  2:893432/1‑149 (MQ=255)
                                                aTACCTCCCGTTTCACCGGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAAtt                                                                                                    <  2:1218778/149‑1 (MQ=255)
                                                       ccGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTgc                                                                                             <  2:816560/149‑1 (MQ=255)
                                                            tCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATg                                                                                        <  1:893432/149‑1 (MQ=255)
                                                            tCACCCGTTTTAATGAGGAGATTGCCGAACGATCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATg                                                                                        <  1:704866/149‑1 (MQ=255)
                                                               cccGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGAGTAATGACTATTGCGCAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTg                                                                                     <  2:478582/149‑1 (MQ=255)
                                                               cccGATTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAATCCAATAAAAATGCTATGCGCAttcag                                                                                     >  1:1090833/1‑145 (MQ=25)
                                                                   gTTTTATTGAGGAGATTGCCGATCGTTCACCTCTACAAAAGGATAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAg                                                                                  <  1:422810/149‑1 (MQ=255)
                                                                   gTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAg                                                                                  <  2:358969/149‑1 (MQ=255)
                                                                      ttATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAg                                                                               <  1:14906/149‑1 (MQ=255)
                                                                        aTTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAggg                                                                             <  2:1184621/149‑1 (MQ=255)
                                                                        aTTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAggg                                                                             <  1:636142/149‑1 (MQ=255)
                                                                         ttGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTTCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTaaat                                                                                                                        <  2:330397/105‑4 (MQ=255)
                                                                         ttGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTaaat                                                                                                                        >  1:330397/1‑102 (MQ=255)
                                                                          tGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCt                                                                           >  1:904135/1‑149 (MQ=255)
                                                                           gaggagATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGAGCATGCTGGAGCAGGGCtt                                                                          >  1:543917/1‑149 (MQ=255)
                                                                                gaTTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGAAGGAAGGTATGCATGAACAGGGCATTTAAATTGCTTTTCGCATGCGGGAGCTGGGcggtgatc                                                                     >  2:309839/1‑142 (MQ=25)
                                                                                  ttGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGt                                                                   <  2:904135/149‑1 (MQ=255)
                                                                                    gCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTga                                                                 <  1:345451/149‑1 (MQ=255)
                                                                                     ccGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAGATAGAACAAGAGATTAAAATTGCAATGCGCAATCATTTGCAGGTCGTTGAACGTgcc                                                                >  1:1450496/1‑147 (MQ=17)
                                                                                         aCGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTa                                                                                                                           >  1:1488176/1‑86 (MQ=255)
                                                                                         aCGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTa                                                                                                                           <  2:1488176/86‑1 (MQ=255)
                                                                                                    tACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCAAGAACAAGCCATTAAAATTGCTTTGAGAATTCTGGAGCAGGGCGTTGAACGTGAGATTGTGCTTGCgccc                                                 >  2:962750/1‑146 (MQ=25)
                                                                                                      cAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAAc                                               <  2:415623/149‑1 (MQ=255)
                                                                                                      aaaaaGGAGAGATTATTGACTATTGATGAACGGCTACGGCAGGAAGGGCATCAAATTCGCGTGCAGGAAGGTATGCATGAACAAGCCATTATAATTGCTTTGCGCATGCTGGACCAGGGCGTTGAACGTGAGATTGTGCTGGCGACAAc                                               <  2:320623/148‑1 (MQ=25)
                                                                                                        aaaGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAAccc                                             >  1:765236/1‑149 (MQ=255)
                                                                                                         aaGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCGATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGATCGTGACATGGTGCTCGCGGCCACCCa                                            >  2:1204950/1‑149 (MQ=18)
                                                                                                          acgagagaTTAATGTCTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTTTGCATGAACAAGCCATTAAACTTGCTTTGCGCATCCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCaa                                           <  1:117117/147‑1 (MQ=255)
                                                                                                          aGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCaa                                           <  2:821162/149‑1 (MQ=255)
                                                                                                          aGGAGAGATAAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTTGCAGGAAGGTATGCATGAACAAGCCATTAAAATTCCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGATAACCCaa                                           <  2:350971/149‑1 (MQ=255)
                                                                                                            gagagaTTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACt                                         <  1:563038/149‑1 (MQ=255)
                                                                                                              gagaTTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCa                                       >  1:611057/1‑149 (MQ=255)
                                                                                                              gagaTTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGAATGAACAATACATTAAAATTGCTTTGAGAATTCTGTAGCAGGGCATTGAACGTGACATTGTGCTGGCGCCCACCCAGCtag                                       >  2:38765/1‑147 (MQ=21)
                                                                                                                 aTTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTg                                    >  1:769451/1‑149 (MQ=255)
                                                                                                                 aTTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTg                                    <  1:312185/149‑1 (MQ=255)
                                                                                                                  ttAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGa                                   <  1:1333865/149‑1 (MQ=255)
                                                                                                                     aTGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGc                                <  2:1356918/149‑1 (MQ=255)
                                                                                                                       gACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACACACTGATGCTg                              <  1:90683/149‑1 (MQ=255)
                                                                                                                       gACTATTCCTGAACGGCGACGGCAGGAAGGGCATCAAATTGGCTCGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTg                              <  1:85622/149‑1 (MQ=255)
                                                                                                                                tGAACGGCTACGGCAGGAAGGTCATCAAATTGGCTGGAAGGAAGGTATGAATGAACAAGACATTAAAATTGATTTGATAATGCTGGAACAGGGCTTTGATCGTGACAATGTGCTCGCGTCCACCCAGCTAAGAGAAGCCGATATAccac                     >  2:1233211/1‑147 (MQ=9)
                                                                                                                                     ggCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGt                <  1:1368827/149‑1 (MQ=255)
                                                                                                                                        tACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCAt             <  1:646962/149‑1 (MQ=255)
                                                                                                                                            gCAGGAAGGGCATCAATTTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGATATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTaaa         <  1:829262/149‑1 (MQ=255)
                                                                                                                                             cAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATTCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTaaaa        <  1:387514/149‑1 (MQ=255)
                                                                                                                                             cAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTaaaa        <  2:765236/149‑1 (MQ=255)
                                                                                                                                              aGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAg       <  2:769451/149‑1 (MQ=255)
                                                                                                                                                gAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGga     >  1:1571529/1‑149 (MQ=255)
                                                                                                                                                gAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGga     <  1:1330571/149‑1 (MQ=255)
                                                                                                                                                   tggCATCAAATTGGCTGGCAGGAAGGTATGCATTAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGgagag  <  2:688219/148‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
ATGACAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAG  >  NC_000913/2357517‑2357811

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