Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,263,409 C→T 100% intergenic (+448/+166) yhaC → / ← rnpB pentapetide repeats‑related protein/RNase P, M1 RNA component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,263,4090CT84.2% 229.0 / 41.5 113intergenic (+448/+166)yhaC/rnpBpentapetide repeats‑related protein/RNase P, M1 RNA component
Reads supporting (aligned to +/- strand):  ref base C (5/12);  new base T (34/58);  total (42/71)
Fisher's exact test for biased strand distribution p-value = 7.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Polymorphic base substitution creates a homopolymer stretch.

Only 100 of 118 total aligned reads displayed.

TTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGC  >  CP009273/3263265‑3263551
                                                                                                                                                |                                                                                                                                              
tttGGGTATTGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATtttt                                                                                                                                             <  2:1626354/148‑1 (MQ=255)
tttGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAACAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATtttt                                                                                                                                             <  1:594355/148‑1 (MQ=255)
tttGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATtttt                                                                                                                                             <  1:594313/148‑1 (MQ=37)
tttGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATtttt                                                                                                                                             <  2:651137/148‑1 (MQ=37)
        taGTGTCGTGGACAGTCATTCATCTTTATGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAACCCGCGTTCTTTTGTGGATATTGTCGTGGAAAGTCATAAAATTTTATGCCCCt                                                                                                                                     >  1:1886964/1‑148 (MQ=11)
         aGTGTCGTGGACAGTCAGTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTc                                                                                                                                    <  2:2111289/148‑1 (MQ=255)
             tCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCaaa                                                                                                                                <  1:2217100/148‑1 (MQ=37)
                 ggACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCa                                                                                                                            <  1:1182155/148‑1 (MQ=37)
                 ggACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCa                                                                                                                            <  1:1600644/148‑1 (MQ=37)
                 ggACAGTCATTCATCTTTATGCCACTCCAAAATTAAAAACCCGCCGTAGCGGTTTATTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCa                                                                                                                            <  1:1601117/148‑1 (MQ=11)
                   aCAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCaaa                                                                                                                          <  2:2989414/148‑1 (MQ=37)
                   aCAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCaaa                                                                                                                          <  2:408386/148‑1 (MQ=37)
                   aCAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCaaa                                                                                                                          <  1:1797055/148‑1 (MQ=37)
                    cAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCaaaa                                                                                                                         <  1:1146747/148‑1 (MQ=37)
                    cAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCaaaa                                                                                                                         <  2:1820666/148‑1 (MQ=37)
                    cAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCaaaa                                                                                                                         <  2:561633/148‑1 (MQ=37)
                     aGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCaaaaa                                                                                                                        <  2:1306299/148‑1 (MQ=37)
                              tCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTAGGGTAAAGAGTAGTGGAAAGTCATTCAATTTTAAGCACCTGCAAAACCAAAAACCCCCCaaa                                                                                                               >  2:148264/1‑148 (MQ=11)
                               cTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAg                                                                                                              >  2:1156221/1‑148 (MQ=37)
                               cTTTCTGCCCCACCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAg                                                                                                              <  1:2789468/148‑1 (MQ=255)
                                 ttCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCg                                                                                                            <  1:166201/148‑1 (MQ=37)
                                 ttCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCg                                                                                                            <  1:166170/148‑1 (MQ=37)
                                   cTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCggg                                                                                                          <  1:2767524/148‑1 (MQ=37)
                                     gCCCCTCCACAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACCGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGtt                                                                                                        <  2:521871/148‑1 (MQ=12)
                                          tCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTAc                                                                                                   >  1:2137497/1‑148 (MQ=37)
                                          tCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTAc                                                                                                   >  1:2134554/1‑148 (MQ=37)
                                              acaggaaaaaCCCGCCGCAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTaa                                                                                               <  2:1526359/143‑1 (MQ=255)
                                               aactAAAAACCCGCCGATGCGTGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGAATGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTaaa                                                                                              <  1:406441/145‑1 (MQ=11)
                                               aaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTaaa                                                                                              <  2:778228/148‑1 (MQ=37)
                                               aaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCAGTCATGTTTCTGGCCCGCCAAAAGCAAAAACCCGCCGAAGCGGGTGTTTACGTaaa                                                                                              >  2:2998435/1‑148 (MQ=11)
                                                aGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAc                                                                                             >  2:1788151/1‑148 (MQ=38)
                                                 gTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCACTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAcc                                                                                            <  2:2148357/148‑1 (MQ=255)
                                                  tAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCa                                                                                           >  2:2812289/1‑148 (MQ=38)
                                                  tAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCa                                                                                           <  1:1154773/148‑1 (MQ=38)
                                                  tAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCa                                                                                           <  1:1154664/148‑1 (MQ=38)
                                                  tAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTATGCCCCTCCAAAAGCAAAAAACCGCCGAAGCGGGTTTTTAAGTAAAACa                                                                                           >  2:945677/1‑148 (MQ=11)
                                                      aaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTg                                                                                       <  2:2841431/148‑1 (MQ=37)
                                                       aCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGGCCTTCATCGTTGTGCCCCGCCCAAAGCCAAAACCCGCCGCAGCGGGTTTTTAAGGAAACCAGGtta                                                                                      >  1:2323223/1‑146 (MQ=11)
                                                         ccgccgAAGCGGTTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGaaa                                                                                    <  2:3106046/148‑1 (MQ=255)
                                                          cgccgAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTAATTTTTCTGCCCCTCCCAAATCAAAAACCCGCCGAAGCGGGTTTTTTCGTAAAACAGGTGAAAc                                                                                   >  2:1149668/1‑148 (MQ=11)
                                                          cgccgAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGGATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACACGTGaaaa                                                                                   >  2:59220/1‑147 (MQ=255)
                                                            ccgAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTg                                                                                 <  2:2134455/148‑1 (MQ=37)
                                                            ccgAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTg                                                                                 <  2:2137497/148‑1 (MQ=37)
                                                            ccgAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTg                                                                                 <  2:2134554/148‑1 (MQ=37)
                                                             cttatcGGTTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTTTAGTTTCGTGGACATTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGa                                                                                <  2:3244716/143‑1 (MQ=11)
                                                             cgAATCGGGTTATTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGa                                                                                <  2:1437624/148‑1 (MQ=255)
                                                             cgAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGa                                                                                <  2:3246170/148‑1 (MQ=37)
                                                              gAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGGGTCGTGGACAGGCATTCATATTTCTGCCCCAACAAAAGCAAATAACCGCCGAAGGGGGTTTTTACGTAAAAAAGGTGAAACTGAc                                                                               >  1:140698/1‑148 (MQ=11)
                                                                   gggTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATa                                                                          <  1:3086588/148‑1 (MQ=37)
                                                                     gTTTTTACGTACAACAGGTTAAACTGACCGATAAGCCGATTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCACAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAg                                                                        <  2:923593/148‑1 (MQ=11)
                                                                     gTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAg                                                                        <  2:3100171/148‑1 (MQ=37)
                                                                     gTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAg                                                                        <  2:1350482/148‑1 (MQ=17)
                                                                      tttttACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGc                                                                       <  2:2109887/148‑1 (MQ=255)
                                                                        tttACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGGGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCg                                                                     <  1:914737/148‑1 (MQ=255)
                                                                         ttACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGc                                                                    <  1:1046401/148‑1 (MQ=38)
                                                                         ttACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTTAAACTGACCGATAAGCCGc                                                                    <  2:239716/148‑1 (MQ=38)
                                                                              tacaaCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTAGAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTTAAACTGACCGATAAGCCGctttct                                                               <  1:2719551/145‑1 (MQ=255)
                                                                              tAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGctttct                                                               <  1:2471147/148‑1 (MQ=255)
                                                                              tAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTTAAACTGACCGATAAGCCGctttct                                                               <  1:2719517/148‑1 (MQ=255)
                                                                               aaaaCAGGTGAAACTGACCTATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGctttctt                                                              <  1:2763601/148‑1 (MQ=255)
                                                                                 aaCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGAACGATAAGCCGGGTTCtttt                                                            >  2:3122993/1‑148 (MQ=255)
                                                                                 aaCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGGTGTTACGTAAAACAGGTGAAACTGTCCGATAAGGCGAGTTCtttt                                                            >  2:1049549/1‑148 (MQ=11)
                                                                                  aCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAAACCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTg                                                           >  1:2475064/1‑148 (MQ=12)
                                                                                  aCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTATGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTTCGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTg                                                           >  1:2961026/1‑148 (MQ=255)
                                                                                   cAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCGTTCTGTTgg                                                          >  2:2448103/1‑148 (MQ=255)
                                                                                   cAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGt                                                                                                         <  1:2853042/101‑1 (MQ=32)
                                                                                    cgctcaaaCTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCGGCCCCTCCAAAATCAAAAACCCGCCGCAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTggg                                                         <  2:667824/143‑1 (MQ=11)
                                                                                    aGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTggg                                                         >  2:911497/1‑148 (MQ=37)
                                                                                        aaaaCTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGtata                                                     <  1:1395645/147‑1 (MQ=18)
                                                                                         aaaCTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAg                                                    <  2:2926039/148‑1 (MQ=21)
                                                                                             tGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCAGTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGcgtc                                                <  2:3301348/148‑1 (MQ=255)
                                                                                              gACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGcgtcg                                               <  1:1287451/148‑1 (MQ=37)
                                                                                                ccGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACCGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAg                                                                                                              >  2:323640/1‑83 (MQ=18)
                                                                                                ccGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAg                                                                                                              <  1:323640/83‑1 (MQ=33)
                                                                                                  gATAAGCCGCTTTCTTTTGGGTATAGTGTCGTTGAAAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAATCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGa                                           <  1:931328/148‑1 (MQ=255)
                                                                                                  gATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGa                                           <  1:930890/148‑1 (MQ=37)
                                                                                                  gATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCCGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTGTTTTGGGGATAGCGTCGTggg                                           >  2:2777248/1‑147 (MQ=255)
                                                                                                   aTAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCGTTCTTTTGGGTCTAGagttggggcc                                          >  2:3083705/1‑141 (MQ=11)
                                                                                                   aTAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTAATTTTTATTCCCCAACAAAAGCAAAAAACCCCCCAAGCGGGTTTTTACGTAAAACAGGTGAAACAGAACCATAAGCCGCGTTCTTTTGGGTATAGCGTCGTggcc                                          >  1:1220168/1‑146 (MQ=11)
                                                                                                   aTAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCAATGAATTTTCTTCCCCACCAAAAGCAAAAACCCCCCGAAGCGGGTTTTTACGTAAACCAGGTGTAACAGACCGATAAGCCCCGTTGTTTTTGGGGTAGCGtggtgggc                                          >  1:3187272/1‑146 (MQ=11)
                                                                                                    tAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTGATTTTTCTGCCCCACCAAAAGCAAAAAACCCCCGAAGCGGGGTTTTAAGTAAACCACGTGAAACTGACCGATAATCCCCTTTCGTTTGGGGATAGCGGCggggccc                                         >  1:2411681/1‑145 (MQ=11)
                                                                                                      aGCCGCTTTCTTTTGGGTATAGTGTCGTGGACCGTCATTCATTTTTATGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTAAGTAAACCACGTGAAACAGAACGAAAAGCCGCGTTATTTTGGGTATAGCGTCGGGGGCAGt                                       >  1:2099698/1‑148 (MQ=11)
                                                                                                        ccGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTATTCCCCTCCAAAAGCAAAAACCCGCCGAACCGGGTTTTTACCTAAACCAGGTGAAAATGACCGATAAGCCCCGTTCTTTTGGGTAACGCGTCGTGGACAGtca                                     >  2:1016146/1‑148 (MQ=11)
                                                                                                        ccGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTAATTTTTCTGCCCCTCCAAAAGCAAAAACCCCCCCAAGCGGGTTTTTACGTAAACCAGGTGAAAATGACCCATAAGCCGCATTCTTTTGGGTATAGCGTCGTGGAAAGgca                                     >  1:2185092/1‑148 (MQ=11)
                                                                                                           ctttctttTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGtcattc                                  <  2:2459375/148‑1 (MQ=37)
                                                                                                           ctttctttTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGtcattc                                  <  1:2812289/148‑1 (MQ=37)
                                                                                                            tttctttTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGtcaataa                                 >  2:2419059/1‑144 (MQ=255)
                                                                                                                     ggTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACCGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGGCAGTCATCTTTGTGc                        >  1:2638136/1‑148 (MQ=11)
                                                                                                                       tataGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGccc                      <  1:1149668/148‑1 (MQ=37)
                                                                                                                               cGTGGACAGTCATTCATTTTTCTGCCCCTCCAATATCATAAACCCGCCGTATCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCaaaa              <  2:2360166/148‑1 (MQ=11)
                                                                                                                               cGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCaaaa              >  1:2694047/1‑148 (MQ=37)
                                                                                                                                   gACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCaa          >  1:222871/1‑148 (MQ=37)
                                                                                                                                   gACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCaa          <  2:2694047/148‑1 (MQ=37)
                                                                                                                                   gACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCGCCAAAAGCaa          >  1:222878/1‑148 (MQ=255)
                                                                                                                                   gACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGGATAGCGTCGTGGACAGTCCTTCAACGTTCTGCCCCTACAAAAGccc          >  1:2965019/1‑146 (MQ=11)
                                                                                                                                       gTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGAAAGTCATTCAAAGTTCTGCCCCACCAAAAGCAAAAAc      >  2:856801/1‑148 (MQ=11)
                                                                                                                                        tcattcatTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAAcc     <  1:911497/148‑1 (MQ=37)
                                                                                                                                         cattcattTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTAATCCCCGCCAAAAACAAAAAacc    >  2:688652/1‑148 (MQ=11)
                                                                                                                                         cattcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTACACCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAAccc    >  2:1109708/1‑148 (MQ=255)
                                                                                                                                           ttcatTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTATGCCCCTCCAAAAGCAAAAACccgc  >  1:2283832/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                              
TTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGC  >  CP009273/3263265‑3263551

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