Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,771,180 G→T 7.9% 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.9% 310.3 / 10.3 95L11L (CTG→CTTyfjTCP4‑57 prophage; uncharacterized protein YfjT
Reads supporting (aligned to +/- strand):  ref base G (44/43);  new base T (2/6);  total (46/49)
Fisher's exact test for biased strand distribution p-value = 2.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.42e-01

CCCGTATACGAAAGAAAAATACTTCGAATATTCCCATGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGG  >  NC_000913/2771031‑2771329
                                                                                                                                                     |                                                                                                                                                     
cccGTATACGAAAGAAAAATACTTCGAATATTCCCATGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTg                                                                                                                                                       >  2:1530555/1‑150 (MQ=255)
        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                                                                                                                                          <  2:896914/131‑4 (MQ=255)
                                cccATGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTTGCCAGTACaaacc                                                                                                                                          >  1:896914/1‑128 (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:1527045/150‑1 (MQ=255)
                                                   gggACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTgg                                                                                                    >  1:766391/1‑150 (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                                                                                                                         >  1:50568/1‑122 (MQ=255)
                                                          ctgctgATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGcc                                                                                                                         <  2:50568/122‑1 (MQ=255)
                                                                aTACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACg                                                                                       >  2:576048/1‑150 (MQ=255)
                                                                 tACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACgg                                                                                      <  2:228590/150‑1 (MQ=255)
                                                                 tACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACgg                                                                                      <  2:1079953/150‑1 (MQ=255)
                                                                   ccGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTc                                                                                    >  1:184766/1‑150 (MQ=255)
                                                                    cGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCa                                                                                   >  2:43444/1‑150 (MQ=255)
                                                                     gTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTg                                                                                                     <  2:936340/131‑1 (MQ=255)
                                                                     gTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTg                                                                                                     >  1:936340/1‑131 (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:723417/1‑150 (MQ=255)
                                                                                                                                           ggTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCgg            <  1:897548/150‑1 (MQ=255)
                                                                                                                                           ggTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCgg            <  2:28675/150‑1 (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)
                                                                                                                                                     |                                                                                                                                                     
CCCGTATACGAAAGAAAAATACTTCGAATATTCCCATGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGG  >  NC_000913/2771031‑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: