Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,272,541 G→T 6.8% C12F (TGT→TTT) ‡ chaC → glutathione‑specific gamma‑glutamylcyclotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,272,5410GT 6.8% 364.1 / 10.6 109C12F (TGT→TTT) ‡chaCglutathione‑specific gamma‑glutamylcyclotransferase
Reads supporting (aligned to +/- strand):  ref base G (51/50);  new base T (5/3);  total (56/53)
Fisher's exact test for biased strand distribution p-value = 7.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

Only 100 of 111 total aligned reads displayed.

AAGTGGTCCGTGGATTGCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATCCGGCACTGG  >  NC_000913/1272393‑1272690
                                                                                                                                                    |                                                                                                                                                     
aaGTGGTCCGTGGATTGCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGt                                                                                                                                                      <  1:369566/150‑1 (MQ=255)
  gTGGTCCGTGGATTGCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTaa                                                                                                                                                    >  2:480686/1‑150 (MQ=255)
            gATTGCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCAttt                                                                                                                                          >  1:15693/1‑150 (MQ=255)
             aTTGCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGTTTTACAATCGGCAttca                                                                                                                                         <  2:329286/150‑3 (MQ=255)
             aTTGCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTg                                                                                                                                         >  1:997056/1‑150 (MQ=255)
               tGCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGt                                                                                                                                       >  1:1149637/1‑150 (MQ=255)
               tGCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATCTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGTTTTACAATCGGCAttcatc                                                                                                                                       <  2:532413/150‑5 (MQ=255)
                gCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTg                                                                                                                                      >  2:1039836/1‑150 (MQ=255)
                gCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTg                                                                                                                                      >  2:123960/1‑150 (MQ=255)
                  atatTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGTTTTACAATCGGCAttcatcaag                                                                                                                                    <  1:686601/150‑8 (MQ=255)
                  atatTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGcc                                                                                                                                    >  1:568851/1‑150 (MQ=255)
                      tGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATATGGTGCCATTg                                                                                                                                <  2:834728/150‑1 (MQ=255)
                       gTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTga                                                                                                                               <  1:1299614/150‑1 (MQ=255)
                        tCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTgaa                                                                                                                              <  2:657612/150‑1 (MQ=255)
                          ccGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTgaaga                                                                                                                            >  1:28659/1‑150 (MQ=255)
                           cGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTgaagaa                                                                                                                           >  1:900085/1‑150 (MQ=255)
                            gTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAAt                                                                                                                          <  2:997056/150‑1 (MQ=255)
                             ttAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATc                                                                                                                         <  1:90213/150‑1 (MQ=255)
                             ttAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATc                                                                                                                         <  1:1400061/150‑1 (MQ=255)
                                 tGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCActc                                                                                                                     <  1:1055790/150‑1 (MQ=255)
                                    tttCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTa                                                                                                                  >  1:875796/1‑150 (MQ=255)
                                    tttCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTa                                                                                                                  >  1:77765/1‑150 (MQ=255)
                                       cAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATgg                                                                                                               >  2:1495854/1‑150 (MQ=255)
                                       aaaaaTGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGGGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATgg                                                                                                               <  1:653211/149‑1 (MQ=255)
                                        aaaaTGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGTTTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGt                                                                                                              <  2:981994/150‑1 (MQ=255)
                                        aaaaTGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGt                                                                                                              >  1:884082/1‑150 (MQ=255)
                                        aaaaTGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGt                                                                                                              >  2:1050887/1‑150 (MQ=255)
                                        aaaaTGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGt                                                                                                              >  1:938566/1‑150 (MQ=255)
                                           aTGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTcag                                                                                                           >  2:1492786/1‑150 (MQ=255)
                                             gAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTcagca                                                                                                         <  2:633516/150‑1 (MQ=255)
                                              aGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTcagcag                                                                                                        <  2:15693/150‑1 (MQ=255)
                                               gCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGa                                                                                                       >  1:332202/1‑150 (MQ=255)
                                                  gTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGaaca                                                                                                    >  1:1396666/1‑150 (MQ=255)
                                                  gTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGaaca                                                                                                    <  2:318720/150‑1 (MQ=255)
                                                   tAAAAATGTGCGGATGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTATGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATGACTCTTATGGTCAGCAGaacaa                                                                                                   <  2:865279/150‑1 (MQ=255)
                                                   tAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGaacaa                                                                                                   <  2:243259/150‑1 (MQ=255)
                                                       aaTGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACggg                                                                                               <  2:1149637/150‑1 (MQ=255)
                                                         tGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGggcg                                                                                             <  1:1492786/150‑1 (MQ=255)
                                                          gTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGggcgg                                                                                            >  1:914271/1‑150 (MQ=255)
                                                            ccGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCt                                                                                          >  2:770157/1‑150 (MQ=255)
                                                             cGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCtt                                                                                         >  2:790111/1‑150 (MQ=255)
                                                              ggAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTc                                                                                        >  1:1535805/1‑150 (MQ=255)
                                                               gAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCg                                                                                       >  1:312421/1‑150 (MQ=255)
                                                                 aGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCt                                                                                     <  2:568851/150‑1 (MQ=255)
                                                                    cacCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTggc                                                                                  <  1:123960/150‑1 (MQ=255)
                                                                    cacCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTggc                                                                                  >  1:1466428/1‑150 (MQ=255)
                                                                     acCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTggcg                                                                                 <  1:222654/150‑1 (MQ=255)
                                                                      ccAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTggcgg                                                                                <  2:916840/150‑1 (MQ=255)
                                                                      ccAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTggcgg                                                                                >  1:388975/1‑150 (MQ=255)
                                                                       cAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTggcggc                                                                               >  1:742450/1‑150 (MQ=255)
                                                                       cAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTggcggc                                                                               >  2:1412422/1‑150 (MQ=255)
                                                                        aaaaaGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGgcggcg                                                                              >  2:21351/1‑150 (MQ=255)
                                                                                 tCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCt                                                                     <  2:1241664/150‑1 (MQ=255)
                                                                                   gCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTg                                                                   >  1:485501/1‑150 (MQ=255)
                                                                                    cAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGACGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGCCGCTGGCTTGt                                                                  >  2:1063453/1‑150 (MQ=255)
                                                                                     aGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTc                                                                 <  2:1466428/150‑1 (MQ=255)
                                                                                      gggAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCg                                                                <  1:770157/150‑1 (MQ=255)
                                                                                      gggAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCg                                                                <  2:875796/150‑1 (MQ=255)
                                                                                            tATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCtga                                                          <  2:119645/150‑1 (MQ=255)
                                                                                             aTGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTggcggc                                                                               >  2:416907/1‑128 (MQ=255)
                                                                                             aTGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTggcggc                                                                               <  1:416907/128‑1 (MQ=255)
                                                                                              tGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCtgatg                                                        <  2:232157/150‑1 (MQ=255)
                                                                                                ccggggCGGAAGTGTAAGGGGATAACGCGTGACTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGACGAATCACTCTTATGGGCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCGTGTCGACCTGATGAg                                                      <  2:965898/146‑1 (MQ=255)
                                                                                                cAGTGGCTGAAGTGTAAGGTGATAACGCTTGATTTCTTGATGAATGCGGATTGTAAATCGGCATTTGGTGCCATTGAAGAATCAGTCTTATGGTCAGCAGTAGAACGGTCGGCGTCGGTGGCGGCGACGCTGGCTTGTCGACCTGATGAg                                                      <  2:443917/150‑1 (MQ=255)
                                                                                                      cGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCg                                                >  2:424478/1‑150 (MQ=255)
                                                                                                       ggAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGCCGCTGGCTTGTCGACCTGATGAGGGACCgg                                               >  1:1217691/1‑150 (MQ=255)
                                                                                                           gtgtAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGtgtg                                           >  2:10312/1‑150 (MQ=255)
                                                                                                            tgtAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCTCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTgg                                          >  2:181892/1‑150 (MQ=255)
                                                                                                             gtAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGa                                         >  2:986696/1‑150 (MQ=255)
                                                                                                             gtAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGa                                         >  2:838299/1‑150 (MQ=255)
                                                                                                             gtAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGa                                         >  1:1026531/1‑150 (MQ=255)
                                                                                                             gtAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGa                                         >  1:278611/1‑150 (MQ=255)
                                                                                                             gtAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGa                                         >  1:446640/1‑150 (MQ=255)
                                                                                                                  gTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTc                                    <  1:1412422/150‑1 (MQ=255)
                                                                                                                  gTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTc                                    >  1:957966/1‑150 (MQ=255)
                                                                                                                     atatCGCGTGCTTTGTTGATGAATGCCGTTTGTAAAACCGCATTTGGTTCAATTGAAGAATAACTCTTATTGTCTGCTGAACAACGGGCGGCTTCGGTGGCGGCGACGCTGGCTGGTCGACATGATGAGGGACCGGTGTGGATCTTCGGc                                 <  2:1358320/150‑1 (MQ=255)
                                                                                                                      tAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCt                                <  2:253004/150‑1 (MQ=255)
                                                                                                                        aCGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTAt                              <  2:884082/150‑1 (MQ=255)
                                                                                                                         cgcgTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATg                             <  2:938566/150‑1 (MQ=255)
                                                                                                                         cgcgTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATg                             <  1:872717/150‑1 (MQ=255)
                                                                                                                         cgcgTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATg                             >  1:262776/1‑150 (MQ=255)
                                                                                                                           cgTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGa                           <  1:504792/150‑1 (MQ=255)
                                                                                                                                 ttCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTg                     <  2:77765/150‑1 (MQ=255)
                                                                                                                                 ttCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTg                     <  1:1039836/150‑1 (MQ=255)
                                                                                                                                    ttGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGa                                         <  2:1285460/127‑1 (MQ=255)
                                                                                                                                     caccaaaTGCCGTTTTACAATCGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGAtgt                 >  2:607970/6‑150 (MQ=255)
                                                                                                                                     caccaaaTGCCGTTTTACAATCGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGAtgt                 >  2:848606/6‑150 (MQ=255)
                                                                                                                                      accaaaTGCCGTTTTACAATCGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGAtgtg                >  1:1106949/5‑150 (MQ=255)
                                                                                                                                       ccaaaTGCCGTTTTACAATCGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGATGTgg               >  1:646284/4‑150 (MQ=255)
                                                                                                                                       atgaATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGATGTgg               <  2:76815/150‑1 (MQ=255)
                                                                                                                                       atgaATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGATGTgg               <  1:870560/150‑1 (MQ=255)
                                                                                                                                         gaATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGaa             >  2:823626/1‑150 (MQ=255)
                                                                                                                                         gaATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGaa             >  2:195490/1‑150 (MQ=255)
                                                                                                                                         aaaTGCCGTTTTACAATCGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGaa             >  2:195360/2‑150 (MQ=255)
                                                                                                                                           aTGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATc           <  1:1050887/150‑1 (MQ=255)
                                                                                                                                            cgatctttGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATcc          <  1:935391/144‑1 (MQ=255)
                                                                                                                                                 aTTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATCCGGCAc     <  1:1194748/150‑1 (MQ=255)
                                                                                                                                                  ttGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTAAATATGAAATcaagatcg    >  2:935391/1‑143 (MQ=255)
                                                                                                                                                    gTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATCCGGCACTgg  >  1:966233/1‑150 (MQ=255)
                                                                                                                                                    gTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATCCGGCACTgg  >  1:1108708/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AAGTGGTCCGTGGATTGCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATCCGGCACTGG  >  NC_000913/1272393‑1272690

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