Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,761,623 C→A 7.2% V640L (GTA→TTA)  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,6230CA 7.2% 190.2 / 4.9 69V640L (GTA→TTA) abpBCP4‑57 prophage; putative helicase YfjK
Reads supporting (aligned to +/- strand):  ref base C (31/33);  new base A (2/3);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.42e-01

CAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATAC‑CACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGT  >  NC_000913/2761483‑2761762
                                                                                                                                             |                                                                                                                                           
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:301665/143‑1 (MQ=255)
                                                                            tGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATAC‑CACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCtt                                                               >  1:78066/1‑143 (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)
                                                                                                                                             |                                                                                                                                           
CAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATAC‑CACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGT  >  NC_000913/2761483‑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: ‑