Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,771,180 G→T 7.4% L11L (CTG→CTT)  yfjT → CP4‑57 prophage; uncharacterized protein YfjT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,771,1800GT 7.4% 294.3 / 8.3 92L11L (CTG→CTT) yfjTCP4‑57 prophage; uncharacterized protein YfjT
Reads supporting (aligned to +/- strand):  ref base G (41/44);  new base T (3/4);  total (44/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.10e-01

CGTATACGAAAGAAAAATACTTCGAATATTCCCATGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGG  >  NC_000913/2771033‑2771329
                                                                                                                                                   |                                                                                                                                                     
cGTATACGAAAGAAAAATACTTCGAATATTCCCATGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGc                                                                                                                                                     <  2:1617470/150‑1 (MQ=255)
 gTATACGAAAGAAAAATACTTCGAATATTCCCATGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCa                                                                                                                                                    <  2:499726/150‑1 (MQ=255)
   ataCGAAAGAAAAATACTTCGAATATTCCCATGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAg                                                                                                                                                  <  1:926902/150‑1 (MQ=255)
      cGAAAGAAAAATACTTCGAATATTCCCATGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTAc                                                                                                                                               >  2:179706/1‑150 (MQ=255)
        aaagaaaAATACTTCGAATATTCCCATGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACaa                                                                                                                                             >  2:312238/1‑150 (MQ=255)
           gaaaAATACTTCGAATATTCCCATGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTTGCCAGTACaaacc                                                                                                                                          <  2:311713/150‑4 (MQ=255)
            aaaaaTACTTCGAATATTCCCATGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTTGCCAGTACaaaccg                                                                                                                                         <  1:746964/150‑5 (MQ=255)
            aaaaaTACTTCGAATATTCCCATGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTTGCCAGTACaaaccg                                                                                                                                         <  2:49027/150‑5 (MQ=255)
            aaaaaTACTTCGAATATTCCCATGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGtt                                                                                                                                         >  2:610609/1‑150 (MQ=255)
             aaaaTACTTCGAATATTCCCATGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTTGCCAGTACaaaccgg                                                                                                                                        <  1:803547/150‑6 (MQ=255)
              aaaTACTTCGAATATTCCCATGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCg                                                                                                                                       <  1:859835/150‑1 (MQ=255)
              aaaTACTTCGAATATTCCCATGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCg                                                                                                                                       >  2:1058012/1‑150 (MQ=255)
                aTACTTCGAATATTCCCATGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGcc                                                                                                                                     <  1:94808/150‑1 (MQ=255)
                 tACTTCGAATAATCCCATGACGGCAAACAGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGCGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCa                                                                                                                                    >  2:1124162/1‑150 (MQ=255)
                   cTTCGAATATTCCCATGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGc                                                                                                                                  >  2:1009601/1‑150 (MQ=255)
                   cTTCGAATATTCCCATGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGc                                                                                                                                  <  1:435287/150‑1 (MQ=255)
                   cTTCGAATATTCCCATGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGc                                                                                                                                  >  2:1066731/1‑150 (MQ=255)
                      cGAATATTCCCATGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCttt                                                                                                                               <  1:288323/150‑1 (MQ=255)
                       gAATATTCCCATGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTa                                                                                                                              >  1:1392542/1‑150 (MQ=255)
                         atatTCCCATGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACa                                                                                                                            <  2:1531511/150‑1 (MQ=255)
                              cccATGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCctc                                                                                                                       >  2:248573/1‑150 (MQ=255)
                                 aTGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGc                                                                                                                    <  1:975291/150‑1 (MQ=255)
                                 aTGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAAACGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGc                                                                                                                    <  2:1243269/150‑1 (MQ=255)
                                  tGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCg                                                                                                                   <  1:550320/150‑1 (MQ=255)
                                     cGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGAAATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTc                                                                                                                >  2:1331661/1‑150 (MQ=255)
                                        cAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTcccc                                                                                                             <  1:248573/150‑1 (MQ=255)
                                        cAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTcccc                                                                                                             >  1:413672/1‑150 (MQ=255)
                                        cAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTcccc                                                                                                             <  1:179706/150‑1 (MQ=255)
                                          aaCTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAAGGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCgg                                                                                                           <  1:309196/150‑1 (MQ=255)
                                          aaCTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCgg                                                                                                           <  1:1062223/150‑1 (MQ=255)
                                            cTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGAc                                                                                                         <  2:1473446/150‑1 (MQ=255)
                                                tGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTg                                                                                                     >  1:1479183/1‑150 (MQ=255)
                                                tGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTg                                                                                                     >  1:1473374/1‑150 (MQ=255)
                                                tGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTg                                                                                                     >  2:659352/1‑150 (MQ=255)
                                                tGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTg                                                                                                     >  1:1041746/1‑150 (MQ=255)
                                                tGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTg                                                                                                     <  1:1009601/150‑1 (MQ=255)
                                                  ggACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGc                                                                                                   <  2:561449/150‑1 (MQ=255)
                                                              aTACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACg                                                                                       >  1:1443782/1‑150 (MQ=255)
                                                              aTACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACg                                                                                       <  1:1352800/150‑1 (MQ=255)
                                                               tACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACgg                                                                                      <  2:1278797/150‑1 (MQ=255)
                                                               tACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACgg                                                                                      <  1:543328/150‑1 (MQ=255)
                                                                 ccGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTc                                                                                    >  1:1036545/1‑150 (MQ=255)
                                                                 ccGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTc                                                                                    >  1:174706/1‑150 (MQ=255)
                                                                 ccGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTc                                                                                    >  1:537141/1‑150 (MQ=255)
                                                                      aaaCCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGc                                                                               >  2:1358144/1‑150 (MQ=255)
                                                                        aCCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCC‑GACTCTGGCAACAGGGGTACGGTCAGGGCAAt                                                                            >  1:452128/1‑150 (MQ=255)
                                                                            tGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCg                                                                         <  1:659352/150‑1 (MQ=255)
                                                                             gATACCCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCgg                                                                        >  2:602805/1‑150 (MQ=255)
                                                                               tAACCGGGAATAGAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGaa                                                                      <  2:291703/150‑1 (MQ=255)
                                                                                 aCCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATa                                                                    >  1:230843/1‑150 (MQ=255)
                                                                                    gggAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACaa                                                                 >  1:726768/1‑149 (MQ=255)
                                                                                    gggAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAg                                                                 <  1:759146/150‑1 (MQ=255)
                                                                                     ggAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAgt                                                                <  1:1440351/150‑1 (MQ=255)
                                                                                       aaTAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAgtgt                                                              <  2:1345862/150‑1 (MQ=255)
                                                                                           aaaaCAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACg                                                          >  2:633191/1‑150 (MQ=255)
                                                                                           aaaaCAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACg                                                          >  1:770069/1‑150 (MQ=255)
                                                                                           aaaaCAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACg                                                          >  2:1559958/1‑150 (MQ=255)
                                                                                             aacaacGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGa                                                        <  1:1358144/150‑1 (MQ=255)
                                                                                                 acGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGcc                                                    <  1:312238/150‑1 (MQ=255)
                                                                                                  cGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGc                                                                                                   <  2:1275291/102‑1 (MQ=255)
                                                                                                  cGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGc                                                                                                   >  1:1275291/1‑102 (MQ=255)
                                                                                                    agaAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGc                                                 >  1:17810/1‑150 (MQ=255)
                                                                                                     gaAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCg                                                <  1:1124162/150‑1 (MQ=255)
                                                                                                     gaAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCg                                                >  2:54097/1‑150 (MQ=255)
                                                                                                        aGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGaa                                             <  1:1058012/150‑1 (MQ=255)
                                                                                                          gagaCATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGTTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGTCAATTCGGAATACAGTGTGACGTATGCCAGCGGAAAg                                           <  2:252746/150‑1 (MQ=255)
                                                                                                             aCATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATg                                        <  2:1036545/150‑1 (MQ=255)
                                                                                                                 aaaTGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTa                                    >  2:1238223/1‑150 (MQ=255)
                                                                                                                  aaTGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTAc                                   <  1:22224/150‑1 (MQ=255)
                                                                                                                  aaTGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTAc                                   >  1:756304/1‑150 (MQ=255)
                                                                                                                       aaaTACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATcg                              <  2:468537/150‑2 (MQ=255)
                                                                                                                           aCGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTg                          <  2:1443782/150‑1 (MQ=255)
                                                                                                                           aCGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTg                          <  1:1269702/150‑1 (MQ=255)
                                                                                                                             gTTCACTAAGCCGGTTTGTACTGGAAAGTACAATGTTCGCCAGCGTTACAGTCTCTGCGGTCCCCGGACTCTTTCAACAGGGGTACGGTCCGGGCAATGCTTAATACAGGGGGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCa                        <  2:1380701/150‑1 (MQ=255)
                                                                                                                              ttCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCAc                       <  1:619940/150‑1 (MQ=255)
                                                                                                                               tCACTAAGCCGGTTTGTAATGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACg                      <  2:1046651/150‑1 (MQ=255)
                                                                                                                                 aCTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACggg                    >  1:1158259/1‑150 (MQ=255)
                                                                                                                                 aCTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCACCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACggg                    >  2:1572808/1‑150 (MQ=255)
                                                                                                                                  cTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGa                   >  1:433970/1‑150 (MQ=255)
                                                                                                                                     gcgaacattGTACTTGCCAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAAt                >  2:4743/8‑150 (MQ=255)
                                                                                                                                      gCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATc               >  1:332453/1‑150 (MQ=255)
                                                                                                                                      cgaacattGTACTTGCCAGTACAATGTTCGCCAGGTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATc               >  1:684631/7‑150 (MQ=255)
                                                                                                                                       ccGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATcc              >  1:762596/1‑150 (MQ=255)
                                                                                                                                         ggTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCgg            <  2:174706/150‑1 (MQ=255)
                                                                                                                                         acattGTACTTGCCAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGCCTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCgg            >  2:1593695/4‑150 (MQ=255)
                                                                                                                                           tttGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGAc          <  2:1479183/150‑1 (MQ=255)
                                                                                                                                           tttGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGAc          <  2:1473374/150‑1 (MQ=255)
                                                                                                                                             tGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCa        >  2:972918/1‑150 (MQ=255)
                                                                                                                                              gTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCTGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGTGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAg       <  2:556062/150‑1 (MQ=255)
                                                                                                                                              gTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACAATCAACTGCACGGGAAATCCGGACCac       >  2:554122/1‑149 (MQ=255)
                                                                                                                                                 cTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAAc    >  2:138823/1‑150 (MQ=255)
                                                                                                                                                 cTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAAc    >  2:1224488/1‑150 (MQ=255)
                                                                                                                                                   ggCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACgg  <  1:633191/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CGTATACGAAAGAAAAATACTTCGAATATTCCCATGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGG  >  NC_000913/2771033‑2771329

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