Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,419,592 T→C 7.4% V537A (GTC→GCC)  PP_3026 → recombinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,419,5920TC 7.4% 842.9 / 36.7 271V537A (GTC→GCC) PP_3026recombinase
Reads supporting (aligned to +/- strand):  ref base T (105/146);  new base C (9/11);  total (114/157)
Fisher's exact test for biased strand distribution p-value = 8.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01

Only 100 of 276 total aligned reads displayed.

CGATACTTTCGATCGCATCGTGATTTACATGCGTGGTGTAGTGCCAGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGTCCCTCCCAGAAATTTTCCTGTACCTCGACGCATACACGCCAGGCGGCCCATGCGACATCTTTT  >  NC_002947/3419443‑3419741
                                                                                                                                                     |                                                                                                                                                     
cGATACTTTCGATCGCATCGTGATTTACATGCGTGGTGTAGTGCCAGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGt                                                                                                                                                       >  2:7323673/1‑150 (MQ=255)
 gATACTTTCGATCGCATCGTGATTTACATGCGTGGTGTAGTGCCAGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTc                                                                                                                                                      >  2:8473122/1‑150 (MQ=255)
 gATACTTTCGATCGCATCGTGATTTACATGCGTGGTGTAGTGCCAGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTc                                                                                                                                                      <  2:2281950/150‑1 (MQ=255)
 gATACTTTCGATCGCATCGTGATTTACATGCGTGGTGTAGTGCCAGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTc                                                                                                                                                      <  1:5692426/150‑1 (MQ=255)
     cTTTCGATCGCATCGTGATTTACATGCGTGGTGTAGTGCCAGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGAcc                                                                                                                                                  <  2:5670040/150‑1 (MQ=255)
         cGATCGCATCGTGATTTACATGCGTGGTGTAGTGCCAGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGTCGACGCCGGCCGccc                                                                                                                                              <  1:3862352/150‑1 (MQ=255)
          gATCGCATCGTGATTTACATGCGTGGTGTAGTGCCAGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGTCGACGCCGGCCGCCCa                                                                                                                                             <  2:2096683/150‑1 (MQ=255)
          gATCGCATCGTGATTTACATGCGTGGTGTAGTGCCAGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGTCGACGCCGGCCGCCCa                                                                                                                                             <  2:5824430/150‑1 (MQ=255)
          gATCGCATCGTGATTTACATGCGTGGTGTAGTGCCAGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCa                                                                                                                                             <  2:6346135/150‑1 (MQ=255)
          gATCGCATCGTGATTTACATGCGTGGTGTAGTGCCAGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCa                                                                                                                                             <  2:3837049/150‑1 (MQ=255)
           aTCGCATCGTGATTTACATGCGTGGTGTAGTGCCAGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCTAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAg                                                                                                                                            >  2:2930438/1‑150 (MQ=255)
           aTCGCATCGTGATTTACATGCGTGGTGTAGTGCCAGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAg                                                                                                                                            >  2:8778726/1‑150 (MQ=255)
              gCATCGTGATTTACATGCGTGGTGTAGTGCCAGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACt                                                                                                                                         >  1:6868624/1‑150 (MQ=255)
                aTCGTGATTTACATGCGTGGTGTAGTGCCAGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGa                                                                                                                                       <  2:1358641/150‑1 (MQ=255)
                    tGATTTACATGCGTGGTGTAGTGCCAGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGtt                                                                                                                                   >  1:86200/1‑150 (MQ=255)
                    tGATTTACATGCGTGGTGTAGTGCCAGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGtt                                                                                                                                   >  1:2993493/1‑150 (MQ=255)
                       tttACATGCGTGGTGTAGTGCCAGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTAc                                                                                                                                <  1:525685/150‑1 (MQ=255)
                        ttACATGCGTGGTGTAGTGCCAGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACg                                                                                                                               >  2:160226/1‑150 (MQ=255)
                         tACATGCGTGGTGTAGTGCCAGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGc                                                                                                                              <  2:5977384/150‑1 (MQ=255)
                            aTGCGTGGTGTAGTGCCAGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGa                                                                                                                           <  2:5913957/150‑1 (MQ=255)
                             tGCGTGGTGTAGTGCCAGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGAc                                                                                                                          <  2:4532052/150‑1 (MQ=255)
                              gCGTGGTGTAGTGCCAGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGAct                                                                                                                         >  1:7327166/1‑150 (MQ=255)
                                  ggtgTAGTGCCAGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGActtctt                                                                                                                     <  1:6529119/150‑1 (MQ=255)
                                      tAGTGCCAGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTgtcg                                                                                                                 >  2:2580952/1‑150 (MQ=255)
                                       aGTGCCAGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTgtcgt                                                                                                                <  2:450666/150‑1 (MQ=255)
                                       aGTGCCAGCGGGGCGCAGGGCCAGGTAAATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCTGCGATGGCTCAGGGTAGAAAGTAAGACCGGAGAGTTGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACTCTGACTTCTTgtcgt                                                                                                                <  2:374568/150‑1 (MQ=255)
                                        gTGCCAGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGa                                                                                                                                           <  1:3803170/122‑1 (MQ=255)
                                          gCCAGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTa                                                                                                             >  1:1783641/1‑150 (MQ=255)
                                            cAGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTAtt                                                                                                           >  1:3259548/1‑150 (MQ=255)
                                             aGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTg                                                                                                          >  1:1035323/1‑150 (MQ=255)
                                             aGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTg                                                                                                          <  1:4495616/150‑1 (MQ=255)
                                             aGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTg                                                                                                          <  1:2584020/150‑1 (MQ=255)
                                             aGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTg                                                                                                          <  2:8456871/150‑1 (MQ=255)
                                             aGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACAGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTg                                                                                                          >  1:5811007/1‑150 (MQ=255)
                                              gCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTgg                                                                                                         <  2:5202741/150‑1 (MQ=255)
                                                 gggCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAg                                                                                                      <  1:8178421/150‑1 (MQ=255)
                                                    cgcAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGcc                                                                                                   >  1:7830130/1‑150 (MQ=255)
                                                    cgcAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGcc                                                                                                   >  1:7128801/1‑150 (MQ=255)
                                                     gcAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCg                                                                                                  >  1:5702390/1‑150 (MQ=255)
                                                     gcAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCg                                                                                                  >  2:697226/1‑150 (MQ=255)
                                                      cAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCgg                                                                                                 <  2:5781497/150‑1 (MQ=255)
                                                      cAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCgg                                                                                                 <  2:6164630/150‑1 (MQ=255)
                                                      cAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCgg                                                                                                 <  1:1156047/150‑1 (MQ=255)
                                                        gggCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGt                                                                                               <  1:4600209/150‑1 (MQ=255)
                                                         ggCCAGGTACATCGATGTCTTGCTGATTTCCTAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTg                                                                                              <  1:2930438/150‑1 (MQ=255)
                                                          gCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTgg                                                                                             <  1:8346928/150‑1 (MQ=255)
                                                           ccAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTggg                                                                                            <  1:7769874/150‑1 (MQ=255)
                                                            cAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGc                                                                                           <  1:1379885/150‑1 (MQ=255)
                                                                tACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAgt                                                                                       >  2:4067494/1‑150 (MQ=255)
                                                                tACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAgt                                                                                       >  2:1831454/1‑150 (MQ=255)
                                                                    tCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGc                                                                                   <  2:8256966/150‑1 (MQ=255)
                                                                    tCGATGTCTTGCTGATTCCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGc                                                                                   <  2:1655721/150‑1 (MQ=255)
                                                                      gATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGcac                                                                                 <  2:1171219/150‑1 (MQ=255)
                                                                           cTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCAccacca                                                                            <  1:2435106/150‑1 (MQ=255)
                                                                               cTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTg                                                                        <  2:7514500/150‑1 (MQ=255)
                                                                               cTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTg                                                                        >  1:5357801/1‑150 (MQ=255)
                                                                                      ccAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTc                                                                                                                                                      <  2:7660949/65‑1 (MQ=255)
                                                                                      ccAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGt                                                                 <  1:2832394/150‑1 (MQ=255)
                                                                                       cAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGtc                                                                <  1:4752757/150‑1 (MQ=255)
                                                                                        aaGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGtcc                                                               >  1:4997089/1‑150 (MQ=255)
                                                                                        aaGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGtcc                                                               >  2:4834095/1‑150 (MQ=255)
                                                                                        aaGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGAACCACCAGGTGGCCTTGtcc                                                               >  2:6542825/1‑150 (MQ=255)
                                                                                          ggCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGtccct                                                             >  1:6052375/1‑150 (MQ=255)
                                                                                          ggCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGtccct                                                             >  1:4256417/1‑150 (MQ=255)
                                                                                           gCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGtccctc                                                            >  1:4392497/1‑150 (MQ=255)
                                                                                           gCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGtccctc                                                            >  1:4530333/1‑150 (MQ=255)
                                                                                            cAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGtccctcc                                                           <  1:7757925/150‑1 (MQ=255)
                                                                                             aGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGtccctccc                                                          >  1:8455533/1‑150 (MQ=255)
                                                                                             aGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGtccctccc                                                          >  2:3307388/1‑150 (MQ=255)
                                                                                             aGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGtccctccc                                                          >  2:8062397/1‑150 (MQ=255)
                                                                                                  aGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGTCCCTCCCAGaaa                                                     <  2:1568068/150‑1 (MQ=255)
                                                                                                     ggcgATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGTCCCTCCCAGAAAttt                                                  >  1:3105709/1‑150 (MQ=255)
                                                                                                      gcgATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGTCCCTCCCAGAAAtttt                                                 >  2:8576823/1‑150 (MQ=255)
                                                                                                      gcgATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGTCCCTCCCAGAAAtttt                                                 >  1:6286470/1‑150 (MQ=255)
                                                                                                       cgATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGTCCCTCCCAGAAATTTTc                                                <  2:5062087/150‑1 (MQ=255)
                                                                                                       cgATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGTCCCTCCCAGAAATTTTc                                                <  1:3321099/150‑1 (MQ=255)
                                                                                                       cgATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGTCCCTCCCAGAAATTTTc                                                <  2:8690911/150‑1 (MQ=255)
                                                                                                        gATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGTCCCTCCCAGAAATTTTcc                                               <  2:3160267/150‑1 (MQ=255)
                                                                                                        gATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGTCCCTCCCAGAAATTTTcc                                               <  1:7977978/150‑1 (MQ=255)
                                                                                                         aTGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGTCCCTCCCAGAAATTTTCCt                                              <  1:2365289/150‑1 (MQ=255)
                                                                                                            gCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGTCCCTCCCAGAAATTTTCCTGTa                                           >  1:4815510/1‑150 (MQ=255)
                                                                                                             cTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGTCCCTCCCAGAAATTTTCCTGTAc                                          <  1:6190571/150‑1 (MQ=255)
                                                                                                                        aaGGAAGACCGGAGAGTGGGCGGCTGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGTCCCTCCCAGAAATTTTCCTGTACCTCGACGCATa                               >  2:7250102/1‑150 (MQ=255)
                                                                                                                        aaGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGTCCCTCCCAGAAATTTTCCTGTACCTCGACGCATa                               <  1:2122140/150‑1 (MQ=255)
                                                                                                                        aaGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGTCCCTCCCAGAAATTTTCCTGTACCTCGACGCATa                               <  1:5449471/150‑1 (MQ=255)
                                                                                                                           gAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGTCCCTCCCAGAAATTTTCCTGTACCTCGACGCATacac                            <  1:1559052/150‑1 (MQ=255)
                                                                                                                             aGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGTCCCTCCCAGAAATTTTCCTGTACCTCGACGCATACACGc                          <  2:2447814/150‑1 (MQ=255)
                                                                                                                               aCCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGTCCCTCCCAGAAATTTTCCTGTACCTCGACGCATACACGCCa                        <  2:6466505/150‑1 (MQ=255)
                                                                                                                                 cGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGTCCCTCCCAGAAATTTTCCTGTACCTCGACGCATACACGCCAgg                      <  1:7113736/150‑1 (MQ=255)
                                                                                                                                 cGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGTCCCTCCCAGAAATTTTCCTGTACCTCGACGCATACACGCCAgg                      >  1:1528268/1‑150 (MQ=255)
                                                                                                                                  ggAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCAGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGTCCCTCCCAGAAATTTTCCTGTACCTCGACGCATACACGCCAggc                     >  1:5138447/1‑150 (MQ=255)
                                                                                                                                  ggAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCAGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGTCCCTCCCAGAAATTTTCCTGTACCTCGACGCATACACGCCAggc                     >  1:746596/1‑150 (MQ=255)
                                                                                                                                  ggAGAGTGGGCGGCCGGCATCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGTCCCTCCCAGAAATTTTCCTGTACCTCGACGCATACACGCCAggc                     <  2:4441511/150‑1 (MQ=255)
                                                                                                                                     gagTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGTCCCTCCCAGAAATTTTCCTGTACCTCGACGCATACACGCCAggcggc                  <  1:4023824/150‑1 (MQ=255)
                                                                                                                                         gggCGGTCGACGCCGGCCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGTCCCTCCCAGAAATTTTCCTGTACCTCGACGCATACACGCCAGGCGGCCCAt              >  1:972424/1‑150 (MQ=255)
                                                                                                                                          ggcggtCGACGCCGGCCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGTCCCTCCCAGAAATTTTCCTGTACCTCGACGCATACACGCCAGGCGGCCCATg             >  1:5470458/1‑150 (MQ=255)
                                                                                                                                          ggcggcCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGTCCCTCCCAGAAATTTTCCTGTACCTCGACGCATACACGCCAGGCGGCCCATg             <  2:2993494/150‑1 (MQ=255)
                                                                                                                                           gcggcCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGTCCCTCCCAGAAATTTTCCTGTACCTCGACGCATACACGCCAGGCGGCCCATGc            >  1:1337762/1‑150 (MQ=255)
                                                                                                                                           gcggcCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGTCCCTCCCAGAAATTTTCCTGTACCTCGACGCATACACGCCAGGCGGCCCATGc            <  2:4154098/150‑1 (MQ=255)
                                                                                                                                                     tCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGTCCCTCCCAGAAATTTTCCTGTACCTCGACGCATACACGCCAGGCGGCCCATGCGACATCtttt  <  2:6868624/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGATACTTTCGATCGCATCGTGATTTACATGCGTGGTGTAGTGCCAGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGTCCCTCCCAGAAATTTTCCTGTACCTCGACGCATACACGCCAGGCGGCCCATGCGACATCTTTT  >  NC_002947/3419443‑3419741

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