Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,771,183 A→C 7.0% A12A (GCA→GCC)  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,1830AC 7.0% 310.7 / 13.4 100A12A (GCA→GCC) yfjTCP4‑57 prophage; uncharacterized protein YfjT
Reads supporting (aligned to +/- strand):  ref base A (46/47);  new base C (2/5);  total (48/52)
Fisher's exact test for biased strand distribution p-value = 4.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

CGAAAGAAAAATACTTCGAATATTCCCATGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTA  >  NC_000913/2771039‑2771331
                                                                                                                                                |                                                                                                                                                    
cGAAAGAAAAATACTTCGAATATTCCCATGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTAc                                                                                                                                                 >  1:588246/1‑150 (MQ=255)
cGAAAGAAAAATACTTCGAATATTCCCATGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTAc                                                                                                                                                 <  2:328185/150‑1 (MQ=255)
cGAAAGAAAAATACTTCGAATATTCCCATGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTAc                                                                                                                                                 >  1:180961/1‑150 (MQ=255)
   aagaaaAATACTTCGAATATTCCCATGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAAt                                                                                                                                              <  2:273817/150‑1 (MQ=255)
    agaaaAATACTTCGAATATTCCCATGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTTGCCAGTACaaac                                                                                                                                             <  1:394471/150‑3 (MQ=255)
     gaaaAATACTTCGAATATTCCCATGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGt                                                                                                                                            >  2:897548/1‑150 (MQ=255)
      aaaaaTACTTCGAATATTCCCATGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTTGCCAGTACaaaccg                                                                                                                                           <  2:1071795/150‑5 (MQ=255)
      aaaaaTACTTCGAATATTCCCATGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTTGCCAGTACaaaccg                                                                                                                                           <  2:1331896/150‑5 (MQ=255)
      aaaaaTACTTCGAATATTCCCATGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGtt                                                                                                                                           <  1:1414693/150‑1 (MQ=255)
       aaaaTACTTCGAATATTCCCATGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTTGCCAGTACaaaccgg                                                                                                                                          <  2:1124602/150‑6 (MQ=255)
        aaaTACTTCGAATATTCCCATGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCg                                                                                                                                         <  2:785459/150‑1 (MQ=255)
             catCGAATATTCCCATGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAACCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGc                                                                                                                                    <  2:694904/148‑1 (MQ=255)
             cTTCGAATATTCCCATGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGc                                                                                                                                    >  1:1274783/1‑150 (MQ=255)
              ttCGAATATTCCCATGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCt                                                                                                                                   >  2:919715/1‑150 (MQ=255)
                      ttCCCATGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGcc                                                                                                                           <  1:377685/150‑1 (MQ=255)
                        cccATGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTTGCCAGTACaaacc                                                                                                                                            >  1:896914/1‑128 (MQ=255)
                        cccATGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTTGCCAGTACaaacc                                                                                                                                            <  2:896914/131‑4 (MQ=255)
                        cccATGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCctc                                                                                                                         >  2:1437327/1‑150 (MQ=255)
                          cATGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTg                                                                                                                       <  1:802654/150‑1 (MQ=255)
                           aTGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGc                                                                                                                      <  2:1381617/150‑1 (MQ=255)
                            tGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCg                                                                                                                     >  2:43830/1‑150 (MQ=255)
                             gACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCgg                                                                                                                    <  1:689912/150‑1 (MQ=255)
                             gACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCgg                                                                                                                    <  1:847679/150‑1 (MQ=255)
                               cGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTc                                                                                                                  >  2:76596/1‑150 (MQ=255)
                               cGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCAATAAGCCGGTTTGTACTGGCAAGTACAATGTTCGACAGCTTTACAGCCTCTGCGGTc                                                                                                                  >  2:735715/1‑150 (MQ=255)
                                  cAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTcccc                                                                                                               >  2:757607/1‑150 (MQ=255)
                                        gATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGActc                                                                                                         >  2:201321/1‑150 (MQ=255)
                                          tGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTg                                                                                                       >  2:251751/1‑150 (MQ=255)
                                           gggACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATGAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTgg                                                                                                      <  1:683605/150‑1 (MQ=255)
                                           gggACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTgg                                                                                                      >  1:766391/1‑150 (MQ=255)
                                           gggACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTgg                                                                                                      <  1:1527045/150‑1 (MQ=255)
                                             gACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCa                                                                                                    <  1:438524/150‑1 (MQ=255)
                                             gACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCa                                                                                                    <  1:622776/150‑1 (MQ=255)
                                             gACGGCTGCTGATACCCTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACAAAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCa                                                                                                    >  2:470425/1‑150 (MQ=255)
                                               cGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAAc                                                                                                  <  1:1140145/150‑1 (MQ=255)
                                                  ctgctgATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGcc                                                                                                                           <  2:50568/122‑1 (MQ=255)
                                                  ctgctgATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGcc                                                                                                                           >  1:50568/1‑122 (MQ=255)
                                                        aTACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACg                                                                                         >  2:576048/1‑150 (MQ=255)
                                                         tACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACgg                                                                                        <  2:1079953/150‑1 (MQ=255)
                                                         tACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACgg                                                                                        <  2:228590/150‑1 (MQ=255)
                                                           ccGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTc                                                                                      >  1:184766/1‑150 (MQ=255)
                                                            cGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCa                                                                                     >  2:43444/1‑150 (MQ=255)
                                                             gTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTg                                                                                                       >  1:936340/1‑131 (MQ=255)
                                                             gTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTg                                                                                                       <  2:936340/131‑1 (MQ=255)
                                                             gTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAg                                                                                    >  1:829586/1‑150 (MQ=255)
                                                                aaaCCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGc                                                                                 <  2:1178377/150‑1 (MQ=255)
                                                                    ccTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATg                                                                             >  1:760216/1‑150 (MQ=255)
                                                                        aTAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCAATAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGTGGa                                                                         <  1:735715/150‑1 (MQ=255)
                                                                         tAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTAATGGCAAGTACAATGTTCGCCAGCATTCCAGCATCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGACAGGGCAATGCGGaa                                                                        >  2:72295/1‑150 (MQ=255)
                                                                           aCCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCg                                                                                                              >  1:635621/1‑110 (MQ=255)
                                                                           aCCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCg                                                                                                              <  2:635621/110‑1 (MQ=255)
                                                                           aCCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATa                                                                      <  1:251751/150‑1 (MQ=255)
                                                                           aCCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATa                                                                      >  2:991599/1‑150 (MQ=255)
                                                                            ccGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATAc                                                                     >  2:316411/1‑150 (MQ=255)
                                                                              gggAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAg                                                                   <  1:1046521/150‑1 (MQ=255)
                                                                                gAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCTCTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAgtg                                                                 <  1:974316/150‑1 (MQ=255)
                                                                                 aaTAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAgtgt                                                                <  2:812871/150‑1 (MQ=255)
                                                                                 aaTAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGGCAGGGCAATGCGGAATACAgtgt                                                                <  1:1420650/150‑1 (MQ=255)
                                                                                   tAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGa                                                              >  2:1104893/1‑150 (MQ=255)
                                                                                     aaaaCAACGAGAAAGGAGACATAAATGATAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAgtg                                                                 <  1:1161012/145‑1 (MQ=255)
                                                                                     aaaaCAACGAGAAAGGAGACATAAATGATAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAgtg                                                                 >  2:1161012/1‑145 (MQ=255)
                                                                                     aaaaCAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGCTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACg                                                            >  2:991671/1‑150 (MQ=255)
                                                                                      aaaCAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACgg                                                           >  1:186823/1‑150 (MQ=255)
                                                                                         caacGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATg                                                        >  1:971219/1‑150 (MQ=255)
                                                                                              agaAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGc                                                   >  1:941960/1‑150 (MQ=255)
                                                                                               tataGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGATTACAGTGTGACGGATGCCAGCg                                                  <  2:1515215/147‑1 (MQ=255)
                                                                                               gaAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCg                                                  <  2:1274783/150‑1 (MQ=255)
                                                                                                   ggAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCa                                                                                                    >  2:469772/1‑96 (MQ=255)
                                                                                                   ggAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCa                                                                                                    <  1:469772/96‑1 (MQ=255)
                                                                                                        cATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGt                                         <  1:874846/150‑1 (MQ=255)
                                                                                                         aTAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGtt                                        <  1:1437327/150‑1 (MQ=255)
                                                                                                           aaaTGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTa                                      <  2:766391/150‑1 (MQ=255)
                                                                                                            aaTGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTAc                                     <  1:1530555/150‑1 (MQ=255)
                                                                                                             aTGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTAcc                                    <  1:186406/150‑1 (MQ=255)
                                                                                                             aTGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTAcc                                    <  1:76596/150‑1 (MQ=255)
                                                                                                                aaaaTACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATc                                 <  2:145434/150‑1 (MQ=255)
                                                                                                                 aaaTACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGTACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCa                                <  2:588246/150‑1 (MQ=255)
                                                                                                                    tACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTAGCAACGGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACt                             >  2:412676/1‑150 (MQ=255)
                                                                                                                      cGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGc                           >  1:274465/1‑150 (MQ=255)
                                                                                                                      cGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGc                           >  2:568418/1‑150 (MQ=255)
                                                                                                                        ttCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCAc                         <  1:316411/150‑1 (MQ=255)
                                                                                                                         tCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACg                        <  1:43830/150‑1 (MQ=255)
                                                                                                                          cccTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATAGAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACgg                       >  2:482277/3‑150 (MQ=255)
                                                                                                                           aCTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACggg                      >  1:192634/1‑150 (MQ=255)
                                                                                                                           aCTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACggg                      >  1:152958/1‑150 (MQ=255)
                                                                                                                            cTAAGCCGGTTTGTACTTGGAAGTACAATGTTCGCGAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGa                     <  2:1255128/150‑1 (MQ=255)
                                                                                                                                  aacattGTACTTGCCAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGGTGCCAGCGGAAAGATGTTTACCATCAACTGCACGGTAAATCCg               >  1:800456/5‑150 (MQ=255)
                                                                                                                                   ggTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCgg              <  2:28675/150‑1 (MQ=255)
                                                                                                                                   ggTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCgg              <  1:897548/150‑1 (MQ=255)
                                                                                                                                   ggTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCgg              >  2:723417/1‑150 (MQ=255)
                                                                                                                                         tACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCTGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGa        >  1:126909/1‑150 (MQ=255)
                                                                                                                                          aCTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGaa       >  1:450892/1‑150 (MQ=255)
                                                                                                                                           cTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAAc      >  2:568310/1‑150 (MQ=255)
                                                                                                                                           cTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAAc      >  2:563018/1‑150 (MQ=255)
                                                                                                                                           cTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGAAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCAAGGGAAATCCGGACCAGAAc      >  2:246274/1‑150 (MQ=255)
                                                                                                                                             ggCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACgg    >  1:528927/1‑150 (MQ=255)
                                                                                                                                             ggCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACgg    <  2:186823/150‑1 (MQ=255)
                                                                                                                                               cAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTa  >  1:1239269/1‑150 (MQ=255)
                                                                                                                                               cAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGGACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTa  <  1:1040992/150‑1 (MQ=255)
                                                                                                                                               cAAGTACAATATTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTa  <  1:911558/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
CGAAAGAAAAATACTTCGAATATTCCCATGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTA  >  NC_000913/2771039‑2771331

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: ‑