Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,263,466 G→T 6.9% intergenic (+505/+109) 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,4660GT 6.9% 347.9 / 8.3 138intergenic (+505/+109)yhaC/rnpBpentapetide repeats‑related protein/RNase P, M1 RNA component
Reads supporting (aligned to +/- strand):  ref base G (76/49);  new base T (4/7);  total (80/58)
Fisher's exact test for biased strand distribution p-value = 1.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 140 total aligned reads displayed.

CCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACA  >  CP009273/3263322‑3263613
                                                                                                                                                |                                                                                                                                                   
ccgccgAAGCGGGTTTTTACGTAACACAGGTGAAACTGACCGATCAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAAc                                                                                                                                                 <  2:3273354/149‑1 (MQ=255)
ccgccgAAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAAc                                                                                                                                                 <  2:1027901/149‑1 (MQ=255)
 cgccgAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACt                                                                                                                                                >  1:2396576/1‑149 (MQ=17)
    cgcaGCGGGTTTTTACGTACAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGAc                                                                                                                                             <  2:2274179/146‑1 (MQ=255)
     gacgCGGGTTATTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGAcc                                                                                                                                            <  2:3268452/146‑1 (MQ=255)
      aaGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCAATCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCg                                                                                                                                           >  1:1650347/1‑149 (MQ=255)
       aGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGa                                                                                                                                          <  1:2102970/149‑1 (MQ=37)
         cGGCTTATTACTTACAACAGGTCAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATa                                                                                                                                        <  2:34115/149‑1 (MQ=11)
           ggTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGGCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGGTTTTACGTTAACCAGGTGAAACTGACCGATAAg                                                                                                                                      >  2:2890447/1‑149 (MQ=255)
            tttcttACGTACAACAGGTGAAACTGACCGATACGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGc                                                                                                                                     <  2:3122466/145‑1 (MQ=255)
            gTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGc                                                                                                                                     <  1:1855570/149‑1 (MQ=37)
            gTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGc                                                                                                                                     <  1:3474301/149‑1 (MQ=37)
              ttttACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCCTTAATTTTTGTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACATACCGATAAGCCg                                                                                                                                   >  2:1424090/1‑149 (MQ=11)
                    ttAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGTTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAATCAAAAACCCGCCGAAGCTGGTTTTTACGTAAACCAGGTTAAACTGACCGATAAGCCGctttct                                                                                                                             <  1:129715/148‑1 (MQ=14)
                         aCAGTTTAAACTGACCGATAAGCCGCTTTCTTTTTTTTATAGTGTCGTGGACAGTCATTAATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg                                                                                                                        <  1:501912/149‑1 (MQ=11)
                         aCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg                                                                                                                        <  1:501853/149‑1 (MQ=38)
                           aGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACACTCATTCATTTTTGTGCCCCACCAAAACCAAAAACACGCCCAATCGGGTTTTTACGTAAACCAGGTGAAACTGACCGTTTATCCGCTTGATTTTgggg                                                                                                                      >  1:2274943/1‑148 (MQ=11)
                              tcaaaCTGAGATATAAGCCGCTTACTGTTGGCTATAGTGACGTTCCCAGTCATTCATTTTTATCCCACTCCACAAGCAAAAACCCGCCGCATCGGGTTTTTACGTAAACCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGtata                                                                                                                   <  2:2014152/147‑1 (MQ=11)
                               gAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAg                                                                                                                  <  1:140420/149‑1 (MQ=25)
                                 aaCTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCGTTCTTTTggggacagtc                                                                                                                >  2:550178/1‑142 (MQ=1)
                                  aCTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGcgt                                                                                                               <  2:1097275/149‑1 (MQ=37)
                                  aCTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCAGCCCCTCCAAAAGCAAAAAACCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTGTTTGGGATAGcgt                                                                                                               >  2:1613318/1‑149 (MQ=11)
                                        ctacaaGCCGCGTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGa                                                                                                         <  2:2654568/145‑1 (MQ=255)
                                        cGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGa                                                                                                         <  2:2654595/149‑1 (MQ=37)
                                         gATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGCCAGTCAATCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGAc                                                                                                        >  2:918537/1‑149 (MQ=255)
                                          aTAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGGATAGCGTAGTGGACa                                                                                                       >  2:197499/1‑149 (MQ=255)
                                          aTAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAAATGACCGATAAGCCGCTTTCTTTTGGGTATAGAGACGTGGCCa                                                                                                       >  2:2273089/1‑149 (MQ=11)
                                          aTAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTAATTTTTCTGCCCCTCCAAAATCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAAATGACCGATAAGCCGCGTTCTTTAGGGTAAAGCGTAGGGGaaa                                                                                                       >  1:2010371/1‑147 (MQ=11)
                                           tAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACACTCAATAATTTTTATGCCCCTCCAAAATCAAAAAACCCCCGAAGCGGGTTTTTTCGTAAAACAGGTGAAACAGACCGCTAAGCCGCTTTCTTTTGGGGATAGCGACGTggcccg                                                                                                      >  1:1647449/1‑145 (MQ=11)
                                             aGCCGGTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGtc                                                                                                    <  1:700055/149‑1 (MQ=255)
                                             aGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTTCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGCTTTTTACGTAAACCAGGTGAAACTTACCGATAAGCCGCGTTCTTTTGGGTATAGAGTCGTGGGCAGtc                                                                                                    >  2:279084/1‑149 (MQ=11)
                                             aGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAATCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTATTTTGGGTATAGCGTCGTGGAAAGtc                                                                                                    >  1:3775993/1‑149 (MQ=255)
                                             aGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATATTTCTGCCCCTCCAACAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGtc                                                                                                    <  2:2550509/149‑1 (MQ=255)
                                             aGCCGCTTTCTTTTGGGTATAGTGCCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGtc                                                                                                    <  1:700100/149‑1 (MQ=255)
                                             aGCCGCTTTCTTTTGGGTATAGTGACGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGtc                                                                                                    <  1:700144/149‑1 (MQ=12)
                                                cGCTATCTTTTGGGTATTGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGtcatt                                                                                                 <  2:3857891/149‑1 (MQ=255)
                                                 gCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGtcaatc                                                                                                >  2:1503670/1‑149 (MQ=255)
                                                 gCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTAATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCGTTCTTTTGGGTATAGCGTCGTGGACAGtcatta                                                                                                >  2:3552388/1‑148 (MQ=255)
                                                   tttctttTGTGTATAGTGTCGTGGACAGTCATTCAATTTTCGGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGtcattcat                                                                                              <  2:940231/149‑1 (MQ=255)
                                                       ttttGGGTATAGTGTCGTGGACAGTCATTAATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCttt                                                                                          <  1:4058014/149‑1 (MQ=255)
                                                        ttagggTTTCGTGTCGTCGACAGTCTATCTTTTTTCTGCCCCTCCACAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTc                                                                                         <  2:2177627/146‑1 (MQ=11)
                                                             gTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGGGGGCAGTGCTTCATCGTTGTGccc                                                                                    >  1:1920985/1‑149 (MQ=11)
                                                              tataGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATATTTATTcccc                                                                                   >  2:1019467/1‑149 (MQ=11)
                                                                taGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTGTTGGGTATAGCGTCGTGGACAGTCAATCATCATTATGccccac                                                                                 >  2:2622966/1‑147 (MQ=11)
                                                                              gTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTTCCCCTCCAAAAGCAAAAAcc                                                                   >  2:3115446/1‑149 (MQ=255)
                                                                              gTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAAcc                                                                   >  1:2070754/1‑149 (MQ=37)
                                                                                 attcatTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGTGTATAGCGTCGTGGACAGTCCTTCATCGTTCAGCCCCGCCAAAAGCAAAAACccgc                                                                >  2:1031538/1‑149 (MQ=11)
                                                                                 attcatTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACccgc                                                                <  1:1439795/149‑1 (MQ=37)
                                                                                 attcatTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTAGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACccgc                                                                >  2:519539/1‑149 (MQ=255)
                                                                                 attcatTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTGTTCGGGATAGCGGCGTGGACAGGCAATCATATTACTTCCCCACCATAAGAAAAAACccga                                                                >  2:2838569/1‑148 (MQ=11)
                                                                                  ttcatTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTGTTGGGTATAGCGTCGTGGACAGTCATTCATATTTCTGCCCCTCCAAAAGCAAAAACccgcc                                                               >  2:3311418/1‑149 (MQ=255)
                                                                                    catTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGa                                                             >  1:2826521/1‑149 (MQ=37)
                                                                                     atttttCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGaa                                                            >  2:778917/4‑149 (MQ=37)
                                                                                     atttttCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCATTCTGCCCCTCCAAAAGCAAAAACCCGCCGaa                                                            >  2:810981/4‑149 (MQ=255)
                                                                                         ttCTGCCCCTCCAAAAGCAAACACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCgg                                                        <  2:1688150/149‑1 (MQ=37)
                                                                                         ttCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTTAAACTGAAAGATAAGCCGCTTTCTTTTTTGTATAGCGTCGTGGACAGTCAGTCAACGTTCTTCCACTCCAATAGAAAAAACCCGCCGAAGCgg                                                        >  2:3583150/1‑149 (MQ=11)
                                                                                         ttCGGCCCCTCCAAAAGGAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGACACTGACCGATAACGCGCGTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCgg                                                        <  2:1688266/149‑1 (MQ=11)
                                                                                          tCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCAATCATATTTCAGCCCCTCCAAAAGCAAAAAACCGCCGAAGCggg                                                       >  2:909031/1‑149 (MQ=255)
                                                                                                    cAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTa                                             >  2:3816766/1‑149 (MQ=37)
                                                                                                     aaaaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGTCAGTCATTCATCTTTGTGCCCCTCCCAAAGCAAAAAACCGCCGAAGCGGGGTTTTTCGTaa                                            >  1:2505536/1‑149 (MQ=11)
                                                                                                     aaaaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTaa                                            >  2:791026/1‑149 (MQ=37)
                                                                                                      aaaaCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACCGTCATTCATCTGTCTGGCCCGCCAAAATCCCAAACACGCCGAAGCGGGGGTTTAAGTaaa                                           >  1:3107701/1‑149 (MQ=11)
                                                                                                      aaaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTaaa                                           >  2:2516606/1‑149 (MQ=37)
                                                                                                      aaaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCACCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTaaa                                           >  2:752115/1‑149 (MQ=37)
                                                                                                      aaaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTATGCCCCTCCAAAATCAAAAACCCGCCGAAGCGGGTTTTTACGTaaa                                           >  1:3360145/1‑149 (MQ=255)
                                                                                                        aGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGGCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGAGGGTTTTTACGTAAATc                                         >  2:505340/1‑149 (MQ=255)
                                                                                                          cAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAg                                       <  1:434580/149‑1 (MQ=17)
                                                                                                             aaaCCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGa                                                                                                                                              <  1:3364519/43‑1 (MQ=38)
                                                                                                             aaaCCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGTCAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTg                                    >  1:670847/1‑149 (MQ=37)
                                                                                                             aaaCCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTg                                    >  1:2737408/1‑149 (MQ=37)
                                                                                                             aaaCCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCGTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAAACCGCCGAAGCGGGTTTTTACGTAAATCAGGTg                                    >  1:671659/1‑149 (MQ=255)
                                                                                                                   gccgAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCCAAAGCAAAAACCCGCCGAAGCGGGTGTTTACGTAAATCAGGTGAAACTg                              >  1:4105515/1‑149 (MQ=255)
                                                                                                                        aGCGGTTTATGACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGa                         <  2:3586228/149‑1 (MQ=255)
                                                                                                                        aGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGa                         <  1:2745221/149‑1 (MQ=37)
                                                                                                                           gggTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATaa                      <  1:3063113/149‑1 (MQ=37)
                                                                                                                             gTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGc                    >  2:1497956/1‑149 (MQ=37)
                                                                                                                             gTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGc                    >  2:3133888/1‑149 (MQ=37)
                                                                                                                             gTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCGCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGc                    >  2:1492668/1‑149 (MQ=37)
                                                                                                                                 ttACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCggg                <  1:1053394/149‑1 (MQ=38)
                                                                                                                                    cGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTc             <  1:918527/149‑1 (MQ=40)
                                                                                                                                    cGTAAACCAGGGGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTc             <  1:918537/149‑1 (MQ=255)
                                                                                                                                     gTAAACCACGTGAAACTGACCGATCAGCCTCTTTATTTTGGGTTTATCGTCGTGGAAAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCt            <  1:2838569/149‑1 (MQ=14)
                                                                                                                                       acatcAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCACAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGt          <  2:3516361/145‑1 (MQ=255)
                                                                                                                                          cccggcgAAACTGACCGATACGCCGCGTTCTTGTGGCTATAGCGTCGCGGTAAGGCATTCATGTGTCTGCCCCTCCACAACCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGt       <  1:2941476/143‑1 (MQ=14)
                                                                                                                                          ccAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGt       <  1:2943132/149‑1 (MQ=255)
                                                                                                                                          ccAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGt       <  1:884368/149‑1 (MQ=255)
                                                                                                                                           cAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAATCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTg      >  1:1627657/1‑149 (MQ=255)
                                                                                                                                            aGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATATTTGTGCCCCGCCAAAAACAAAAAACCCCCGAAGCGGGGTTTTAACTAAAAAACGTGAAAATGACCGATAAGCCGGGGTGTGtggtgg     >  1:3122619/1‑149 (MQ=11)
                                                                                                                                            aGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGGCAATAATATTTGTGCCCCACCAAAATCAAAAAACCCCCGAACCGGGTTTTTTCGTTAATAAGGTGAAACTTGCCGATAAACCCGGGGGTGTCGTgg     >  1:3122554/1‑149 (MQ=2)
                                                                                                                                            aGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCaa                                                                        <  1:4082848/82‑1 (MQ=18)
                                                                                                                                             ggTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATATTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGa    <  1:2972764/149‑1 (MQ=255)
                                                                                                                                             ggTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGGGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGa    <  1:2345035/149‑1 (MQ=255)
                                                                                                                                             ggTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGa    >  1:2290527/1‑149 (MQ=255)
                                                                                                                                             ggTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGa    >  1:2290408/1‑149 (MQ=255)
                                                                                                                                             ggTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGa    <  1:1019467/149‑1 (MQ=255)
                                                                                                                                              gTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGCCAGTCATTCATCTTTCTGCCCCTCCCAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTTAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGAc   >  2:1913784/1‑149 (MQ=255)
                                                                                                                                              gTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTATGCCCCTCCAAAAGCAAAAAACCGCCGAAGCGGGTTTTTAAGTAAATCAGGTGAAACAGACCGATAATCCGGGTTCTGTCGTGGAc   >  2:1237469/1‑149 (MQ=14)
                                                                                                                                               tGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACa  >  2:2030245/1‑149 (MQ=255)
                                                                                                                                               tGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACa  >  2:2028165/1‑149 (MQ=255)
                                                                                                                                                gAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGcgtcgt                                                                                                            <  2:3149718/42‑1 (MQ=38)
                                                                                                                                                |                                                                                                                                                   
CCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACA  >  CP009273/3263322‑3263613

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