Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,455,018 T→A 5.1% T2T (ACA→ACT) ‡ PP_1268 ← HtrA suppressor protein SohA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,455,0180TA 5.1% 1002.0 / 27.0 294T2T (ACA→ACT) ‡PP_1268HtrA suppressor protein SohA
Reads supporting (aligned to +/- strand):  ref base T (135/144);  new base A (5/10);  total (140/154)
Fisher's exact test for biased strand distribution p-value = 2.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

Only 100 of 304 total aligned reads displayed.

TTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAGGTGGCGGTGGCCCATGATGATCAG  >  NC_002947/1454869‑1455166
                                                                                                                                                     |                                                                                                                                                    
tttttCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCtt                                                                                                                                                      <  1:8083922/150‑1 (MQ=255)
tttttCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCtt                                                                                                                                                      <  2:5952108/150‑1 (MQ=255)
 ttttCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTg                                                                                                                                                     >  2:6422519/1‑150 (MQ=255)
  tttCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCAATCGCATTACCTTga                                                                                                                                                    <  1:2113592/150‑2 (MQ=255)
   ttCAAGCCCTAGGCGCGCCGCTCCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTc                                                                                                                                                   <  2:5868335/150‑1 (MQ=255)
   ttCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTc                                                                                                                                                   >  2:4720154/1‑150 (MQ=255)
   ttCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGCAATGCAATCGCATTACCTTgag                                                                                                                                                   <  1:8524268/150‑3 (MQ=255)
      aaGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTTGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAAt                                                                                                                                                <  2:7455716/150‑1 (MQ=255)
      aaGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAAt                                                                                                                                                <  1:5115463/150‑1 (MQ=255)
       aGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATg                                                                                                                                               >  1:3680939/1‑150 (MQ=255)
           cAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAgg                                                                                                                                           >  1:8363927/1‑150 (MQ=255)
           cAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAgg                                                                                                                                           >  2:7682590/1‑150 (MQ=255)
            aaGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGTTACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGt                                                                                                                                          >  1:8415495/1‑150 (MQ=255)
            aaGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGt                                                                                                                                          >  1:6186962/1‑150 (MQ=255)
            aaGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGt                                                                                                                                          >  1:3982881/1‑150 (MQ=255)
            aaGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGACAAGGCAATGCaatc                                                                                                                                          <  2:8234114/150‑4 (MQ=255)
            aaGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGACAAGGCAATGCaatc                                                                                                                                          <  2:7741350/150‑4 (MQ=255)
            aaGGCGCGCCGCACCGGTTCGGGACCAGGGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGGCATGAGTAATTTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGACAAGGCAATGCaatc                                                                                                                                          <  2:8544825/150‑4 (MQ=255)
              ggCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGc                                                                                                                                        <  1:2381224/150‑1 (MQ=255)
                   gccgcACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTAc                                                                                                                                   >  1:7322035/1‑150 (MQ=255)
                    ccgcACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTAcc                                                                                                                                  >  2:2429670/1‑150 (MQ=255)
                    ccgcACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTAcc                                                                                                                                  <  2:1114308/150‑1 (MQ=255)
                     cgcACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAAGGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCa                                                                                                                                 <  1:904730/150‑1 (MQ=255)
                       cACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACg                                                                                                                               >  1:446784/1‑150 (MQ=255)
                        aCCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACgg                                                                                                                              <  2:279681/150‑1 (MQ=255)
                         ccGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGt                                                                                                                             <  1:2672196/150‑1 (MQ=255)
                            tttCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTTGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGtgat                                                                                                                          <  2:5988942/150‑1 (MQ=255)
                            tttCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGtgat                                                                                                                          <  2:6732428/150‑1 (MQ=255)
                             ttCGGGAACAGTTGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGtgatg                                                                                                                         <  2:898801/150‑1 (MQ=255)
                                gggAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGc                                                                                                                      <  1:3260314/150‑1 (MQ=255)
                                 ggAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCg                                                                                                                     >  2:1340916/1‑150 (MQ=255)
                                       gTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAgg                                                                                                               <  2:7540372/150‑1 (MQ=255)
                                            ttGATTGCGATCGGTCAGGGTGGACTCGACGTCGATGGTGGCAGCCATGAGTAATGTCGTGTGGGTGACAGCGATACATCTCTTGGTTATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCagag                                                                                                          <  2:8625171/150‑1 (MQ=255)
                                             tGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGc                                                                                                         <  2:5590594/150‑1 (MQ=255)
                                              gATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGcc                                                                                                        <  1:7793795/150‑1 (MQ=255)
                                               aTAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGccc                                                                                                       >  1:645222/1‑150 (MQ=255)
                                               aTAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGccc                                                                                                       >  2:196616/1‑150 (MQ=255)
                                                 aGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCccac                                                                                                     >  2:3975651/1‑150 (MQ=255)
                                                    gATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGc                                                                                                                      >  2:1827541/1‑130 (MQ=255)
                                                    gATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGc                                                                                                                      <  1:1827541/130‑1 (MQ=255)
                                                    gATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCaa                                                                                                  >  1:5865965/1‑150 (MQ=255)
                                                          tCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTgg                                                                                            >  1:6967370/1‑150 (MQ=255)
                                                          tCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTgg                                                                                            <  2:2293569/150‑1 (MQ=255)
                                                             gCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCAt                                                                                         <  1:357962/150‑1 (MQ=255)
                                                             gCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCAt                                                                                         <  2:3138972/150‑1 (MQ=255)
                                                                  gACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTa                                                                                    >  2:6756433/1‑150 (MQ=255)
                                                                     tCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAAt                                                                                                                                                <  1:2516978/87‑1 (MQ=255)
                                                                     tCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAAt                                                                                                                                                >  2:2516978/1‑87 (MQ=255)
                                                                     tCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTAccc                                                                                 >  1:7622303/1‑150 (MQ=255)
                                                                      cGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCg                                                                                >  1:5404069/1‑150 (MQ=255)
                                                                       gACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCgg                                                                               <  1:67003/150‑1 (MQ=255)
                                                                       gACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCgg                                                                               >  2:7056511/1‑150 (MQ=255)
                                                                       gACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCgg                                                                               >  1:301817/1‑150 (MQ=255)
                                                                       gACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCgg                                                                               <  1:5266189/150‑1 (MQ=255)
                                                                             gAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATc                                                                         >  1:100612/1‑150 (MQ=255)
                                                                             gAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATc                                                                         >  1:2850346/1‑150 (MQ=255)
                                                                                   ggCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTc                                                                   >  1:5945478/1‑150 (MQ=255)
                                                                                    gCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTcg                                                                  >  1:2067635/1‑150 (MQ=255)
                                                                                          aTGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGa                                                            >  1:4667167/1‑150 (MQ=255)
                                                                                            gAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACCCGGGCATGTCTACCCGGTCAATCACCTTCGCGCTGAtt                                                          <  2:4915757/150‑1 (MQ=255)
                                                                                            gAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGAtt                                                          <  2:5815050/150‑1 (MQ=255)
                                                                                                  tGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACggg                                                    >  2:5169935/1‑150 (MQ=255)
                                                                                                    tCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTc                                                  >  1:1250198/1‑150 (MQ=255)
                                                                                                      cTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCa                                                <  2:3446450/150‑1 (MQ=255)
                                                                                                       tgtgGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGTATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAg                                               <  1:8117749/150‑1 (MQ=255)
                                                                                                       tgtgGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAg                                               >  1:3787502/1‑150 (MQ=255)
                                                                                                       tgtgGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAg                                               >  2:5908646/1‑150 (MQ=255)
                                                                                                          gggTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCaa                                            >  2:6231534/1‑150 (MQ=255)
                                                                                                          gggTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACAGGTCCAGCaa                                            <  2:8415261/150‑1 (MQ=255)
                                                                                                          gggTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCATTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCTCCTTCGCGCTGATTGACGGGTCCAGCaa                                            <  2:2578063/150‑1 (MQ=255)
                                                                                                             tGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAAAGCAGGATCATTAACACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATCTCTACCCGGACAATCAACTTCGCGCTGATAGACGGGTCCAGCAACCg                                         >  2:2932575/1‑150 (MQ=255)
                                                                                                              gACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGc                                        <  1:6422778/150‑1 (MQ=255)
                                                                                                              gACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGc                                        <  2:1227674/150‑1 (MQ=255)
                                                                                                               aCAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGcc                                       >  1:5639950/1‑150 (MQ=255)
                                                                                                                 aGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGAGTGACGGGTCCAGCAACCGcccc                                     >  2:6932042/1‑150 (MQ=255)
                                                                                                                    gATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAAc                                  >  2:17920/1‑150 (MQ=255)
                                                                                                                    gATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGCCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAAc                                  >  2:423732/1‑150 (MQ=255)
                                                                                                                      tACAACTCAAGGTAATGCGATTGCATTGCCTTGTCTATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACcg                                <  1:7704325/150‑1 (MQ=255)
                                                                                                                      tACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACcg                                >  1:4233751/1‑150 (MQ=255)
                                                                                                                         aaCTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGa                             >  1:6930020/1‑150 (MQ=255)
                                                                                                                         aaCTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGa                             <  1:7116800/150‑1 (MQ=255)
                                                                                                                          aCTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGAc                            >  1:3828086/1‑150 (MQ=255)
                                                                                                                           cTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACa                           <  2:7322035/150‑1 (MQ=255)
                                                                                                                            tCAAGGTAATGCGATTGCATGGCCTTGTCAATGCAGGGGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGAGTGACGGGTCCAGCACCCGCCCCCACcgcgcccg                          >  2:8148943/1‑146 (MQ=255)
                                                                                                                                tgTAATGCGATTGCATTGCATTGTCAATGCAGGTGGATTACCACGTTGATGAGCGGTCAGGCAGAGCCCTCCAACACTGGCATGTCTTCCCGGTCAATCTCCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAggtgg                      <  2:5305088/149‑1 (MQ=255)
                                                                                                                                        gATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAGGTGGCGGTGGcc              >  2:6605415/1‑150 (MQ=255)
                                                                                                                                        gATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAGGTGGCGGTGGcc              >  1:244506/1‑150 (MQ=255)
                                                                                                                                         cctgcattgaCAAGGCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAGGTGGCGGTGGccc             >  1:1945027/3‑150 (MQ=255)
                                                                                                                                         cctgcattgaCAAGGCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAGGTGGCGGTGGccc             >  2:4049369/3‑150 (MQ=255)
                                                                                                                                         attgcattgcCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAGGTGGCGGTGGccc             >  2:15204/1‑150 (MQ=255)
                                                                                                                                          ttgcattgcCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAGGTGGCGGTGGCCCa            <  1:6301206/150‑1 (MQ=255)
                                                                                                                                             cattgcCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAGGTGGCGGTGGCCCatga         >  2:3692131/1‑150 (MQ=255)
                                                                                                                                             cattgcCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAGGTGGCGGTGGCCCatga         <  2:1881716/150‑1 (MQ=255)
                                                                                                                                              attgcCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAGGTGGCGGTGGCCCatgat        >  2:3803102/1‑150 (MQ=255)
                                                                                                                                                tgcCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAGGTGGCGGTGGCCCAtgatga      <  1:8627988/150‑1 (MQ=255)
                                                                                                                                                tgcCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAGGTGGCGGTGGCCCAtgatga      <  2:446784/150‑1 (MQ=255)
                                                                                                                                                    ttGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAGGTGGCGGTGGCCCATGATGATCAg  >  1:6248861/1‑150 (MQ=255)
                                                                                                                                                    ttGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAGGTGGCGGTGGCCCATGATGATCAg  >  1:684200/1‑150 (MQ=255)
                                                                                                                                                    ttGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAGGTGGCGGTGGCCCATGATGATCAg  >  2:8285404/1‑150 (MQ=255)
                                                                                                                                                    ttGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAGGTGGCGGTGGCCCATGATGATCAg  >  1:2344009/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAGGTGGCGGTGGCCCATGATGATCAG  >  NC_002947/1454869‑1455166

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