Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,761,620 T→G 8.1% K641Q (AAG→CAG)  abpB ← CP4‑57 prophage; putative helicase YfjK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,761,6200TG 8.1% 160.9 / 7.0 67K641Q (AAG→CAG) abpBCP4‑57 prophage; putative helicase YfjK
Reads supporting (aligned to +/- strand):  ref base T (29/31);  new base G (2/4);  total (31/36)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.29e-01

AATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATAC‑CACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGT  >  NC_000913/2761478‑2761762
                                                                                                                                               |                                                                                                                                              
aaTTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACATGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATAC‑CACGCtt                                                                                                                                                >  2:147175/1‑143 (MQ=255)
  ttCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATGGACTTTGCTGCGGATAA‑CAGCCtttg                                                                                                                                              >  2:184772/1‑142 (MQ=255)
     cAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTGTAC‑CACGCTTTACAg                                                                                                                                           >  1:390724/1‑143 (MQ=255)
        cgaGTCTGCTAGTCGGTGTAACGTCGGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATAC‑CACGCTTTACAGCgt                                                                                                                                        <  2:4046/141‑1 (MQ=255)
         taGTCTACCAGTCGGTGTAACGTCTGAATGGTAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATAC‑CACGCTGTAAAGCgtg                                                                                                                                       <  2:143968/142‑1 (MQ=255)
         gagccTACCAGTAGGTGGAGCGTCTGAGTGGGAATGGCAAGCACTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCACCTTGACATGTGTGTAATCGACTTTGCTGCGTATAC‑CACGCTTTACAGCgtg                                                                                                                                       <  2:97889/139‑1 (MQ=255)
         gaGTCGACCAGTAGGTGCGACGTCTGGATGGGAATGCCAATCTCTTCTAATGCAGCGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGCCATGTGTTTACTCGACTGGGCTGCGTATAC‑CACGCTGTAAAGCgtg                                                                                                                                       <  1:816781/142‑1 (MQ=255)
         aaGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATAC‑CACGCTTTACAGCgtg                                                                                                                                       <  1:137016/143‑1 (MQ=255)
         aaGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATAC‑CACGCTGTAAAGCgtg                                                                                                                                       <  1:724974/143‑1 (MQ=255)
           ggcTACCAGTCGGTGTATCGTCTGAATGGGAATGCCAATCTCTTCTAAGGCTTTGACACCTGGTGGTAGATGGCAACTTTCGAACGCCAATTTGACATGTGTGTAATCGGCTTTGCTGCGTATAC‑CACGCTTTACAGCgtggt                                                                                                                                     <  2:682742/141‑1 (MQ=255)
            tCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGCTGATTTCGCTCGCCAACCTGACATGTGTGTAATCGACTTTGTTACGTATAA‑CAAGCTTTACAGCgtggtg                                                                                                                                    >  1:207384/1‑143 (MQ=255)
             cTAGCAGTCGGTGTAACGTCGCAATGGGAATGCCGATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAATCTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATAC‑CACGCTTTACAGCGTGgtgt                                                                                                                                   <  2:407835/143‑1 (MQ=255)
                    tCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATAC‑CACGCTTTACAGCGTGGTGTTTAACTa                                                                                                                            <  2:125939/143‑1 (MQ=255)
                       gtgtAACGTCTGAATGGAAATGCCAATCTATTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTCGCTGCGTATAC‑CACGCTTTACAGCGTGGTGTCGTACTACAt                                                                                                                         >  2:782735/1‑143 (MQ=255)
                            aCGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATAC‑CACGCTTTACAGCGTGGTGTTTAACTACAtcctcc                                                                                                                    <  2:746672/143‑1 (MQ=255)
                             cGTCTAAATGGGCATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATAC‑CACGCTTTACAGCGTGGTGTTTAACTACATCCTCCa                                                                                                                   <  2:756510/143‑1 (MQ=255)
                               tCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATAC‑CACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATc                                                                                                                 >  2:334094/1‑143 (MQ=255)
                                  gAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATAC‑CACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAg                                                                                                              <  1:143068/143‑1 (MQ=255)
                                             ccAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATAC‑CACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCtt                                                                                                   <  1:270716/143‑1 (MQ=255)
                                              ctaTCTCTTCTAAGGCATTTACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATAC‑CACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTta                                                                                                  <  2:454872/141‑1 (MQ=255)
                                              cAATCTCTTTCAAGGCATTGACACCTGGCGGTAGAGGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATAC‑CACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTta                                                                                                  <  1:389474/143‑1 (MQ=255)
                                              cAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGTAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATAC‑CACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTta                                                                                                  >  1:425940/1‑143 (MQ=255)
                                                  ctctTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATAC‑CACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTtaggta                                                                                              >  1:779595/1‑143 (MQ=255)
                                                    tttgttgggCATCGACACGGGGGGGGAGGTGGAAGCCTTCGACCGCCAATTTGACATGTGTGTAATCGTCTTTGCTGCGTATAC‑CACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTtaggtagg                                                                                            <  1:1220/136‑1 (MQ=255)
                                                    tttctAAGGCCTTGACACCGGGTTGTAGATGGAAACGTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATAC‑CACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTtaggtagg                                                                                            <  1:46976/142‑1 (MQ=255)
                                                     ttctAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATAC‑CACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAggg                                                                                           <  2:558729/143‑1 (MQ=255)
                                                       ctAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATAC‑CACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGta                                                                                         <  2:21844/143‑1 (MQ=255)
                                                        tAAGGCATTGACACCTGGTGGTAGATGGAAAATTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATAC‑CACGCTTTACAGCGTGGTGTTTAACTACATCCCCCATCAAGCTGGTCACCTTCGGCAgggaag                                                                                        >  1:867407/1‑140 (MQ=255)
                                                          aGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATAC‑CACGCTTTACAGCGTGGTGTCTAACTACATCCTCCATCACGCTGAGCACCTTAGGCAGGGtatac                                                                                      >  2:160991/1‑142 (MQ=255)
                                                              aTTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATAC‑CACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTCGGCAGGGTCTAGCCGa                                                                                  >  2:711933/1‑143 (MQ=255)
                                                              aTTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGAGTGTAAACGACTATGCTGCGGATAC‑CACGATTTACAGCGTGGTGTTTAACTACATCCTCGATCGGGCTGAGCACCATAGGACGGGTATAGCCGa                                                                                  >  2:250329/1‑143 (MQ=255)
                                                               ttGACACCTGGTGGTAGATGGAAACTTGCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATAC‑CACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGaa                                                                                 <  1:713899/143‑1 (MQ=255)
                                                                    acCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATAC‑CACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTAt                                                                            <  2:514928/143‑1 (MQ=255)
                                                                     ccTGGTGGTAGAGGGAAACCTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATAC‑CACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTAtt                                                                           <  1:275508/143‑1 (MQ=255)
                                                                      cTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATAC‑CACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTg                                                                          >  2:781402/1‑143 (MQ=255)
                                                                      cTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATAC‑CACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTg                                                                          >  2:207442/1‑143 (MQ=255)
                                                                      cTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGATGCGTATAC‑CACGCTTTACAGCGTGGGGTTTAACTACATCCTCGATCACACTGACAACCTCAGGTAGGGTATACCCGAAGCTAtag                                                                          >  2:222204/1‑141 (MQ=255)
                                                                           ggtAGATGGAAACTTCCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATAC‑CACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTAt                                                                     <  1:542452/143‑1 (MQ=255)
                                                                             tAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATAC‑CACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGa                                                                   <  1:334094/143‑1 (MQ=255)
                                                                                 tGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATAC‑CACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCtt                                                               >  1:78066/1‑143 (MQ=255)
                                                                                 tGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATAC‑CACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCtt                                                               <  1:301665/143‑1 (MQ=255)
                                                                                       cTTGCGAACGCCAATTTGACATGTGTGTAATCGACGTTGCTGCGTATACTGGCGCCAT‑CAACTTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATc                                                         <  1:418354/143‑1 (MQ=255)
                                                                                        tttCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATAC‑CACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCaacttg                                                        >  1:411022/1‑139 (MQ=255)
                                                                                        tttCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATAC‑CACGCTTTACAGCGTGGTGTTTAACTACATCCTCAATCCAGCTGATGACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCa                                                        <  1:768232/143‑1 (MQ=255)
                                                                                            gAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATAC‑CACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGTGCGCCTTAGGTTGGGTATATCCTAAGGTATTTGTTATTATCTTCAAATAAttat                                                    >  1:125585/1‑141 (MQ=255)
                                                                                             aaCGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATAC‑CACGCTTTACAGCGTGGTGTTTAACTACATCCTCAATAAAGCAAAACAAAATAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTGAAATCAttgtt                                                   >  1:186048/1‑143 (MQ=255)
                                                                                               cGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATAC‑CACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGa                                                 <  1:446000/143‑1 (MQ=255)
                                                                                                gCCAAGTTGACATGTGTGTAATCGACTTTGCTGCGTATAC‑CACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAg                                                <  1:259112/143‑1 (MQ=255)
                                                                                                 ccAATTTGACATGTGTGTAATCGACTTTGCTGCGTATAC‑CACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAgg                                               <  1:194591/143‑1 (MQ=255)
                                                                                                            tgtgtgTAATCGACTTTGCTGCGTATAC‑CACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAag                                    >  1:10673/1‑143 (MQ=255)
                                                                                                             gttggtAATCGACTTTGCTGCGTCTAC‑CACGCGTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTCAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAaga                                   <  2:490184/139‑1 (MQ=255)
                                                                                                             gtgtgtAATCGACTTTGCTGCGTATAC‑CACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAaga                                   <  2:500848/143‑1 (MQ=255)
                                                                                                                tgtAATCGACTTTGCTGCGTATAC‑CACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTa                                >  1:658269/1‑143 (MQ=255)
                                                                                                                tgtAATCGACTTTGCTGCGTATAC‑CACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTGCACATCATTGTTCAGGGCAATGGAAAGAGTa                                >  2:542328/1‑143 (MQ=255)
                                                                                                                  tAATCGACTGTGCTGCGTATAC‑CACGCTGGCCAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTGAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACg                              <  1:77675/143‑1 (MQ=255)
                                                                                                                      cGACTTTGCTGCGTATAC‑CACGCTTTACAGCGTGGTGTTTAAGTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTtc                          <  2:246839/143‑1 (MQ=255)
                                                                                                                      cGACTTTGCTGCGTATAC‑CACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGGACGTTtc                          >  1:855230/1‑143 (MQ=255)
                                                                                                                       gACTTTGCTGCGTATAC‑CACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGGGGGCTATGTAAAGAGTATGTTgct                         >  1:245509/1‑143 (MQ=255)
                                                                                                                         cTTTGCTGCGTATAC‑CACGCTTTACAGCGTGGTGTTTAACTACATCCTACATCAAGATGAGCACCTTAGGTAGGGTATAGAAGAAGGTATTGGTTATGAGCTTCAAATGATTGTTGAGGGAAATGGAAAGAGTAAGtttataa                       >  2:244110/1‑139 (MQ=255)
                                                                                                                          tttGCTGCGTATAC‑CACGCTTTTCAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTccc                      <  2:912961/143‑1 (MQ=255)
                                                                                                                          tttGCTGCGTATAC‑CACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGACGGTATTGGTTATGAGCTTCAAATCATTGTTGTGGCGAGTGGAATGAGTACGATTTTccc                      >  1:702801/1‑143 (MQ=255)
                                                                                                                          tttGCTGCGTATAC‑CACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGCGGCCACTGGAAAGGGTACGTTTCTccc                      >  1:280128/1‑143 (MQ=255)
                                                                                                                                 cGTATAC‑CACGCTTTACAGCGTGGTGTTTAACTACACCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGa               <  1:147175/143‑1 (MQ=255)
                                                                                                                                     caC‑CACGCTGTAAAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGACGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATc           >  2:306231/2‑143 (MQ=255)
                                                                                                                                     caC‑CACGCTGTAAAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCACCTTCAGAAATc           >  2:30326/2‑143 (MQ=255)
                                                                                                                                      aC‑CACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATcc          >  1:570731/1‑143 (MQ=255)
                                                                                                                                       c‑cACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCGATGGAAAGAGTACGTTTCTCCCCTTCAGAAATccc         >  2:611425/1‑142 (MQ=255)
                                                                                                                                         cACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTCAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTGAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCAc        >  2:188615/1‑143 (MQ=255)
                                                                                                                                          aCGCTTTACAGCGTGGTGTTTAAATACATACACCAGTAAGCTGAACACGTTAGGTAGGGTATAGGCAAAGAAATTAGGTATGAGCATCAAATCAGTGTTGAGGGCAGTGGAAATAGTACGTTTCTACCCTTCAGAAATCCACt       >  2:272825/1‑143 (MQ=255)
                                                                                                                                              tttACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCACCTTCAGAAATCCACTgacc   >  2:666644/1‑140 (MQ=255)
                                                                                                                                               ttACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATGGTTGAGGGCAATGGAAAGAGTGCGTTTCACGCCTTCAGAAATCCACTGTCGt  >  1:27490/1‑143 (MQ=255)
                                                                                                                                               |                                                                                                                                              
AATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATAC‑CACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGT  >  NC_000913/2761478‑2761762

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