Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,933,401 A→C 5.5% V96V (GTT→GTG)  fucA ← L‑fuculose‑1‑phosphate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,933,4010AC 5.5% 291.0 / 4.0 99V96V (GTT→GTG) fucAL‑fuculose‑1‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base A (50/43);  new base C (2/4);  total (52/47)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

TGAGAGCCAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAGTTTTTCATA  >  NC_000913/2933256‑2933550
                                                                                                                                                 |                                                                                                                                                     
tGAGAGCCAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACgg                                                                                                                                                   <  2:183993/150‑1 (MQ=255)
 gagaGCCAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGc                                                                                                                                                  <  1:904994/150‑1 (MQ=255)
   gagCCAGCGCAACATGTTCAGAAAGTTCGCGGGTCCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGAGCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCAt                                                                                                                                                <  1:733925/150‑1 (MQ=255)
    agCCAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTTCATTGCACGGCagt                                                                                                                                               <  1:1684598/150‑3 (MQ=255)
    acacAGCGCAACATGTTCAGAAACTTCGCGTGTTCCAAAGGTCTCATAAGGCGAGATAGGAATATAATTACCTCGATAAGCCTCAATCATGTTGTGAATTTCGGGGATCGATCGTTTTATAATGGAAACTTCCGTTCATTGCACGGCAtt                                                                                                                                               <  2:1403565/147‑2 (MQ=255)
     cccAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGa                                                                                                                                              <  1:31978/149‑1 (MQ=255)
      ccAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTTCATTGCACGGCagttt                                                                                                                                             <  2:481747/150‑5 (MQ=255)
        aGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTCCCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATtg                                                                                                                                           <  1:1524152/150‑1 (MQ=255)
         gcgcAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATtgt                                                                                                                                          >  1:1258252/1‑150 (MQ=255)
         gcgcAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATtgt                                                                                                                                          <  2:1371847/150‑1 (MQ=255)
           gcAACATGTTCAGAAAGTTCGCGTGTTCCAAGGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGa                                                                                                                                        <  1:1475143/150‑1 (MQ=255)
           gcAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGa                                                                                                                                        <  2:1003866/150‑1 (MQ=255)
           gcAACATGTTCAGAAAGATCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGTGGATCTAACCGTAAAGAATGGAAACTGCCGTGCAATGAACGGCATCATTGTGa                                                                                                                                        >  2:1214815/1‑150 (MQ=255)
            cAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGaa                                                                                                                                       >  1:64746/1‑150 (MQ=255)
            cAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGaa                                                                                                                                       >  1:1531570/1‑150 (MQ=255)
               cATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGTCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTTTGaacaa                                                                                                                                    >  2:1512643/1‑150 (MQ=255)
                aTGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCACTGAACGGCATGATTGTGaacaaa                                                                                                                                   >  2:852904/1‑149 (MQ=255)
                 tGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAAcc                                                                                                                                  >  1:1360653/1‑150 (MQ=255)
                     cAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGt                                                                                                                              >  2:144703/1‑150 (MQ=255)
                     cAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGt                                                                                                                              >  2:1546395/1‑150 (MQ=255)
                     cAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAAAAATCGCGt                                                                                                                              >  2:40612/1‑150 (MQ=255)
                      aGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGtt                                                                                                                             >  1:1558149/1‑150 (MQ=255)
                       gAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTg                                                                                                                            >  2:1486963/1‑150 (MQ=255)
                              cgcgTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAAACGCGTTGGCATCCg                                                                                                                     >  1:116386/1‑150 (MQ=255)
                                 gtgtTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTc                                                                                                                  <  1:158959/150‑1 (MQ=255)
                                  tgtTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCt                                                                                                                 >  1:390615/1‑150 (MQ=255)
                                  tgtTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCt                                                                                                                 >  1:1471582/1‑150 (MQ=255)
                                  tgtTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCt                                                                                                                 >  1:1455424/1‑150 (MQ=255)
                                   gtTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTg                                                                                                                >  1:1546557/1‑150 (MQ=255)
                                      ccAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGACCGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCtt                                                                                                             >  1:1058955/1‑150 (MQ=255)
                                      ccAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATGCGGTCTGCtt                                                                                                             >  1:1064080/1‑150 (MQ=255)
                                      ccAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCtt                                                                                                             >  2:1113315/1‑150 (MQ=255)
                                      ccAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCtt                                                                                                             >  2:48532/1‑150 (MQ=255)
                                      ccAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCtt                                                                                                             >  2:359313/1‑150 (MQ=255)
                                      ccAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCtt                                                                                                             <  1:487425/150‑1 (MQ=255)
                                      ccAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCtt                                                                                                             >  1:1269683/1‑150 (MQ=255)
                                         aaGGTCGCATAATGCGCGCAAGGAATTGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATTGCGGGGATAGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGa                                                                                                          <  1:947219/150‑1 (MQ=255)
                                          aGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGAt                                                                                                         <  2:1210969/150‑1 (MQ=255)
                                          aGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGACGTGCAATGAACTGCATGATTGTGAACAACCGCGTTGGCATACGGTCTGCTTTGAt                                                                                                         >  2:947230/1‑150 (MQ=255)
                                           ggTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATa                                                                                                        >  1:649759/1‑150 (MQ=255)
                                           ggTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATa                                                                                                        <  1:1542841/150‑1 (MQ=255)
                                                cATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTg                                                                                                   >  1:1553381/1‑150 (MQ=255)
                                                  tAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGcc                                                                                                 >  1:1124073/1‑150 (MQ=255)
                                                          gcAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGaa                                                                                         >  1:1261410/1‑150 (MQ=255)
                                                            aagctaTAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGTATCGATCGGTAAAGAATCGAAACTGCCGTGCAATGCACTTCATGATTGTTAACAACCTCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAAc                                                                                       <  1:205666/145‑1 (MQ=255)
                                                            aaGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAAc                                                                                       >  1:503877/1‑150 (MQ=255)
                                                             aGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACg                                                                                      <  2:1449489/150‑1 (MQ=255)
                                                             aGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACg                                                                                      >  1:981828/1‑150 (MQ=255)
                                                             aGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACg                                                                                      <  2:917010/150‑1 (MQ=255)
                                                              ggAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGc                                                                                     >  1:1648867/1‑150 (MQ=255)
                                                               gAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGcc                                                                                    <  1:1688829/150‑1 (MQ=255)
                                                                aaTAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCa                                                                                   >  2:584670/1‑150 (MQ=255)
                                                                aaTAGAATTACCGCAAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCa                                                                                   >  1:1481025/1‑150 (MQ=255)
                                                                     aaTTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCTTTCTGTTAAGAATGTAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGc                                                                              <  1:829655/150‑1 (MQ=255)
                                                                     aaTTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGTAACGCCATTCGc                                                                              <  1:918400/150‑1 (MQ=255)
                                                                        tcccGCCAGCAGCCGCAATCATGTACTGAATAGCGGGGATCTATAGGTTAAGAATGGAAACTGCCGGGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTg                                                                           <  2:321369/148‑1 (MQ=255)
                                                                           cGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGCATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTATGATAGGCTGCCATATGGAAACGCCATTCGCTTGAgg                                                                        >  1:256524/1‑150 (MQ=255)
                                                                           cGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCtt                                                                            >  2:1597120/1‑146 (MQ=255)
                                                                           cGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCtt                                                                            <  1:1597120/146‑1 (MQ=255)
                                                                             ccAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAgggg                                                                      <  1:1586341/150‑1 (MQ=255)
                                                                                gccgccGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGATGTGATATGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAga                                                                   >  2:683621/1‑149 (MQ=255)
                                                                                   gccgcAATCATGTAGTGAATAGCGGGTATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACTTCATGATTGTGAACAACCGCGTTGGCATCCGGTTTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCttt                                                                <  1:72980/150‑1 (MQ=255)
                                                                                   gccgcAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCttt                                                                >  1:1252130/1‑150 (MQ=255)
                                                                                      gcAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTttcct                                                             >  1:1627747/1‑150 (MQ=255)
                                                                                        aaTCATGTATAGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACCGCATGATTGTGAAGAACCGCGTTGGCATACGGTCTGCTTTTATAGGCAGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTatacttc                                                           <  2:612999/150‑1 (MQ=255)
                                                                                        aaTCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTttccttc                                                           <  2:528605/150‑1 (MQ=255)
                                                                                          tCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGGGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTtccttcct                                                         <  2:952409/150‑1 (MQ=255)
                                                                                           cATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTc                                                        <  2:23980/150‑1 (MQ=255)
                                                                                           cATGTACTGAATACAGTGGATCGATCGGTTAAGACTGGTAAATGCCGTGCAATGACCGGCATGATTTTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGTCCTCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTc                                                        <  1:330611/150‑1 (MQ=255)
                                                                                            aTGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCa                                                       <  1:1258543/150‑1 (MQ=255)
                                                                                                aGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGtt                                                   <  1:1327228/150‑1 (MQ=255)
                                                                                                   gAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTAc                                                >  1:604619/1‑150 (MQ=255)
                                                                                                             gatcgatcGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCTTGATTGTGAACAACCGCGTTGGCATCCGGTCTGCGTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTGACCGTTGCCATc                                      >  2:280485/1‑150 (MQ=255)
                                                                                                             gatcgatcGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATc                                      <  1:1571340/150‑1 (MQ=255)
                                                                                                                cgatcgGTTAAGAATGGAAACTGCCGTGCAATTAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAAt                                   <  2:1272125/150‑1 (MQ=255)
                                                                                                                cgatcgGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAAt                                   <  1:144703/150‑1 (MQ=255)
                                                                                                                cgatcgGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAAt                                   <  1:1546395/150‑1 (MQ=255)
                                                                                                                 gatcgGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATa                                  >  1:144099/1‑150 (MQ=255)
                                                                                                                 gatcgGTTAAGAATCGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTGGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCGTTCCTTCCTCATGTTTACCGTTGCCATCAATa                                  <  2:853070/150‑1 (MQ=255)
                                                                                                                    cgGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACACCCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAg                               >  1:759369/1‑150 (MQ=255)
                                                                                                                     ggTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTATACCGTTGCCATCAATAAAGa                              >  1:717109/1‑150 (MQ=255)
                                                                                                                        tAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACaa                           >  1:879814/1‑150 (MQ=255)
                                                                                                                         aaGATTGGAAACTGCCGTGCAATGAACGTCATTTTTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAat                          <  1:1604167/150‑1 (MQ=255)
                                                                                                                           gAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAatat                        >  1:893093/1‑150 (MQ=255)
                                                                                                                           gAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAatat                        <  2:1124073/150‑1 (MQ=255)
                                                                                                                            aaTGGAAACTGCCGGGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATg                       <  1:1432493/150‑1 (MQ=255)
                                                                                                                             aTGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATCGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTGTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGc                      <  1:1609722/150‑1 (MQ=255)
                                                                                                                              acaatcatGCCGTTCATTGCACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGCg                     >  2:456565/8‑150 (MQ=255)
                                                                                                                                gAAACTGCCCTGCAATGCACGGCAAGATTGTGAACACCCGCCTTGGCATCCTGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGCAGCTGTCCTTCCTCATGTTTACCGTTGCCATCAAAAAAGACAAGATGCGAc                   >  2:232187/1‑150 (MQ=255)
                                                                                                                                    cTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCg               <  1:1113315/150‑1 (MQ=255)
                                                                                                                                        cGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATTTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAg           >  2:860811/1‑150 (MQ=255)
                                                                                                                                          tGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAGtt         >  2:1010995/1‑150 (MQ=255)
                                                                                                                                           gCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAGttt        >  1:1384870/1‑150 (MQ=255)
                                                                                                                                            cAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAGtttt       >  1:245574/1‑150 (MQ=255)
                                                                                                                                             aaTGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAGttttt      <  1:359313/150‑1 (MQ=255)
                                                                                                                                              aTGAACGGCATGATTGTGAACAATCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCTCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAGTTTTTc     <  1:40612/150‑1 (MQ=255)
                                                                                                                                              aTGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATATGCTGCCATATGGAAACGCCATTCGCTGGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAGTTTTTc     <  1:1512643/150‑1 (MQ=255)
                                                                                                                                              aTGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAGTTTTTc     <  2:1546557/150‑1 (MQ=255)
                                                                                                                                              aTGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAGTTTTTc     >  1:1142258/1‑150 (MQ=255)
                                                                                                                                                 aaCGGCATGATTGTGAACGACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAGTTTTTCata  >  2:1725736/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
TGAGAGCCAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAGTTTTTCATA  >  NC_000913/2933256‑2933550

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