Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,357,653 A→G 8.9% L250L (CTA→CTG yfaD → transposase_31 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,357,6530AG 8.9% 208.7 / 9.3 81L250L (CTA→CTGyfaDtransposase_31 family protein
Reads supporting (aligned to +/- strand):  ref base A (33/39);  new base G (2/5);  total (36/45)
Fisher's exact test for biased strand distribution p-value = 4.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.81e-01

TCACTAATGACAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCAT  >  NC_000913/2357511‑2357800
                                                                                                                                               |                                                                                                                                                   
tCACTAATGACAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCTCCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAg                                                                                                                                               >  1:415623/1‑149 (MQ=255)
tCACTAATGACAGCCAGCTACAAACACTGTTTAATTATATGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATAAATGACTATTGCTGAACGGCTACGGCAg                                                                                                                                               >  1:320623/1‑149 (MQ=255)
   cTAATGACAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCCCCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAggca                                                                                                                                            >  1:1047950/1‑147 (MQ=255)
     aaTGACAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGGTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGCGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAgg                                                                                                                                          <  2:839737/149‑1 (MQ=255)
      aTGACAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAggg                                                                                                                                         >  1:827273/1‑149 (MQ=255)
       tGACAGCAATCTACAAACCCTGTTTAATTATCTGCGGCATTCCGGCGATACCTCCCGTGTGACCCGTTAAATTGAGGAGATTGCCGAACGTTCACCTCTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAggac                                                                                                                                        <  2:248074/149‑3 (MQ=25)
         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)
                                                                                            cGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAgga                                                                                                                                         <  2:1148561/64‑2 (MQ=255)
                                                                                            cGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAgga                                                                                                                                         >  1:1148561/1‑63 (MQ=255)
                                                                                               aCGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTa                                                                                                                >  1:1488176/1‑86 (MQ=255)
                                                                                               aCGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTa                                                                                                                <  2:1488176/86‑1 (MQ=255)
                                                                                                     accacTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAg                                                                                                                                           <  2:344268/53‑1 (MQ=255)
                                                                                                     accacTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAg                                                                                                                                           >  1:344268/1‑53 (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)
                                                                                                                                               |                                                                                                                                                   
TCACTAATGACAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCAT  >  NC_000913/2357511‑2357800

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: