Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,035,205 G→T 5.2% V13V (GTG→GTT PP_0893 → ThiJ/PfpI family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,035,2050GT 5.2% 301.6 / 3.3 92V13V (GTG→GTTPP_0893ThiJ/PfpI family protein
Reads supporting (aligned to +/- strand):  ref base G (44/43);  new base T (2/3);  total (46/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

TGAAGCAACGGATCTGGGGCTGGCTCGATCGGCTGGTGGTCAACTTCCTCGACCGCCGTGACTGACAGCGGCTATCGTGCAGTGACTATTCACAGTTAACCGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCTTGAGCATGGTGGGCCACACGGTGCATGCGGTGTGCCCGGAAAAGCTGGCGGGGCAGACCGTGCGTACCGCGATTCATGACTTCGAGGGGGAGCAGACCTACA  >  NC_002947/1035060‑1035353
                                                                                                                                                 |                                                                                                                                                     
tGAAGCAACGGATCTGGGGCTGGCTCGATCGGCTGGTGGTCAACTTCCTCGACCGCCGTGACTGACAGCGGCTATCGTGCAGTGACTATTCACAGTTAACCGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGCg                                                                                                                                                   >  1:1148877/1‑150 (MQ=255)
  aaGCAACGGATCTGGGGCTGGCTCGATCGGCTGGTGGTCAACTTCCTCGACCGCCGTGACTGACAGCGGCTATCGTGCAGTGACTATTCACAGTTAACCGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGCGAc                                                                                                                                                 <  2:467377/150‑1 (MQ=255)
       aCGGATCTGGGGCTGGCTCGATCGGCTGGTGGTCAACTTCCTCGACCGCCGTGACTGACAGCGGCTATCGTGCAGTGACTATTCACAGTTAACCGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGt                                                                                                                                            >  2:569811/1‑150 (MQ=255)
            tCTGGGGCTGGCTCGATCGGCTGGTGGTCAACTTCCTCGACCGCCGTGACTGACAGCGGCTATCGTGCAGTGACTATTCACAGTTAACCGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAg                                                                                                                                       <  2:1595331/150‑1 (MQ=255)
               ggggCTGGCTCGATCGGCTGGTGGTCAACTTCCTCGACCGCCGTGACTGACAGCGGCTATCGTGCAGTGACTATTCACAGTTAACCGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGAtt                                                                                                                                    <  1:569811/150‑1 (MQ=255)
                gggCTGGCTCGATCGGCTGGTGGTCAACTTCCTCGACCGCCGTGACTGACAGCGGCTATCGTGCAGTGACTATTCACAGTTAACCGAAGGAGTGGGCGAGATGACGGCGCAGAAAATTCTCATGCTGGTGGGCGACTATGTCTAAGATTa                                                                                                                                   <  2:919408/150‑1 (MQ=255)
                 ggctggctCGATCGGCTGGTGGTCAACTTCCTCGACCGCCGTGACTGACAGCGGCTATCGTGCAGTGACTATTCACAGTTAACCGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTAc                                                                                                                                  <  2:321586/150‑1 (MQ=255)
                  gctggctCGATCGGCTGGTGGTCAACTTCCTCGACCGCCGTGACTGACAGCGGCTATCGTGCAGTGACTATTCACAGTTAACCGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACg                                                                                                                                 >  1:746462/1‑150 (MQ=255)
                      gctCGATCGGCTGGTGGTCAACTTCCTCGACCGCCGTGACTGACAGCGGCTATCGTGCAGTGAATATTCACCATCCAGCGAAGGAGTGTGCGATATGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGc                                                                                                                             <  1:198342/150‑1 (MQ=255)
                       ctCGATCGGCTGGTGGTCAACTTCCTCGACCGCCGTGACTGACAGCGGCTATCGTGCAGTGACTATTCACAGTTAACCGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCg                                                                                                                            >  1:136193/1‑150 (MQ=255)
                          gATCGGCTGGTGGTCAACTTCCTCGACCGCCGTGACTGACAGCGGCTATCGTGCAGTGACTATTCACAGTTAACCGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCTGGTTGGCGACTACGTCGAAGATTACGAAGCGATg                                                                                                                         >  1:406307/1‑150 (MQ=255)
                            tcggacggtggtTAACTTCCTCGACCGCCGTGACTGACAGCGGCTATCGTGCAGTGACTATTCACAGTTAACCGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGt                                                                                                                       <  2:1325372/144‑1 (MQ=255)
                            tCGGCTGGTGGTCAACTTCCTCGACCGCCGTGACTGACAGCGGCTATCGTGCAGTGACGAGTCACAGTTAGCCGAAGGGGTGGGGGAGATGGCGGGGAAGAAAATTCGCATGCTGGTGGGCGACTACGTAGAAGAGTACAAAGCGATggg                                                                                                                       >  2:40240/1‑149 (MQ=255)
                                 tggtggTCAACTTCCTCGACCGCCGTGACTGACAGCGGCTATCGTGCAGTGACTATTCACAGTTAACCGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGc                                                                                                                                                    <  1:145636/116‑1 (MQ=255)
                                 tggtggTCAACTTCCTCGACCGCCGTGACTGACAGCGGCTATCGTGCAGTGACTATTCACAGTTAACCGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGc                                                                                                                                                    >  2:145636/1‑116 (MQ=255)
                                         aaCTTCCTCGACCGCCGTGACTGACAGCGGCTATCGTGCAGTGACTATTCACAGTTAACCGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGcc                                                                                                          >  2:1313150/1‑150 (MQ=255)
                                          aCTTCCTCGACCGCCGTGACTGACAGCGGCTATCGTGCAGTGACTATTCACAGTTAACCGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCt                                                                                                         <  2:136193/150‑1 (MQ=255)
                                           cTTCCTCGACCGCCGTGACTGACAGCGGCTATCGTGCAGTGACTATTCACAGTTAACCGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGcatc                                                                                                        >  2:780444/1‑147 (MQ=255)
                                           cTTCCTCGACCGCCGTGACTGACAGCGGCTATCGTGCAGTGACTATTCACAGTTAACCGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCtt                                                                                                        >  2:250592/1‑150 (MQ=255)
                                           cTTCCTCGACCGCCGTGACTGACAGCGGCTATCGTGCAGTGACTATTCACAGTTAACCGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCtt                                                                                                        <  1:50143/150‑1 (MQ=255)
                                             tCCTCGACCGCCGTGACTGACAGCGGCTATCGTGCAGTGACTATTCACAGTTAACCGAAGGAGTGGACGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCTTGa                                                                                                      <  1:510690/150‑1 (MQ=255)
                                                tCGACCGCCGTGACTGACAGCGGCTATCGTGCAGTGACTATTCACAGTTAACCGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCtggtgg                                                                                                                                                      >  1:464573/1‑99 (MQ=255)
                                                tCGACCGCCGTGACTGACAGCGGCTATCGTGCAGTGACTATTCACAGTTAACCGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCtggtgg                                                                                                                                                      <  2:464573/99‑1 (MQ=255)
                                                      gccgTGACTGACAGCGGCTATCGTGCAGTGACTATTCACAGTTAACCGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCTTGAGCAtggtgg                                                                                             <  2:306308/150‑1 (MQ=255)
                                                       ccgTGACTGACAGCGGCTATCGTGCAGTGACTATTCACAGTTAACCGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCTTGAGCATGGTggg                                                                                            <  2:746462/150‑1 (MQ=255)
                                                        cgTGACTGACAGCGGCTATCGTGCAGTGACTATTCACAGTTAACCGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCTTGAGCATGGTGGGc                                                                                           <  2:1148877/150‑1 (MQ=255)
                                                           gactgacAGCGGCTATCGTGCAGTGACTATTCACAGTTAACCGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCTTGAGCATGGTGGGCcac                                                                                        >  1:217812/1‑150 (MQ=255)
                                                            actgacAGCGGCTATCGTGCAGTGACTATTCACAGTTAACCGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGaa                                                                                                                               >  2:1233880/1‑110 (MQ=255)
                                                            actgacAGCGGCTATCGTGCAGTGACTATTCACAGTTAACCGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGaa                                                                                                                               <  1:1233880/110‑1 (MQ=255)
                                                               gacAGCGGCTATCGTGCAGTGACTATTCACAGTTAACCGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCTTGAGCATGGTGGGCCACACgg                                                                                    >  2:1086094/1‑150 (MQ=255)
                                                                  aGCGGCTATCGTGCAGTGACTATTCACAGTTAACCGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCTTGAGCAtggtg                                                                                              >  1:669816/1‑137 (MQ=255)
                                                                  aGCGGCTATCGTGCAGTGACTATTCACAGTTAACCGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCTTGAGCAtggtg                                                                                              <  2:669816/137‑1 (MQ=255)
                                                                      gCTATCGTGCAGTGACTATTCACAGTTAACCGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCTTGAGCATGGTGGGCCACACGGTGCATGc                                                                             <  1:1509077/150‑1 (MQ=255)
                                                                        tATCGTGCAGTGACTATTCACAGTTAACCGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCTTGAGCATGGTGGGCCACACGGTGCATGCgg                                                                           >  1:755769/1‑150 (MQ=255)
                                                                         aTCGTGCAGTGACTATTCACAGTTAACCGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCTTGAGCATGGTGGGCCACACGGTGCATGCGgt                                                                          <  2:527596/150‑1 (MQ=255)
                                                                            gTGCAGTGACTATTCACAGTTAACCGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCtggtgg                                                                                                                                                      <  2:369032/71‑1 (MQ=255)
                                                                            gTGCAGTGACTATTCACAGTTAACCGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCtggtgg                                                                                                                                                      >  1:369032/1‑71 (MQ=255)
                                                                            gTGCAGTGACTATTCACAGTTAACCGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCTTGAg                                                                                                     <  1:1396247/120‑1 (MQ=255)
                                                                            gTGCAGTGACTATTCACAGTTAACCGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCTTGAg                                                                                                     >  2:1396247/1‑120 (MQ=255)
                                                                              gCAGTGACTATTCACAGTTAACCGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCTTGAGCATGGTGGGCCACACGGTGCATGCGGTGTGcc                                                                     <  2:217812/150‑1 (MQ=255)
                                                                               cAGTGACTATTCACAGTTAACCGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCTGGTTGGCGACTACGTCGAAGATTACGAAg                                                                                                                              <  1:221439/92‑1 (MQ=255)
                                                                               cAGTGACTATTCACAGTTAACCGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAg                                                                                                                              >  2:221439/1‑92 (MQ=255)
                                                                                  tGACTATTCACAGTTAACCGAAGGAGAGGGCGAGCTGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCTTGAGCATGGTGGGCCACACGGTGCATGCGGTGTGCCCGGa                                                                 >  2:1253102/1‑150 (MQ=255)
                                                                                      tATTCACAGTTAACCGAAGGAGTGGGCGAGATGACGGCGACGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCTTGAGCATGGTGGGCCACACGGTGCATGCGGTGTGCCCGGAAAAg                                                             >  2:189326/1‑150 (MQ=255)
                                                                                        ttCACAGTTAACCGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCTTGAGCATGGTGGGCCACACGGTGCATGCg                                                                            <  2:265896/133‑1 (MQ=255)
                                                                                        ttCACAGTTAACCGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCTTGAGCATGGTGGGCCACACGGTGCATGCg                                                                            >  1:265896/1‑133 (MQ=255)
                                                                                         tCACAGTTAACCGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCTTGAGCATGGTGGGCCACACGGTGCATGCGGTGTGCCCGGAAAAGCTg                                                          >  1:498177/1‑150 (MQ=255)
                                                                                           acaGTTAACCGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCTGGTTGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCTTGAGCATGGTGGGCCACACGGTGCATGCGGTGTGCCCGGAAAAGCTGGc                                                        >  1:330412/1‑150 (MQ=255)
                                                                                           acaGTTAACCGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCTTGAGCATGGTGGGCCACACGGTGCATGCGGTGTGCCCGGAAAAGCTGGc                                                        >  2:737288/1‑150 (MQ=255)
                                                                                           acaGTTAACCGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCTTGAGCATGGTGGGCCACACGGTGCATGCGGTGTGCCCGGAAAAGCTGGc                                                        >  2:44049/1‑150 (MQ=255)
                                                                                              gTTAACCGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCTTg                                                                                                       >  1:184708/1‑100 (MQ=255)
                                                                                              gTTAACCGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCTTg                                                                                                       <  2:184708/100‑1 (MQ=255)
                                                                                              gTTAACCGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCTTGAGCATGGTGGGCCACACGGTGCATGCGGTGTGCCCGGAAAAGCTGGCggg                                                     >  2:548871/1‑150 (MQ=255)
                                                                                               ttAACCGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGCCTATGCAGATCCAGGCCTTGAGCATGGTGGGCCAAACGCTGCATGCGGTGTGCCCGTAAAAGCTGGCgggg                                                    >  1:1264105/1‑150 (MQ=255)
                                                                                               ttAACCGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCTTGAGCATGGTGGGCCACACGGTGCATGCGGTGTGCCCGGAAAAGCTGGCgggg                                                    <  1:44049/150‑1 (MQ=255)
                                                                                               ttAACCGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCTTGAGCATGGTGGGCCACACGGTGCATGCGGTGTGCCCGGAAAAGCTGGCgggg                                                    <  2:82369/150‑1 (MQ=255)
                                                                                                tAACCGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCTTGAGCATGGTGGGCCACACGGTGCATGCGGTGTGCCCGGAAAAGCTGGCGGGGc                                                   <  1:1086094/150‑1 (MQ=255)
                                                                                                tAACCGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCTTGAGCATGGTGGGCCACACGGTGCATGCGGTGTGCCCGGAAAAGCTGGCGGGGc                                                   <  1:250592/150‑1 (MQ=255)
                                                                                                 aaCCGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGcc                                                                                                                    >  1:85986/1‑84 (MQ=255)
                                                                                                 aaCCGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGcc                                                                                                                    <  2:85986/84‑1 (MQ=255)
                                                                                                  aCCGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCt                                                                                                         <  2:422262/94‑1 (MQ=255)
                                                                                                  aCCGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCt                                                                                                         >  1:422262/1‑94 (MQ=255)
                                                                                                   ccGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCTTGAGCATGGTGGGCCACACGGTGCATGCGGTGTGCCCGGAAAAGCTGGCGGGGCAGa                                                <  1:1313150/150‑1 (MQ=255)
                                                                                                    cGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAg                                                                                                             <  1:123967/88‑1 (MQ=255)
                                                                                                    cGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAg                                                                                                             >  2:123967/1‑88 (MQ=255)
                                                                                                    cGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGAGATAATGGCGTTTCAGGCCTTGAGCATGGTGGGCCACACGGTGCATGCGGTGTGCCCGGAAAAGCTGGCGGGGCAGAc                                               >  1:1017457/1‑150 (MQ=255)
                                                                                                     gAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCTTGAGCATGGTGGGCCACACGGTGCATGCGGTGTGCCCGGAAAAGCTGGCGGGGCAGAcc                                              >  1:116932/1‑150 (MQ=255)
                                                                                                     gAAGGAGGGGGCGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCTTGAGCATGGTGGGCCACACGGTGCATGCGGTGTGCCCGGAAAAGCTGGCGGGGCAGAcc                                              >  2:455153/1‑150 (MQ=255)
                                                                                                      aaGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCTTGAGCATGGTGGGCCACACGGTGCATGCGGTGTGCCCGGAAAAGCTGGCGGGGCAGACCg                                             >  2:852632/1‑150 (MQ=255)
                                                                                                              ggCGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCTTGAGCATGGTGGGCCACACGGTGCATGCGGTGTGCCCGGAAAAGCTGGCGGGGCAGACCGTGCGTAcc                                     <  1:1620909/150‑1 (MQ=255)
                                                                                                               gCGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCTTGAGCATGGTGGGCCACACGGTGCATg                                                                              >  2:1356394/1‑108 (MQ=255)
                                                                                                               gCGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCTTGAGCATGGTGGGCCACACGGTGCATg                                                                              <  1:1356394/108‑1 (MQ=255)
                                                                                                                  agaTGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCTTGAGCATGGTGGGCCACACGGTGCATGCGGTGTGCCCGGAAAAGCTGGCGGGGCa                                                  <  2:1539933/133‑1 (MQ=255)
                                                                                                                  agaTGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCTTGAGCATGGTGGGCCACACGGTGCATGCGGTGTGCCCGGAAAAGCTGGCGGGGCa                                                  >  1:1539933/1‑133 (MQ=255)
                                                                                                                  agaTGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCTTGAGCATGGTGGGCCACACGGTGCATGCGGTGTGCCCGGAAAAGCTGGCGGGGCAGACCGTGCGTACCGCGa                                 >  2:983263/1‑150 (MQ=255)
                                                                                                                       aCGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCTTGAGCATGGTGGGCCACACGGTGCATGCGGTGTGCCCGGAAAAGCTGGCGGGGCAGACCGTGCGTACCGCGATTCAt                            <  1:1289317/150‑1 (MQ=255)
                                                                                                                       aCGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCTTGAGCATGGTGGGCCACACGGTGCATGCGGTGTGCCCGGAAAAGCTGGCGGGGCAGACCGTGCGTACCGCGATTCAt                            <  2:498177/150‑1 (MQ=255)
                                                                                                                         ggCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCTTGAGCATGGTGGGCCACACGGTGCATGCGGTGTGCCCGGAAAAGCTGGCgggg                                                    <  2:266912/124‑1 (MQ=255)
                                                                                                                         ggCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCTTGAGCATGGTGGGCCACACGGTGCATGCGGTGTGCCCGGAAAAGCTGGCgggg                                                    >  1:266912/1‑124 (MQ=255)
                                                                                                                         ggCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCTTGAGCATGGTGGGCCACACGGTGCATGCGGTGTGCCCGGAAAAGCTGGCGGGGCAGACCGTGCGTACCGCGATTCATGa                          >  1:1188231/1‑150 (MQ=255)
                                                                                                                             aagaaAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCTTGAGCATGGTGGGCcaca                                                                                       >  1:565915/1‑85 (MQ=255)
                                                                                                                             aagaaAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCTTGAGCATGGTGGGCcaca                                                                                       <  2:565915/85‑1 (MQ=255)
                                                                                                                             aagaaAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCTTGAGCATGGTGGGCCACACGGTGCATGCGGTGTGCCCGGAAAAGCTGGCggg                                                     <  1:276064/119‑1 (MQ=255)
                                                                                                                             aagaaAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCTTGAGCATGGTGGGCCACACGGTGCATGCGGTGTGCCCGGAAAAGCTGGCggg                                                     >  2:276064/1‑119 (MQ=255)
                                                                                                                              agaaAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCTTGAGCATGGTGGGCCACACGGTGCATGCGGTGTGCCCGGAAAAGCTGGCGGGGCAGACCGTGCGTACCGCGATTCATGACTTCg                    >  1:1262993/1‑150 (MQ=255)
                                                                                                                                aaaaTTCTCATGCTGGTTGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCTTGAGCATGGTGGGCCACACGGTGCATGCGGTGTGCCCGGAAAAGCTGGCGGGGCAGACCGTGCGTACCGCGATTCATGACTTCGAg                  <  2:330412/150‑1 (MQ=255)
                                                                                                                                  aaTTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCTGTTCCAGGCCtt                                                                                                        <  2:25284/63‑1 (MQ=255)
                                                                                                                                  aaTTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCTGTTCCAGGCCtt                                                                                                        >  1:25284/1‑63 (MQ=255)
                                                                                                                                  aaTTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCTTGAGCATGGTGGGCCACACGGTGCATGCGGTGTGCCCGGAAAAGCTGGCGGGGCAGACCGTGCGTACCGCGATTCATGACTTCGAggg                >  1:712632/1‑150 (MQ=255)
                                                                                                                                   aTTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCTTGAGCATGGTGGGCCACACGGTGCATGCGGTGTGCCCGGAAAAGCTGGCGGGGCAGACCGTGCGTACCGCGATTCATGACTTCGAgggg               <  2:1262993/150‑1 (MQ=255)
                                                                                                                                     tctcATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCTTGAGCATGGTGGGCCACACGGTGCATGCGGTGTGCCCGGAAAAGCTGGCGGGGCAGACCGTGCGTACCGCGa                                 >  2:1228523/1‑131 (MQ=255)
                                                                                                                                     tctcATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCTTGAGCATGGTGGGCCACACGGTGCATGCGGTGTGCCCGGAAAAGCTGGCGGGGCAGACCGTGCGTACCGCGa                                 <  1:1228523/131‑1 (MQ=255)
                                                                                                                                     tctcATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCTTGAGCATGGTGGGCCACACGGTGCATGCGGTGTGCCCGGAAAAGCTGGCGGGGCAGACCGTGCGTACCGCGATTCATGACTTCGAGGGGGa             <  1:40240/150‑1 (MQ=255)
                                                                                                                                     tctcATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCTTGAGCATGGTGGGCCACACGGTGCATGCGGTGTGCCCGGAAAAGCTGGCGGGGCAGACCGTGCGTACCGCGATTCATGACTTCGAGGGGGa             <  2:712632/150‑1 (MQ=255)
                                                                                                                                     tctcATGCTGGTGGGCGACTACGTCGAAGATTAAGAAGCGATGGTGCCGTTCCAGGCCTTGAGCATGGTGGGCCACACGGTGCATGCGGTGTGCCCGGAAAAGCTGGCGGGGCAGACCGTGCGTACCGCGATTCATGACTTCGAGGGGGa             >  1:1621493/1‑150 (MQ=255)
                                                                                                                                        cATGCTGGTTGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCTTGAGCATGGTGGGCCACACGGTGCATGCGGTGTGCCCGGAAAAGCTGGCGGGGCAGACCGTGCGTACCGCGATTCATGGCATTCGggggggtgg           <  2:406307/150‑4 (MQ=255)
                                                                                                                                                tggGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCTTGAGCATGGTGGGCCACACGGTGCATGCGGTGTGCCCGGAAAAGCTGGCGGGGCAGACCGTGCGTACCGCGATTCATGACTTCGAGGGGGAGCAGACCTACa  <  2:1264105/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
TGAAGCAACGGATCTGGGGCTGGCTCGATCGGCTGGTGGTCAACTTCCTCGACCGCCGTGACTGACAGCGGCTATCGTGCAGTGACTATTCACAGTTAACCGAAGGAGTGGGCGAGATGACGGCGAAGAAAATTCTCATGCTGGTGGGCGACTACGTCGAAGATTACGAAGCGATGGTGCCGTTCCAGGCCTTGAGCATGGTGGGCCACACGGTGCATGCGGTGTGCCCGGAAAAGCTGGCGGGGCAGACCGTGCGTACCGCGATTCATGACTTCGAGGGGGAGCAGACCTACA  >  NC_002947/1035060‑1035353

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: