Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,370,155 A→G intergenic (+152/‑61) pspE → / → ycjM thiosulfate:cyanide sulfurtransferase (rhodanese)/alpha amylase catalytic domain family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,370,1550AG98.9% 336.6 / ‑4.0 89intergenic (+152/‑61)pspE/ycjMthiosulfate:cyanide sulfurtransferase (rhodanese)/alpha amylase catalytic domain family protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (52/36);  total (52/37)
Fisher's exact test for biased strand distribution p-value = 4.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.05e-01

AATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCATTTACAGAAACTTGTAACGCGT  >  NC_000913/1370007‑1370299
                                                                                                                                                    |                                                                                                                                                
aaTTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCa                                                                                                                                                <  2:334847/151‑1 (MQ=255)
  ttCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGt                                                                                                                                              >  2:1403475/1‑151 (MQ=255)
    cACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGttt                                                                                                                                            >  2:104850/1‑151 (MQ=255)
    cACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGtt                                                                                                                                             <  1:1451472/150‑1 (MQ=255)
          gAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTc                                                                                                                                      <  2:52853/151‑1 (MQ=255)
           aGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCt                                                                                                                                     <  1:898556/151‑1 (MQ=255)
             cggcgCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGc                                                                                                                                   <  1:1275143/151‑1 (MQ=255)
              ggcgCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCg                                                                                                                                  >  2:33759/1‑151 (MQ=255)
                 gcTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCgg                                                                                                                                 >  1:510205/1‑149 (MQ=255)
                 gcTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGa                                                                                                                                <  1:535529/150‑1 (MQ=255)
                   ttATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGagt                                                                                                                              >  2:99571/1‑150 (MQ=255)
                   ttATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGagt                                                                                                                              >  2:160029/1‑150 (MQ=255)
                    tATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGagta                                                                                                                             <  1:122557/150‑1 (MQ=255)
                     aTGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGagtagt                                                                                                                           >  2:1609910/1‑151 (MQ=255)
                        cgcCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGtctc                                                                                                                        <  2:1216023/151‑1 (MQ=255)
                         gcCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGtctct                                                                                                                       <  2:262931/151‑1 (MQ=255)
                          ccGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGtctctc                                                                                                                      <  1:216763/151‑1 (MQ=255)
                            ttttttttttCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGt                                                                                                                    >  1:36412/2‑151 (MQ=255)
                            gTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGt                                                                                                                    >  2:934857/1‑151 (MQ=255)
                                  ttttCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATg                                                                                                              >  2:511012/1‑151 (MQ=255)
                                  ttttCTGTCACACCTTATTTACATCACCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTATCTCGTTTCAt                                                                                                               >  2:1307112/1‑150 (MQ=255)
                                    ttCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATggg                                                                                                            >  2:346720/1‑151 (MQ=255)
                                     tCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATggg                                                                                                            >  1:251230/1‑150 (MQ=255)
                                        gTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGa                                                                                                           <  1:737355/148‑1 (MQ=255)
                                        gTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCg                                                                                                        <  2:770165/151‑1 (MQ=255)
                                          cacaCCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCt                                                                                                      <  1:346695/151‑1 (MQ=255)
                                             acCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTAcc                                                                                                   >  1:443736/1‑151 (MQ=255)
                                              ccTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTAcc                                                                                                   >  1:1072239/1‑150 (MQ=255)
                                               cTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCAc                                                                                                 >  2:8676/1‑151 (MQ=255)
                                                 tATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACgg                                                                                               >  1:1064236/1‑151 (MQ=255)
                                                    ttACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGaa                                                                                             <  2:1197846/150‑1 (MQ=255)
                                                     tACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGaaa                                                                                            <  2:366908/150‑1 (MQ=255)
                                                     tACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGaa                                                                                             >  2:733549/1‑149 (MQ=255)
                                                       cATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGc                                                                                         >  1:1600789/1‑151 (MQ=255)
                                                       cATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGc                                                                                         >  2:1235456/1‑151 (MQ=255)
                                                        aTCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCa                                                                                        <  1:1440471/151‑1 (MQ=255)
                                                           cccATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAAc                                                                                      >  2:461346/1‑150 (MQ=255)
                                                             cATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATg                                                                                   <  2:1130007/151‑1 (MQ=255)
                                                                aGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGaaa                                                                                <  2:635850/151‑1 (MQ=255)
                                                                  attattTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAAc                                                                               >  2:1330028/1‑150 (MQ=255)
                                                                  attattTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAAc                                                                               >  1:1103536/1‑150 (MQ=255)
                                                                    tattTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAAc                                                                               >  1:732919/1‑148 (MQ=255)
                                                                     attTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGa                                                                            <  2:1127296/150‑1 (MQ=255)
                                                                       ttGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGaaa                                                                          >  2:480776/1‑150 (MQ=255)
                                                                               gCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGAt                                                                 >  1:189121/1‑151 (MQ=255)
                                                                                cTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGAtt                                                                >  2:74722/1‑151 (MQ=255)
                                                                                  cacaAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATTGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGAtt                                                                <  1:338953/149‑1 (MQ=255)
                                                                                  cacaAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGAtt                                                                <  1:278139/149‑1 (MQ=255)
                                                                                  cacaAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTa                                                               <  2:1692418/150‑1 (MQ=255)
                                                                                  cacaAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTa                                                               <  2:74368/150‑1 (MQ=255)
                                                                                  cacaAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTAc                                                              >  1:567114/1‑151 (MQ=255)
                                                                                     aaaTACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTg                                                           >  1:1642215/1‑151 (MQ=255)
                                                                                      aaTACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTgg                                                          <  1:1605478/151‑1 (MQ=255)
                                                                                       aTACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGa                                                         <  2:191389/151‑1 (MQ=255)
                                                                                          cGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACg                                                       >  2:1066044/1‑150 (MQ=255)
                                                                                              ttttCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGaaa                                                    <  1:1307018/149‑1 (MQ=255)
                                                                                               tttCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCt                                                 >  1:1440734/1‑151 (MQ=255)
                                                                                                ttCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCt                                                 >  1:1328126/1‑150 (MQ=255)
                                                                                                ttCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCt                                                 <  1:145143/150‑1 (MQ=255)
                                                                                                 gcacTGGGAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCACTTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTAc                                               >  1:201255/4‑151 (MQ=255)
                                                                                                  cccTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACg                                              <  1:1643352/151‑1 (MQ=255)
                                                                                                    cTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACgg                                             >  1:1557269/1‑150 (MQ=255)
                                                                                                      ggTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTgg                                          <  1:470046/151‑1 (MQ=255)
                                                                                                         aaaaaaTGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGaa                                        <  2:1103621/150‑1 (MQ=255)
                                                                                                          aaaaaTGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAAc                                       >  2:648389/1‑150 (MQ=255)
                                                                                                                  tttCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGc                              >  2:462471/1‑151 (MQ=255)
                                                                                                                       tGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTc                         >  1:422172/1‑151 (MQ=255)
                                                                                                                        gCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACt                          <  1:787185/149‑1 (MQ=255)
                                                                                                                         cGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCAt                       >  2:1167862/1‑151 (MQ=255)
                                                                                                                         cGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCAt                       >  1:582272/1‑151 (MQ=255)
                                                                                                                          gtgtCTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCAt                       <  2:205551/150‑1 (MQ=255)
                                                                                                                           tGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCAttt                     >  1:1065229/1‑151 (MQ=255)
                                                                                                                               tAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCATTTAc                   >  2:301644/1‑149 (MQ=255)
                                                                                                                                 aaaCCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCATTTACAGa                >  1:1650747/1‑150 (MQ=255)
                                                                                                                                   aCCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCATTTACAGAAAc             >  2:656747/1‑151 (MQ=255)
                                                                                                                                    cccTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCATTTACAGAAACt            <  1:733492/151‑1 (MQ=255)
                                                                                                                                     ccTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCATTTACAGAAACtt           <  2:403936/151‑1 (MQ=255)
                                                                                                                                       ttGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCATTTACAGAAACTTGt         >  2:490691/1‑151 (MQ=255)
                                                                                                                                        tgtgCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCATTTACAGAAACTTg          <  2:272169/149‑1 (MQ=255)
                                                                                                                                        tgtgCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCATTTACAGAAACTTGt         >  1:1295932/1‑150 (MQ=255)
                                                                                                                                           gCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCATTTACAGAAACTTGTAAcg     >  2:898130/1‑151 (MQ=255)
                                                                                                                                            cTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCATTTACAGAAACTTGTAAcgc    >  2:1380462/1‑151 (MQ=255)
                                                                                                                                            cTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCATTTACAGAAACTTGTAAcgc    >  1:719770/1‑151 (MQ=255)
                                                                                                                                             tCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCATTTACAGAAACTTGTAAcgcg   >  2:1097329/1‑151 (MQ=255)
                                                                                                                                             tCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCATTTACAGAAACTTGTAAcgcg   >  2:943863/1‑151 (MQ=255)
                                                                                                                                             tCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCATTTACAGAAACTTGTAAcgcg   >  2:971009/1‑151 (MQ=255)
                                                                                                                                              cAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCATTTACAGAAACTTGTAAcgcg   <  1:1647323/150‑1 (MQ=255)
                                                                                                                                              cAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCATTTACAGAAACTTGTAACGCGt  >  2:1164495/1‑151 (MQ=255)
                                                                                                                                              cAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCATTTACAGAAACTTGTAACGCGt  <  1:648332/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
AATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCATTTACAGAAACTTGTAACGCGT  >  NC_000913/1370007‑1370299

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