Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,226,832 T→G 6.5% F78C (TTT→TGT)  ycgJ → PF05666 family protein YcgJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,226,8320TG 6.5% 293.8 / 10.4 92F78C (TTT→TGT) ycgJPF05666 family protein YcgJ
Reads supporting (aligned to +/- strand):  ref base T (45/41);  new base G (2/4);  total (47/45)
Fisher's exact test for biased strand distribution p-value = 4.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

GACTTCAATCTCCGGCGTGCGGCGTTGTTTGCGACCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACA  >  NC_000913/1226685‑1226981
                                                                                                                                                   |                                                                                                                                                     
gACTTCAATCTCCGGCGTGCGGCGTTGTTTGCGACCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTa                                                                                                                                                     <  2:87218/150‑1 (MQ=255)
  cTTCAATCTCCGGCGTGCGGCGTTGTTTGCGACCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACa                                                                                                                                                   >  1:95246/1‑150 (MQ=255)
     cAATCTCCGGCGTGCGGCGTTGTTTGCGACCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCCTTACAAGGCTACGCCCCCAGCGAATTTACATTc                                                                                                                                                >  2:131724/1‑150 (MQ=255)
     cAATCTCCGGCGTGCGGCGTTGTTTGCGACCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTc                                                                                                                                                <  1:653909/150‑1 (MQ=255)
     cAATCTCCGGCGTGCGGCGTTGTTTGCGACCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTc                                                                                                                                                >  2:864549/1‑150 (MQ=255)
       aTCTCCGGCGTGCGGCGTTGTTTGCGACCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATGTAAATTCGc                                                                                                                                              >  1:1477278/1‑150 (MQ=255)
        tctcCGGCGTGCGGCGTTGTTTGCGACCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTATGATCCCAGCGAATTTACATTCGCt                                                                                                                                             <  1:293472/150‑1 (MQ=255)
        tctcCGGCGTGCGGCGTTGTTTGCGACCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATGTAAATTCGCt                                                                                                                                             <  2:1477278/150‑1 (MQ=255)
         ctcCGGCGTGCGGCGTTGTTTGCGACCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTa                                                                                                                                            >  2:797841/1‑150 (MQ=255)
          tcCGGCGTGCGGCGTTGTTTGCGACCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATGTAAATTCGCtgg                                                                                                                                           <  2:1544647/150‑3 (MQ=255)
           ccGGCGTGCGGCGTTGTTTGCGACCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATGTAAATTCGCtggg                                                                                                                                          <  1:550608/150‑4 (MQ=255)
           ccGGCGTGCGGCGTTGTTTGCGACCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATGTAAATTCGCtggg                                                                                                                                          <  1:825312/150‑4 (MQ=255)
             ggCGTGCGGCGTTGTTTGCGACCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACgg                                                                                                                                        >  2:199038/1‑150 (MQ=255)
               cGTGCGGCGTTGTTTGCGACCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGgtg                                                                                                                                      <  2:1671525/150‑1 (MQ=255)
                    ggcgTTGTTTGCGACCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTAttt                                                                                                                                 >  1:77108/1‑150 (MQ=255)
                     gcgTTGTTTGCGACCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGAGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTAtttt                                                                                                                                >  2:1087659/1‑150 (MQ=255)
                       gttgttTGCGACCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGTAAACTTACAATCATTACAAGGCTACGATCCCAGCTAATTTACATTCGCTAACGGTGTATTTTGc                                                                                                                              <  1:95893/150‑1 (MQ=255)
                          gttTGCGACCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGAt                                                                                                                           >  1:1500366/1‑150 (MQ=255)
                            ttGCGACCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGt                                                                                                                         <  1:1911621/150‑1 (MQ=255)
                            ttGCGACCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGt                                                                                                                         >  2:996007/1‑150 (MQ=255)
                             tGCGACCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGtt                                                                                                                        >  1:854757/1‑150 (MQ=255)
                               cGACCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTaa                                                                                                                      <  1:1069224/150‑1 (MQ=255)
                                 aCCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAg                                                                                                                    <  1:864549/150‑1 (MQ=255)
                                 aCCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAg                                                                                                                    >  1:653364/1‑150 (MQ=255)
                                    cGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAAc                                                                                                                                          >  1:183881/1‑125 (MQ=255)
                                    cGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAAc                                                                                                                                          <  2:183881/125‑1 (MQ=255)
                                        tatTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAAtt                                                                                                             >  1:1084292/1‑150 (MQ=255)
                                          tttGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTAt                                                                                                           <  1:1216282/150‑1 (MQ=255)
                                           ttGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATg                                                                                                          <  1:1676438/150‑1 (MQ=255)
                                            tgtgtgAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGt                                                                                                         <  1:1153990/150‑1 (MQ=255)
                                             gtgtgAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTc                                                                                                        >  2:1874657/1‑150 (MQ=255)
                                                  aaCTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGtgat                                                                                                   <  2:1026911/150‑1 (MQ=255)
                                                  aaCTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGtgat                                                                                                   >  2:17893/1‑150 (MQ=255)
                                                   aCTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGtgatg                                                                                                  >  1:1951195/1‑150 (MQ=255)
                                                       aGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCACAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCg                                                                                              <  2:1509627/150‑1 (MQ=255)
                                                           ggCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTAt                                                                                          >  1:1993231/1‑150 (MQ=255)
                                                           ggCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTAt                                                                                          <  2:1084292/150‑1 (MQ=255)
                                                              aTTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTAtttt                                                                                       <  1:1647852/150‑1 (MQ=255)
                                                                  ccccAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGgtg                                                                                   <  2:1902648/150‑1 (MQ=255)
                                                                  ccccAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGgtg                                                                                   <  2:1943672/150‑1 (MQ=255)
                                                                  ccccAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGgtg                                                                                   <  2:1423261/150‑1 (MQ=255)
                                                                   cccAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGgtgt                                                                                  >  1:113802/1‑150 (MQ=255)
                                                                      agagTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTATATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGa                                                                               >  2:1652857/1‑150 (MQ=255)
                                                                       tagTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGTTAACGGTGTATTTTGCGATGTTATAGAAAAATTATGTCGTGATGATCGTTATTTTGGGGTGGAt                                                                              <  2:415100/149‑1 (MQ=255)
                                                                        agTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATg                                                                             >  1:1904375/1‑150 (MQ=255)
                                                                               aaGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTgatgat                                                                                                >  2:1090792/1‑124 (MQ=255)
                                                                               aaGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTgatgat                                                                                                <  1:1090792/124‑1 (MQ=255)
                                                                                aGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTAtttt                                                                                       <  1:1362836/132‑1 (MQ=255)
                                                                                aGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTAtttt                                                                                       >  2:1362836/1‑132 (MQ=255)
                                                                                 ggAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGa                                                                    >  2:926806/1‑150 (MQ=255)
                                                                                    aaTATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTg                                                                 >  2:1356931/1‑150 (MQ=255)
                                                                                     atatCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTgg                                                                <  1:1708766/150‑1 (MQ=255)
                                                                                     atatCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTgg                                                                <  1:797841/150‑1 (MQ=255)
                                                                                       atCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGaa                                                              >  1:619502/1‑150 (MQ=255)
                                                                                            ggCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAaatc                                                         >  2:1802744/1‑150 (MQ=255)
                                                                                            ggCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAaatc                                                         <  2:95246/150‑1 (MQ=255)
                                                                                              cGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAaatcaa                                                       >  1:613359/1‑150 (MQ=255)
                                                                                              cGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAaatcaa                                                       <  1:420872/150‑1 (MQ=255)
                                                                                               gAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAaatcaat                                                      >  2:1783171/1‑150 (MQ=255)
                                                                                                aaaaaGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAaatcaata                                                     >  1:1338502/1‑149 (MQ=255)
                                                                                                     gCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAacc                                                <  2:1500366/150‑1 (MQ=255)
                                                                                                        gCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAaccacc                                             <  2:854757/150‑1 (MQ=255)
                                                                                                          tGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACcaccac                                           <  2:619502/150‑1 (MQ=255)
                                                                                                          tGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACcaccac                                           >  2:384821/1‑150 (MQ=255)
                                                                                                            aaaaCTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACaa                                         >  2:100322/1‑150 (MQ=255)
                                                                                                             aaaCTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACaaa                                        >  2:1276534/1‑150 (MQ=255)
                                                                                                               aCTTACAATCATTACAAGGCTACGATCCCAGCGAATTTATATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGa                                      <  1:1652857/150‑1 (MQ=255)
                                                                                                                cTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGAt                                     >  2:1675351/1‑150 (MQ=255)
                                                                                                                  tACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGt                                   >  2:657771/1‑150 (MQ=255)
                                                                                                                    cAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGtta                                 <  1:1601463/150‑1 (MQ=255)
                                                                                                                    cAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGtta                                 >  2:498748/1‑150 (MQ=255)
                                                                                                                     aaTCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGttat                                <  1:1485035/150‑1 (MQ=255)
                                                                                                                      aTCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGttatt                               >  1:1169951/1‑150 (MQ=255)
                                                                                                                          ttACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTAtg                           <  1:1830577/150‑1 (MQ=255)
                                                                                                                               aGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAaatc                                                         >  2:1585034/1‑115 (MQ=255)
                                                                                                                               aGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAaatc                                                         <  1:1585034/115‑1 (MQ=255)
                                                                                                                               aGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAacca                                               >  1:688408/1‑125 (MQ=255)
                                                                                                                               aGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAacca                                               <  2:688408/125‑1 (MQ=255)
                                                                                                                               aGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCg                      >  1:542206/1‑150 (MQ=255)
                                                                                                                                ggCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATc                                                                                               <  1:1212116/76‑1 (MQ=255)
                                                                                                                                ggCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATc                                                                                               >  2:1212116/1‑76 (MQ=255)
                                                                                                                                 gCTACGATCCCATCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCgtg                    <  1:901665/150‑1 (MQ=255)
                                                                                                                                   tACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCgtgag                  >  1:512122/1‑150 (MQ=255)
                                                                                                                                     caccgttaGCGAATGTAAATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCgtgagtg                >  2:1624522/8‑150 (MQ=255)
                                                                                                                                     cGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACcaca                                          >  1:538747/1‑124 (MQ=255)
                                                                                                                                     cGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACcaca                                          <  2:538747/124‑1 (MQ=255)
                                                                                                                                        tCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCGTGAGTGATg             <  2:653364/150‑1 (MQ=255)
                                                                                                                                          ccAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGta           >  2:575394/1‑149 (MQ=255)
                                                                                                                                          ccAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGTc           >  2:1384982/1‑150 (MQ=255)
                                                                                                                                           cAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGTCt          >  2:718512/1‑150 (MQ=255)
                                                                                                                                              cGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGt       <  1:100322/150‑1 (MQ=255)
                                                                                                                                               gAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGtt      <  2:113802/150‑1 (MQ=255)
                                                                                                                                                   ttAGATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTCTGTGTCGTGAGTGATGTCTTGTTTACa  <  1:58718/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GACTTCAATCTCCGGCGTGCGGCGTTGTTTGCGACCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACA  >  NC_000913/1226685‑1226981

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