Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 898,208 C→A 5.9% P30Q (CCA→CAA)  PP_0781 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947898,2080CA 5.9% 213.4 / 2.1 68P30Q (CCA→CAA) PP_0781hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (35/29);  new base A (2/2);  total (37/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01

CTACAACTGACTCAGAAACCTCCGATCCGAATCCATGATGGGATTTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGAC  >  NC_002947/898061‑898356
                                                                                                                                                   |                                                                                                                                                    
cTACAACTGACTCAGAAACCTCCGATCCGAATCCATGATGGGATTTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCaa                                                                                                                                                    >  1:1317214/1‑150 (MQ=255)
  aCAACTGACTCAGAAACCTCCGATCCGAATCCATGATGGGATTTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGc                                                                                                                                                  <  2:475823/150‑1 (MQ=40)
    aaCTGACTCAGAAACCTCCGATCCGAATCCATGATGGGATTTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCgg                                                                                                                                                <  1:500496/150‑1 (MQ=255)
           tCAGAAACCTCCGATCCGAATCCATGATGGGATTTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAg                                                                                                                                         <  1:1183628/150‑1 (MQ=35)
           tCAGAAACCTCCGATCCGAATCCATGATGGGATTTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAg                                                                                                                                         >  1:853419/1‑150 (MQ=35)
             aGAAACCTCCGATCCGAATCCATGATGGGATTTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGtt                                                                                                                                       >  2:404050/1‑150 (MQ=35)
               aaaCCTCCGATCCGAATCCATGATGGGATTTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTcc                                                                                                                                     <  1:404050/150‑1 (MQ=35)
               aaaCCTCCGATCCGAATCCATGATGGGATTTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTcc                                                                                                                                     >  1:158704/1‑150 (MQ=35)
                    tccgatccgaATCCATGATGGGATTTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCg                                                                                                                                    <  1:893380/146‑1 (MQ=34)
                    tccgatccgaATCCATGATGGGATTTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCg                                                                                                                                    >  2:893380/1‑146 (MQ=35)
                    tccgatccgaATCCATGATGGGATTTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCAcc                                                                                                                                >  2:1245362/1‑150 (MQ=35)
                    tccgatccgaATCCATGATGGGATTTCCCCCTGAACACGTTGTGAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCAcc                                                                                                                                >  2:747235/1‑150 (MQ=25)
                     ccgatccgaATCCATGATGGGATTTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCa                                                                                                                               <  2:1153844/150‑1 (MQ=35)
                         tccgaATCCATGATGGGATTTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCTCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCt                                                                                                                           <  1:378381/150‑1 (MQ=25)
                               tCCATGATGGGTTGTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATga                                                                                                                     >  2:1096331/1‑150 (MQ=25)
                                   tgatgGGATTTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATgaagaa                                                                                                                 <  1:96655/150‑1 (MQ=37)
                                    gatgGGATTTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGaagaag                                                                                                                <  2:610925/150‑1 (MQ=37)
                                     atgGGATTTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGt                                                                                                               >  1:1413873/1‑150 (MQ=38)
                                      tgGGATTTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGCGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTCCCGCAGCAGCCTCAATGAAGAAGtg                                                                                                              >  1:346328/1‑149 (MQ=255)
                                        ggATTTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGc                                                                                                            >  1:547038/1‑150 (MQ=39)
                                         gATTTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCCACACCTCCTGGCACGAGATGGATTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGcc                                                                                                           >  1:610182/1‑150 (MQ=11)
                                          aTTTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCt                                                                                                          <  1:281210/150‑1 (MQ=39)
                                           tttCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTc                                                                                                         >  2:1132770/1‑150 (MQ=39)
                                                cccTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTgg                                                                                                    <  2:1317214/150‑1 (MQ=255)
                                                  cTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCGATCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGc                                                                                                  <  1:41522/150‑1 (MQ=255)
                                                        acGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAAc                                                                                            >  1:1199819/1‑150 (MQ=255)
                                                          gTTGACAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCGAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCt                                                                                          >  1:1439895/1‑150 (MQ=255)
                                                              tCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTgcgc                                                                                      <  2:547038/150‑1 (MQ=255)
                                                                 gggAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAAt                                                                                                                       <  2:1474351/114‑1 (MQ=38)
                                                                 gggAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAAt                                                                                                                       >  1:1474351/1‑114 (MQ=38)
                                                                    agaCACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATc                                                                                <  1:1245362/150‑1 (MQ=255)
                                                                        acCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATgaa                                                                                                                    >  2:953148/1‑110 (MQ=38)
                                                                        acCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATgaa                                                                                                                    <  1:953148/110‑1 (MQ=38)
                                                                           tCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCtg                                                                         >  1:566324/1‑150 (MQ=255)
                                                                               gTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCACCCTGCGCCGGATCGACGCTGTgcg                                                                     >  1:508732/1‑150 (MQ=255)
                                                                                        ttACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCtt                                                            >  1:883448/1‑150 (MQ=255)
                                                                                        ttACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCt                                                             >  2:72138/1‑149 (MQ=255)
                                                                                        ttACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCt                                                             <  1:72138/149‑1 (MQ=255)
                                                                                          atctcgtgccaggagggcacgaggtCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCt                                                          >  1:65338‑M2/26‑150 (MQ=255)
                                                                                          atctcgtgccaggagggcacgaggtCGAGTCACGTGTCATCGCCATGGTGGGCTACCAAAGc                                                                                                                                                  >  2:999761‑M2/26‑62 (MQ=255)
                                                                                          atctcgtgccaggagggcacgaggtCGAGTCACGTGTCATCGCCATGGTGGGCTACCAAAGc                                                                                                                                                  <  1:999761‑M2/37‑1 (MQ=255)
                                                                                          atctcgtgccaggagggcacgagatCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCt                                                          >  1:231039‑M2/26‑150 (MQ=255)
                                                                                          ACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCt                                                          >  1:747235‑M1/1‑150 (MQ=255)
                                                                                           tctcgtgccaggagggcacgaggtCGAGTCACGTGTCATCGCCATGGTGGGCTACCAAAGCGGTGGc                                                                                                                                            >  2:654438‑M2/25‑67 (MQ=255)
                                                                                           tctcgtgccaggagggcacgaggtCGAGTCACGTGTCATCGCCATGGTGGGCTACCAAAGCGGTGGc                                                                                                                                            <  1:654438‑M2/43‑1 (MQ=255)
                                                                                           tctcgtgccaggagggcacgagatTGAGTCACGTTTCATCGCCATGGTGAGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTa                                                         >  2:60652‑M2/25‑150 (MQ=255)
                                                                                            ctcgtgccaggagggcacgagatCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTgcgc                                                                                      >  1:592450‑M2/24‑120 (MQ=255)
                                                                                            ctcgtgccaggagggcacgagatCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTgcgc                                                                                      <  2:592450‑M2/97‑1 (MQ=255)
                                                                                            cTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTAc                                                        >  1:981438/1‑150 (MQ=255)
                                                                                            cTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTAc                                                        >  2:1495631/1‑150 (MQ=255)
                                                                                              cGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCa                                                      <  2:1199819/150‑1 (MQ=255)
                                                                                                 gccaggagggcacgagatCGAGTAACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGc                                                   <  2:231039‑M2/132‑1 (MQ=255)
                                                                                                 gCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGc                                                   <  2:1413873/150‑1 (MQ=255)
                                                                                                   cctcctGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCgg                                                                                                                                                >  1:1329273/1‑55 (MQ=38)
                                                                                                   cctcctGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCgg                                                                                                                                                <  2:1329273/55‑1 (MQ=38)
                                                                                                      cctGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAAt                                              <  2:981438/150‑1 (MQ=255)
                                                                                                          gCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATgaag                                                                                                                   >  1:254183/1‑77 (MQ=255)
                                                                                                          gCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATgaag                                                                                                                   <  2:254183/77‑1 (MQ=255)
                                                                                                          gCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGc                                                               <  2:167742/129‑1 (MQ=255)
                                                                                                          gCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGc                                                               >  1:167742/1‑129 (MQ=255)
                                                                                                            aCGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCTCTCAGCCTTCTACCAGGCTGAATGGCAGa                                        >  2:884738/1‑150 (MQ=255)
                                                                                                                   cGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGca                                 <  2:853419/150‑1 (MQ=255)
                                                                                                                         aCGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCt                           >  2:643734/1‑150 (MQ=255)
                                                                                                                             gtCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGc                       <  2:1439895/150‑1 (MQ=255)
                                                                                                                                 tCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAgc                   >  2:581846/1‑150 (MQ=255)
                                                                                                                                         gtggGAGACCCAAGCGGTGGCCAGGTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTCAATGGCAGACGGCGCACATGCTACGCCAGCGCTtcgtc           <  2:346328/150‑1 (MQ=255)
                                                                                                                                             gCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGAcc       <  1:643734/150‑1 (MQ=255)
                                                                                                                                              cTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCt      >  2:361707/1‑150 (MQ=255)
                                                                                                                                                  ccAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTtc              <  2:236714/138‑1 (MQ=255)
                                                                                                                                                  ccAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTtc              >  1:236714/1‑138 (MQ=255)
                                                                                                                                                  ccAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGAc  >  2:908863/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CTACAACTGACTCAGAAACCTCCGATCCGAATCCATGATGGGATTTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCTACCAGGCTGAATGGCAGACGGCGCACATGCTACGCCAGCGCTTCGTCGACCTCGAC  >  NC_002947/898061‑898356

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