Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,966,050 C→A 5.9% G48C (GGC→TGC)  flhB ← flagellin export apparatus, substrate specificity protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,966,0500CA 5.9% 233.4 / 2.1 68G48C (GGC→TGC) flhBflagellin export apparatus, substrate specificity protein
Reads supporting (aligned to +/- strand):  ref base C (31/33);  new base A (2/2);  total (33/35)
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.36e-01

GCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGT  >  NC_000913/1965912‑1966184
                                                                                                                                          |                                                                                                                                      
gcagcaGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCTGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCa                                                                                                                                       <  1:412618/140‑1 (MQ=255)
  agcagcGCCAGCATGGCTTCTCTGATCAGCTGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCaaa                                                                                                                                     <  2:1563391/140‑1 (MQ=255)
  agcagcGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCaaa                                                                                                                                     <  2:1091710/140‑1 (MQ=255)
  agcagcGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCaaa                                                                                                                                     <  1:133626/140‑1 (MQ=255)
  agcagcGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCaaa                                                                                                                                     <  1:362538/140‑1 (MQ=255)
   gcagcGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGTCACACCGCCAAAc                                                                                                                                    <  2:1115409/140‑1 (MQ=255)
   gcagcGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAAc                                                                                                                                    <  2:642338/140‑1 (MQ=255)
   gcagcGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAAc                                                                                                                                    <  2:1392555/140‑1 (MQ=255)
   gcagcGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAAc                                                                                                                                    <  2:1257550/140‑1 (MQ=255)
        gcCAGCTTGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGAt                                                                                                                               <  1:678127/140‑1 (MQ=255)
        gcCAGCATTGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGAt                                                                                                                               <  2:915021/140‑1 (MQ=255)
        gcCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGAt                                                                                                                               <  2:262568/140‑1 (MQ=255)
        gcCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGAt                                                                                                                               <  2:262520/140‑1 (MQ=255)
        gcCAGCATGGCTTATCTGATCAGCAGAATATTCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGAt                                                                                                                               <  2:915029/140‑1 (MQ=255)
        gcCAGAATGGCTTCTCTGATCAGCAGAATAATCTGCCCTAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCTTAAATCGACGGGCCAGCGACACACCGCCAAACCAGAt                                                                                                                               <  2:834401/140‑1 (MQ=255)
                  cTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACg                                                                                                                     >  2:1177339/1‑140 (MQ=255)
                  cTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACg                                                                                                                     >  2:276541/1‑140 (MQ=255)
                             agcagAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCaa                                                                                                          <  2:1424001/140‑1 (MQ=255)
                             agcagAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCaa                                                                                                          >  2:221984/1‑140 (MQ=255)
                             agcagAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCaa                                                                                                          >  2:222053/1‑140 (MQ=255)
                               cagAATAATCTGCCCGAGGATCCGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACACTCGACGGGCCCGCGACACACCGCCCACCCAGATAACACTAACGCCCACCAGCaaaa                                                                                                        >  2:52827/1‑140 (MQ=255)
                                    taatCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATcagc                                                                                                   >  1:447940/1‑140 (MQ=255)
                                       tCTGCCCGAGGATCAGATTCGGGTGATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACTATCGACGGGCCAGCGACACACCGCCAAACCCGATAACAATAACGGCCACCAGCAAAATCAGCAGt                                                                                                <  2:628223/140‑1 (MQ=255)
                                         tGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGa                                                                                              >  1:549885/1‑140 (MQ=255)
                                         tGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGa                                                                                              >  2:508730/1‑140 (MQ=255)
                                            ccGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGt                                                                                           <  1:1327092/140‑1 (MQ=255)
                                            ccGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGt                                                                                           >  2:1017718/1‑140 (MQ=255)
                                               aGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAg                                                                                        >  2:834547/1‑140 (MQ=255)
                                                            ggTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAAcgc                                                                           >  1:615524/1‑140 (MQ=255)
                                                            ggTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAAcgc                                                                           >  1:801097/1‑140 (MQ=255)
                                                                       ataCTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGccc                                                                >  1:100184/1‑140 (MQ=255)
                                                                       ataCTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGccc                                                                >  1:1039060/1‑140 (MQ=255)
                                                                       ataCTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGccc                                                                >  1:1038215/1‑140 (MQ=255)
                                                                       ataCTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGccc                                                                >  2:1641812/1‑140 (MQ=255)
                                                                       ataCTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGccc                                                                >  1:1038039/1‑140 (MQ=255)
                                                                       ataCTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGccc                                                                >  1:1267316/1‑140 (MQ=255)
                                                                       ataCTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGccc                                                                >  1:1507489/1‑140 (MQ=255)
                                                                             tgATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTc                                                          >  1:749848/1‑140 (MQ=255)
                                                                             tgATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTc                                                          >  1:342693/1‑140 (MQ=255)
                                                                             tgATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTc                                                          >  1:665869/1‑140 (MQ=255)
                                                                              gATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCAAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTcc                                                         >  1:736680/1‑140 (MQ=255)
                                                                              gATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCAAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTcc                                                         >  1:736584/1‑140 (MQ=255)
                                                                                  aaaaTGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCcgcg                                                     >  2:58412/1‑140 (MQ=255)
                                                                                        cAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCAAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCttttt                                               <  2:736680/140‑1 (MQ=255)
                                                                                        cAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCAAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCttttt                                               <  2:736584/140‑1 (MQ=255)
                                                                                          gCCCAGCGGAGAGCATGCCCTACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCt                                             <  1:1065280/140‑1 (MQ=255)
                                                                                          gCCCAGCGGAGAGCATGCCCTACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCt                                             <  1:1364973/140‑1 (MQ=255)
                                                                                          gCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCt                                             >  1:1216485/1‑140 (MQ=255)
                                                                                          gCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCt                                             <  2:337058/140‑1 (MQ=255)
                                                                                          gCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCt                                             >  1:1216611/1‑140 (MQ=255)
                                                                                          gCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTGTTCt                                             >  1:1216470/1‑140 (MQ=255)
                                                                                          ccccAGCGCCGAGCAGGCCCCACACACGACCGGCCACCTTCACACCCCCAAACCAGATAACACATACGCCCACCAGCAAAATCATCAGTGAGTTCAGTTCACGGGAACGCGGGTTTTGCCCCTCTTCCCGCGCTTTTTCt                                             <  2:1672640/139‑1 (MQ=255)
                                                                                            ccAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAg                                           >  1:3302/1‑140 (MQ=255)
                                                                                            ccAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGAGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAg                                           <  2:1553627/140‑1 (MQ=255)
                                                                                              aGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTc                                         >  1:1493866/1‑140 (MQ=255)
                                                                                              aGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAAAGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTc                                         >  1:1494041/1‑140 (MQ=255)
                                                                                                 ggAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGt                                      <  1:284529/140‑1 (MQ=255)
                                                                                                 ggAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGt                                      <  2:865876/140‑1 (MQ=255)
                                                                                                  gagaGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTg                                     <  1:1388735/140‑1 (MQ=255)
                                                                                                    gagCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAg                                   <  1:1177827/140‑1 (MQ=255)
                                                                                                          gCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTgg                             >  1:151670/1‑140 (MQ=255)
                                                                                                              gACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGc                         <  2:1145697/140‑1 (MQ=255)
                                                                                                               aCAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCt                        >  2:547888/1‑140 (MQ=255)
                                                                                                                cAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCtt                       <  2:330873/140‑1 (MQ=255)
                                                                                                                                gACACACCGCCAACCCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCACTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGct       <  2:848622/140‑1 (MQ=255)
                                                                                                                                gACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGct       <  2:1064603/140‑1 (MQ=255)
                                                                                                                                     acCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGt  <  2:1451818/140‑1 (MQ=255)
                                                                                                                                     acCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGt  <  2:1361955/140‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                      
GCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGT  >  NC_000913/1965912‑1966184

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