Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,226,832 T→G 10.7% 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,8320TG10.7% 257.2 / 22.7 94F78C (TTT→TGT) ycgJPF05666 family protein YcgJ
Reads supporting (aligned to +/- strand):  ref base T (42/42);  new base G (5/5);  total (47/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

CGGACTTCAATCTCCGGCGTGCGGCGTTGTTTGCGACCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTAC  >  NC_000913/1226683‑1226980
                                                                                                                                                     |                                                                                                                                                    
cGGACTTCAATCTCCGGCGTGCGGCGTTGTTTGCGACCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAAtt                                                                                                                                                      >  1:134294/1‑150 (MQ=255)
    cTTCAATCTCCGGCGTGCGGCGTTGTTTGCGACCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACa                                                                                                                                                  <  2:256953/150‑1 (MQ=255)
    cTTCAATCTCCGGCGTGCGGCGTTGTTTGCGACCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACa                                                                                                                                                  >  1:548222/1‑150 (MQ=255)
    cTTCAATCTCCGGCGTGCGGCGTTGTTTGCGACCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACa                                                                                                                                                  >  2:1349341/1‑150 (MQ=255)
    cTTCAATCTCCGGCGTGCGGCGTTGTTTGCGACCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACa                                                                                                                                                  >  1:1499432/1‑150 (MQ=255)
      tCAATCTCCGGCGTGCGGCGTTGTTTGCGACCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACAtt                                                                                                                                                <  1:1397797/150‑1 (MQ=255)
       cAATCTCCGGCGTGCGGCGTTGTTTGCGACCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTc                                                                                                                                               >  2:197117/1‑150 (MQ=255)
          tctcCGGCGTGCGGCGTTGTTTGCGACCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCt                                                                                                                                            <  2:424766/150‑1 (MQ=255)
           ctcCGGCGTGCGGCGTTGTTTGCGACCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATGTAAATTCGCtg                                                                                                                                           <  1:1113202/150‑2 (MQ=255)
           ctcCGGCGTGCGGCGTTGTTTGCGACCCGTATATTTGTGAGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTa                                                                                                                                           >  2:327113/1‑150 (MQ=255)
             ccGGCGTGCGGCGTTGTTTGCGACCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAAc                                                                                                                                         <  1:27994/150‑1 (MQ=255)
             ccGGCGTGCGGCGTTGTTTGCGACCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATGTAAATTCGCtggg                                                                                                                                         <  2:929019/150‑4 (MQ=255)
              cGGCGTGCGGCGTTGTTTGCGACCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATGTAAATTCGCtggga                                                                                                                                        <  2:737353/150‑5 (MQ=255)
                gCGTGCGGCGTTGTTTGCGACCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATGTAAATTCGCtgggatc                                                                                                                                      <  1:103319/150‑7 (MQ=255)
                 cGTGCGGCGTTGTTTGCGACCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGgtg                                                                                                                                     <  2:134294/150‑1 (MQ=255)
                   tGCGGCGTTGTTTGCGACCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTa                                                                                                                                   >  1:1446786/1‑150 (MQ=255)
                      ggcgTTGTTTGCGACCCGTAAATTTGTTTGAACTCAGATGGCATTTCCCCACAGTTAACAACGAAATATCACGGCCAAAAAGCCGCTGAAAACTTAAAATCATTACAAGGCTACGATCCCAGCGAATTTACAATACCTAACGGTGTAatt                                                                                                                                >  2:1396913/1‑150 (MQ=255)
                        cgTTGTTTGCGACCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGGGTATTTTg                                                                                                                              <  2:860922/150‑1 (MQ=255)
                         gttgttTGCGACCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAGAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGc                                                                                                                             >  2:335586/1‑150 (MQ=255)
                         gttgttTGCGACCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGc                                                                                                                             <  1:1165167/150‑1 (MQ=255)
                         gttgttTGCGACCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGc                                                                                                                             <  1:985450/150‑1 (MQ=255)
                         gttgttTGCGACCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGc                                                                                                                             >  2:930915/1‑150 (MQ=255)
                           tgttTGCGACCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGa                                                                                                                           >  2:76769/1‑150 (MQ=255)
                              ttGCGACCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGt                                                                                                                        <  2:198684/150‑1 (MQ=255)
                               tGCGACCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGtt                                                                                                                       <  1:931187/150‑1 (MQ=255)
                                gCGACCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTa                                                                                                                      >  1:480249/1‑150 (MQ=255)
                                   aCCCGTATATTTGTGTGTACTCAGATGGCATTTCCGCAGAGTTAACAATGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACTGTGTATTTTGCGATGTTAAAg                                                                                                                   <  2:338061/150‑1 (MQ=255)
                                      cGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTaa                                                                                                                     <  2:1183473/145‑1 (MQ=255)
                                      cGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTaa                                                                                                                     >  1:1183473/1‑145 (MQ=255)
                                        tataTTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGaaaaa                                                                                                              <  2:1510068/150‑1 (MQ=255)
                                        tataTTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGaaaaa                                                                                                              <  2:224875/150‑1 (MQ=255)
                                          tatTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAAtt                                                                                                            <  2:1052794/150‑1 (MQ=255)
                                               gtgtgAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTc                                                                                                       <  2:577357/150‑1 (MQ=255)
                                                  tgAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGtg                                                                                                    <  2:1174240/150‑1 (MQ=255)
                                                  tgAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGtg                                                                                                    <  1:998498/150‑1 (MQ=255)
                                                  tgAACTAAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAAATTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGtg                                                                                                    >  2:149390/1‑150 (MQ=255)
                                                      cTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGtgatga                                                                                                >  2:160830/1‑150 (MQ=255)
                                                           aTGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTaaa                                                                                                                    >  1:61545/1‑125 (MQ=255)
                                                           aTGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTaaa                                                                                                                    <  2:61545/125‑1 (MQ=255)
                                                              gCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTAtt                                                                                        >  1:1406660/1‑150 (MQ=255)
                                                               cATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTa                                                                                          <  1:622173/147‑1 (MQ=255)
                                                               cATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTa                                                                                          >  2:622173/1‑147 (MQ=255)
                                                                aTTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTAtttt                                                                                      <  2:1057160/150‑1 (MQ=255)
                                                                  ttCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTgg                                                                                    <  1:1396913/150‑1 (MQ=255)
                                                                    ccccAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGgtg                                                                                  >  2:875096/1‑150 (MQ=255)
                                                                     cccAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGgtgt                                                                                 >  1:1173327/1‑150 (MQ=255)
                                                                      ccAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGtgtg                                                                                >  1:631549/1‑150 (MQ=255)
                                                                         gagTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGAt                                                                             <  2:255573/150‑1 (MQ=255)
                                                                          agTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATg                                                                            <  1:841364/150‑1 (MQ=255)
                                                                           gTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATgg                                                                           <  1:509032/150‑1 (MQ=255)
                                                                           gTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATgg                                                                           <  1:1416627/150‑1 (MQ=255)
                                                                           gTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATgg                                                                           <  2:1007786/150‑1 (MQ=255)
                                                                            ttAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGt                                                                          <  2:548222/150‑1 (MQ=255)
                                                                               aCAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAg                                                                       >  1:696399/1‑150 (MQ=255)
                                                                                  aGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACTGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCgg                                                                    <  2:1499432/150‑1 (MQ=255)
                                                                                    gAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAg                                                                  <  2:1043669/150‑1 (MQ=255)
                                                                                        tatCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGa                                                              >  1:1505274/1‑150 (MQ=255)
                                                                                         atCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGaa                                                             >  2:680315/1‑150 (MQ=255)
                                                                                         atCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGaa                                                             >  2:796058/1‑150 (MQ=255)
                                                                                          tctcGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGaaa                                                            >  1:829881/1‑150 (MQ=255)
                                                                                              ggCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAaatc                                                        <  1:891698/150‑1 (MQ=255)
                                                                                                  aagaaGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAaatcaatc                                                    <  1:335586/147‑1 (MQ=255)
                                                                                                  aaaaaGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGAGGGTAAGCGGAGTGGAAAaatcaatc                                                    <  2:944651/150‑1 (MQ=255)
                                                                                                          gCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAaccacc                                            >  2:1356435/1‑150 (MQ=255)
                                                                                                                  cTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGAt                                    >  1:1083475/1‑150 (MQ=255)
                                                                                                                    tACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAATATGt                                  >  1:542568/1‑150 (MQ=255)
                                                                                                                    tACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGt                                  >  1:244343/1‑150 (MQ=255)
                                                                                                                       aaTCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGttat                               <  1:174799/150‑1 (MQ=255)
                                                                                                                        aTCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGttatt                              >  2:496279/1‑150 (MQ=255)
                                                                                                                        aTCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGttatt                              >  2:309599/1‑150 (MQ=255)
                                                                                                                        aTCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGttatt                              >  1:810058/1‑150 (MQ=255)
                                                                                                                          cATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGttattta                            <  2:626286/150‑1 (MQ=255)
                                                                                                                          cATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGttattta                            <  1:357425/150‑1 (MQ=255)
                                                                                                                           aTTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGttatttat                           >  2:1468161/1‑150 (MQ=255)
                                                                                                                           aTTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGttatttat                           <  2:1446786/150‑1 (MQ=255)
                                                                                                                            ttACAAGGCTACGATCCCAGCGAATTTACATTGGCTTACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAATATCAATCAAACCACCACAAAGATGTTATTTAtg                          <  1:882584/150‑1 (MQ=255)
                                                                                                                                aaGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTc                      >  1:130010/1‑150 (MQ=255)
                                                                                                                                 aGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCg                     >  1:1412807/1‑150 (MQ=255)
                                                                                                                                  ggCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCgt                    >  1:1030749/1‑150 (MQ=255)
                                                                                                                                   gCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCgtg                   >  1:354733/1‑150 (MQ=255)
                                                                                                                                    cTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCGATCAAACCACCACAAAGATGTTATTTATGTGTCgtgt                  <  1:1291480/150‑2 (MQ=255)
                                                                                                                                    cTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCgtga                  >  1:1219575/1‑150 (MQ=255)
                                                                                                                                    cTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCgtga                  >  2:722663/1‑150 (MQ=255)
                                                                                                                                        gATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAACAAAAATTATGTCGTGATGATCGTTATATTGTTGTGGATGGTAAGCTGAGAGGAAAAAACAATCAAACAACAAAAAATATTTAATTTTTTGGACgtgagtca              >  1:502518/1‑148 (MQ=255)
                                                                                                                                           gttaGCGAATGTAAATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGt           >  1:1197070/4‑150 (MQ=255)
                                                                                                                                           gttaGCGAATGTAAATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGt           >  1:170229/4‑150 (MQ=255)
                                                                                                                                           gttaGCGAATGTAAATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGt           >  2:318558/4‑150 (MQ=255)
                                                                                                                                           cccAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGt           <  1:327113/150‑1 (MQ=255)
                                                                                                                                           cccAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGt           >  2:854903/1‑150 (MQ=255)
                                                                                                                                             taGCGAATGTAAATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGTCt         <  2:1197070/149‑1 (MQ=255)
                                                                                                                                             cAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGTCt         >  1:887471/1‑150 (MQ=255)
                                                                                                                                              aGCGAATGTAAATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGTCtt        >  2:1133886/1‑150 (MQ=255)
                                                                                                                                              aGCGAATGTAAATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGTCtt        >  2:777771/1‑150 (MQ=255)
                                                                                                                                                  aaTTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGttt    <  2:342832/150‑1 (MQ=255)
                                                                                                                                                  aaTTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGttt    <  2:1401941/150‑1 (MQ=255)
                                                                                                                                                    tttACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTAc  <  2:137048/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGGACTTCAATCTCCGGCGTGCGGCGTTGTTTGCGACCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTAC  >  NC_000913/1226683‑1226980

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: