Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 401,316:1 +C 5.1% coding (899/1059 nt) PP_0334 → transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947401,3161.C 5.1% 261.7 / 10.7 78coding (899/1059 nt)PP_0334transposase
Reads supporting (aligned to +/- strand):  ref base . (37/37);  new base C (2/2);  total (39/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.43e-01

ACAGCAGCGGAAACAAAATCAGATTGGGCAAGATGAGCAAGCGCGGTGACGGCTACATACGCAGCCTGATGATCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCT  >  NC_002947/401190‑401464
                                                                                                                               |                                                                                                                                                    
gccccggtggcattacccggttggtcagaggctcatagccgcctccggtttttcgttcggtttgcgttggtaaacctgatcgctctgtaacagggcccatataatccgcagatttCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGc                                                                                                                               <  1:2441401‑M1/35‑1 (MQ=255)
                        tcagaggctcatagccgcctccggtttttcgttcggtttgcgttggtaaacctgatcgctctgtaacagggcccatataatccgcagattcCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAaggagg                                                                                                       <  1:2465046‑M1/59‑1 (MQ=255)
                        tcagaggctcatagccgcctccggtttttcgttcggtttgcgttggtaaacctgatcgctctgtaacagggcccagataatccgcagattcCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAaggagg                                                                                                       >  1:2613510‑M1/92‑150 (MQ=255)
                            aggctcatagccgcctccggtttttcgttcggtttgcgttggtaaacctgatcgctctgtaacagggcccatataatccgcagattcCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAggcggc                                                                                                   >  2:2465046‑M1/88‑150 (MQ=255)
                                           tccggtttttcgttcggtttgcgttggtaaacctgatcgctctgtaacagggcccagataatccgcagattcCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAaggag                                                                                                        <  1:282920‑M1/58‑1 (MQ=255)
                                           tccggtttttcgttcggtttgcgttggtaaacctgatcgctctgtaacagggcccagataatccgcagattcCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAaggag                                                                                                        >  2:282920‑M1/73‑130 (MQ=255)
                                            ccggtttttcgttcggtttgcgttggtaaacctgatcgctctgtaacagggcccagataatccgcagatttCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAAt                                                                                   >  1:894227‑M1/72‑150 (MQ=255)
                                                  tttcgttcggtttgcgttggtaaacctgatcgctctgtaacagggcccatataatccgcagatttCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAAt                                                                             >  2:2441401‑M1/66‑150 (MQ=255)
                                                               gCCTGATGATCCAAGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCCCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCt                                                                <  1:1121499/150‑1 (MQ=14)
                                                                ccTGATGATCCCGGGCGCGCCTTAAGTTCTGAGACATTTTTGCCCGGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTg                                                               <  2:2335884/150‑1 (MQ=11)
                                                                 cTGATGATCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTgg                                                              <  1:151622/150‑1 (MQ=32)
                                                                 cTGATGATCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTgg                                                              >  1:1822071/1‑150 (MQ=32)
                                                                  ttggtaaacctgatcgctctgtaacagggcccagataatccgcagatttCGATTAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAggcg                                                                                                     <  1:677978‑M1/61‑1 (MQ=255)
                                                                  ttggtaaacctgatcgctctgtaacagggcccagataatccgcagatttCGATTAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAggcg                                                                                                     >  2:677978‑M1/50‑110 (MQ=255)
                                                                  tgatgaTCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTggg                                                             >  2:1964482/1‑150 (MQ=32)
                                                                   gatgatCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGc                                                            >  2:2904624/1‑150 (MQ=32)
                                                                    atgatCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGATCGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGcc                                                           <  1:1463490/150‑1 (MQ=17)
                                                                      taaacctgatcgctctgtaacagggcccatataatccgcagatttCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAggcg                                                                                                     <  1:770337‑M1/61‑1 (MQ=255)
                                                                      taaacctgatcgctctgtaacagggcccatataatccgcagatttCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAggcg                                                                                                     >  2:770337‑M1/46‑106 (MQ=255)
                                                                        tCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGt                                                       >  2:2320780/1‑150 (MQ=32)
                                                                        tCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGt                                                       <  1:351248/150‑1 (MQ=32)
                                                                        tCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGt                                                       >  1:250916/1‑150 (MQ=32)
                                                                        tCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGt                                                       >  2:1064077/1‑150 (MQ=32)
                                                                        tCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGt                                                       >  2:1479001/1‑150 (MQ=32)
                                                                        tCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGt                                                       >  2:1096649/1‑150 (MQ=32)
                                                                        aacctgatcgctctgtgacagggcccagataatccgcagatttCGATTAGCCAGCCTGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAaggagg                                                                                                       <  2:1731921‑M1/60‑1 (MQ=255)
                                                                        aacctgatcgctctgtgacagggcccagataatccgcagatttCGATTAGCCAGCCTGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAaggagg                                                                                                       >  1:1731921‑M1/44‑103 (MQ=255)
                                                                          cAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTa                                                     <  2:1027335/150‑1 (MQ=32)
                                                                            gggCGCGCATTCCGTGTTGAGCCAGCTAGGGCCTGATTCCGATCAGGGAGG‑GGAGCGCCGCCTCTTGCGCTGGCTAGCCCGCCTTGGGCGCAAGGAGGCGGCGGTCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACa                                                   <  2:1352522/150‑1 (MQ=11)
                                                                            gggCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTCGCCTGCCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACa                                                   <  1:1725100/150‑1 (MQ=255)
                                                                                  gcATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGa                                             <  1:1064077/150‑1 (MQ=32)
                                                                                    aTTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATc                                           >  1:1564217/1‑150 (MQ=32)
                                                                                    aTTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATc                                           >  2:240981/1‑150 (MQ=32)
                                                                                     ttCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCa                                          <  1:738898/150‑1 (MQ=32)
                                                                                     ttCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCa                                          <  1:1964482/150‑1 (MQ=32)
                                                                                       ccGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAgg                                        >  1:409072/1‑150 (MQ=32)
                                                                                       ccGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCAAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAgg                                        <  1:941770/150‑1 (MQ=17)
                                                                                        cGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGc                                               <  1:2804375/142‑1 (MQ=32)
                                                                                        cGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGc                                               >  2:2804375/1‑142 (MQ=32)
                                                                                        cGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGt                                       >  1:2035775/1‑150 (MQ=32)
                                                                                          tCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGttt                                     <  2:1822071/150‑1 (MQ=32)
                                                                                          tCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGttt                                     <  2:1603054/150‑1 (MQ=32)
                                                                                          cagggcccatataatccgcagatttCGATTAGCCAGCCTGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAggc                                                                                                      <  2:500986‑M1/61‑1 (MQ=255)
                                                                                          cagggcccatataatccgcagatttCGATTAGCCAGCCTGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAggc                                                                                                      >  1:500986‑M1/26‑86 (MQ=255)
                                                                                           cTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTa                                    <  2:250916/150‑1 (MQ=32)
                                                                                            ttGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTAc                                   <  1:1479001/150‑1 (MQ=32)
                                                                                            ttGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TCATCGCCGCCCCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTAc                                   <  1:2465882/150‑1 (MQ=17)
                                                                                              gAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGACCAGGTTCACCa                                 >  2:794241/1‑150 (MQ=17)
                                                                                               aGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCaa                                <  2:191150/150‑1 (MQ=32)
                                                                                                gCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAAc                               <  1:1061729/150‑1 (MQ=32)
                                                                                                 ccAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACg                              >  1:1042523/1‑150 (MQ=32)
                                                                                                       aaGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAAcc                        >  2:2792061/1‑150 (MQ=32)
                                                                                                        aGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTAc                                   >  2:1085631/1‑138 (MQ=33)
                                                                                                        aGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTAc                                   <  1:1085631/138‑1 (MQ=33)
                                                                                                          gCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACcgaa                     >  2:2369777/1‑150 (MQ=32)
                                                                                                           ccTGATTCCGATCAGCCAGA‑TGATCTCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGATTTACCAACGCAAACcgaac                    >  2:1498951/1‑150 (MQ=17)
                                                                                                               aTTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACcgaac                    >  2:363328/1‑146 (MQ=32)
                                                                                                               aTTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACcgaac                    <  1:363328/146‑1 (MQ=32)
                                                                                                               aTTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGaaa                >  1:989547/1‑150 (MQ=32)
                                                                                                                   CGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAAcc            <  1:2911287‑M1/150‑1 (MQ=255)
                                                                                                                    GATCAGCCAGA‑TGATCGCCGGCTCTTGACCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCg           <  2:372052‑M1/150‑1 (MQ=255)
                                                                                                                    GATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGACCGAAAAACCg           <  1:1667402‑M1/150‑1 (MQ=255)
                                                                                                                    GATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCg           >  1:2828835‑M1/1‑150 (MQ=255)
                                                                                                                     aTCAGCCAGA‑TGACCGCAGGAAATTGCACTGGCTAGCCCGCCTTGAAAGCAACGCGGCGAGTAACGGAAAGTCTGATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCgg          <  1:55965/150‑1 (MQ=2)
                                                                                                                     ATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAAc                         <  1:69297‑M1/135‑1 (MQ=255)
                                                                                                                     ATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAAc                         >  2:69297‑M1/1‑135 (MQ=255)
                                                                                                                     ATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCgg          >  2:2730967‑M1/1‑150 (MQ=255)
                                                                                                                      TCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGa         <  2:2066693‑M1/150‑1 (MQ=255)
                                                                                                                      TCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGa         >  1:1075983‑M1/1‑150 (MQ=255)
                                                                                                                        aGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAgg       <  2:2911287/150‑1 (MQ=32)
                                                                                                                         gCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAggc      >  1:2341991/1‑150 (MQ=32)
                                                                                                                         gCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAggc      >  2:2668327/1‑150 (MQ=32)
                                                                                                                           cAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAggcgg    <  1:2792061/150‑1 (MQ=32)
                                                                                                                            aGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTgg                                                              <  2:1968025/91‑1 (MQ=33)
                                                                                                                            aGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTgg                                                              >  1:1968025/1‑91 (MQ=34)
                                                                                                                            aGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAAc                         >  1:459139/1‑128 (MQ=33)
                                                                                                                            aGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAAc                         <  2:459139/128‑1 (MQ=33)
                                                                                                                            aGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAggcggc   >  2:614556/1‑150 (MQ=32)
                                                                                                                             ga‑tgatCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAggcgg    <  1:725444/148‑1 (MQ=32)
                                                                                                                             ga‑tgatCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAggcgg    >  2:725444/1‑148 (MQ=32)
                                                                                                                             ga‑tgatCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCt  >  1:853877/1‑150 (MQ=32)
                                                                                                                               |                                                                                                                                                    
ACAGCAGCGGAAACAAAATCAGATTGGGCAAGATGAGCAAGCGCGGTGACGGCTACATACGCAGCCTGATGATCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCT  >  NC_002947/401190‑401464

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: ‑