Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,026,881 C→T R119C (CGC→TGC)  fre → NAD(P)H‑flavin reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,026,8810CT96.4% 320.9 / 0.9 83R119C (CGC→TGC) freNAD(P)H‑flavin reductase
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base T (45/35);  total (47/36)
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.57e-01

CAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCCGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAACCGTGATATCACCATTTACTGGGGCGGGCGTGAAGAGCAGCATCTGTATGATCTCTGCGAGCTTGAGGCGCTTTCGTTGAAGCATCCTGGTCTGCAAGTGGTGCCGGTGG  >  NC_000913/4026736‑4027028
                                                                                                                                                 |                                                                                                                                                   
cAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCg                                                                                                                                                <  2:184059/151‑1 (MQ=255)
cAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCg                                                                                                                                                >  1:956077/1‑151 (MQ=255)
cAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCg                                                                                                                                                >  2:1623366/1‑151 (MQ=255)
 aaCCTTTACGCGAAAGCAGTCATGGACAGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCg                                                                                                                                                >  1:603490/1‑150 (MQ=255)
  aCCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGAt                                                                                                                                              >  2:1009190/1‑151 (MQ=255)
  aCCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGAt                                                                                                                                              <  1:1581190/151‑1 (MQ=255)
    cTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGAttt                                                                                                                                            >  1:154054/1‑151 (MQ=255)
     tttACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGAtttt                                                                                                                                           >  2:810162/1‑151 (MQ=255)
      ttACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGATTTTg                                                                                                                                          >  2:943465/1‑151 (MQ=255)
       tACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGATTTTGc                                                                                                                                         >  1:814215/1‑151 (MQ=255)
        aCGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGATTTTGCt                                                                                                                                        <  2:804080/151‑1 (MQ=255)
          gcgAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGATTTTGCTGa                                                                                                                                      <  1:726394/151‑1 (MQ=255)
           cgAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGATTTTGCTGAc                                                                                                                                     >  1:306958/1‑151 (MQ=255)
                   gTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGATTTTGCTGACAGCGTTgg                                                                                                                             >  2:1082354/1‑151 (MQ=255)
                     cATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGATTTTGCTGACAGCGTTGgcg                                                                                                                           <  1:264427/151‑1 (MQ=255)
                        ggACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGATTTTGCTGACAGCGTTGGCGCGt                                                                                                                        >  2:327982/1‑151 (MQ=255)
                             gCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGTAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAAccc                                                                                                                   >  1:1415564/1‑151 (MQ=255)
                             gCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAAccc                                                                                                                   >  1:846733/1‑151 (MQ=255)
                              cATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAAccc                                                                                                                   >  1:925725/1‑150 (MQ=255)
                               aTCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCaa                                                                                                                 <  1:1588479/151‑1 (MQ=255)
                               aTCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCaa                                                                                                                 <  1:64064/151‑1 (MQ=255)
                               aTCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCa                                                                                                                  >  2:1116317/1‑150 (MQ=255)
                                 ccTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAAc                                                                                                               <  1:282232/151‑1 (MQ=255)
                                   tCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAACCg                                                                                                             >  2:1555875/1‑151 (MQ=255)
                                    cAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAAcc                                                                                                              <  2:68953/149‑1 (MQ=255)
                                    cAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAACCg                                                                                                             <  2:187432/150‑1 (MQ=255)
                                    cAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAACCg                                                                                                             <  1:778264/150‑1 (MQ=255)
                                    cAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAACCGt                                                                                                            <  1:604057/151‑1 (MQ=255)
                                      aaGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAACCGTg                                                                                                           >  2:202694/1‑150 (MQ=255)
                                      aaGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAACCGTGa                                                                                                          <  2:1658991/151‑1 (MQ=255)
                                      aaGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAACCGTGa                                                                                                          <  1:863670/151‑1 (MQ=255)
                                       aGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAACCGTGat                                                                                                         >  1:1357857/1‑151 (MQ=255)
                                        gATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAACCGTGata                                                                                                        <  2:1263820/151‑1 (MQ=255)
                                              caAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAACCGTGATATCACCa                                                                                                  <  2:196929/151‑1 (MQ=255)
                                                 aTCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAACCGTGATATCACCAttt                                                                                               <  2:1462085/151‑1 (MQ=255)
                                                     tGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAACCGTGATATCACCATTTACTg                                                                                           >  1:831951/1‑151 (MQ=255)
                                                      ggTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAACCGTGATATCACCATTTACTgg                                                                                          <  2:1063876/151‑1 (MQ=255)
                                                       gTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGTAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAACCGTGATATCACCATTTACTggg                                                                                         <  2:1415657/151‑1 (MQ=255)
                                                         cGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAACCGTGATATCACCATTTACTGgggc                                                                                       >  1:901735/1‑151 (MQ=255)
                                                           aCATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAACCGTGATATCACCATTTACTGgggcgg                                                                                     >  2:707162/1‑151 (MQ=255)
                                                             aTTCCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAACCGTGATATCACCATTTACTGgggcggg                                                                                    >  1:163130/1‑150 (MQ=255)
                                                               tCCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAACCGTGATATCACCATTTACTGGggcgggcg                                                                                  >  2:1213713/1‑150 (MQ=255)
                                                                   cACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAACCGTGATATCACCATTTACTGGGGCGGGCGTGAag                                                                             <  1:303659/151‑1 (MQ=255)
                                                                    aCGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAACCGTGATATCACCATTTACTGGGGCGGGCGTGa                                                                               >  1:366838/1‑148 (MQ=255)
                                                                       gagaAGCGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAACCGTGATATCACCATTTACTGGGGCGGGCGTGAAGagca                                                                         <  1:864583/151‑1 (MQ=255)
                                                                             cGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAACCGTGATATCACCATTTACTGGGGCGGGCGTGAAGAGCAGCATCt                                                                   <  1:1528646/151‑1 (MQ=255)
                                                                             cGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAACCGTGATATCACCATTTACTGGGGCGGGCGTGAAGAGCAGCATCt                                                                   >  1:587804/1‑151 (MQ=255)
                                                                                 gCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAACCGTGATATCACCATTTACTGGGGCGGGCGTGAAGAGCAGCATCTGTAt                                                               >  2:939789/1‑151 (MQ=255)
                                                                                   tGCGCGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAACCGTGATATCACCATTTACTGGGGCGGGCGTGAAGAGCAGCATCTGTATg                                                              >  2:305077/1‑150 (MQ=255)
                                                                                    gcgcGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAACCGTGATATCACCATTTACTGGGGCGGGCGTGAAGAGCAGCATCTGTATGAt                                                            <  2:937892/151‑1 (MQ=255)
                                                                                     cgcgATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAACCGTGATATCACCATTTACTGGGGCGGGCGTGAAGAGCAGCATCTGTATGa                                                             <  1:912298/149‑1 (MQ=255)
                                                                                     cgcgATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAACCGTGATATCACCATTTACTGGGGCGGGCGTGAAGAGCAGCATCTGTATGAtc                                                           >  1:1027748/1‑151 (MQ=255)
                                                                                      gcgATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAACCGTGATATCACCATTTACTGGGGCGGGCGTGAAGAGCAGCATCTGTATGAtct                                                          >  2:1334387/1‑151 (MQ=255)
                                                                                        gatgatGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAACCGTGATATCACCATTTACTGGGGCGGGCGTGAAGAGCAGCATCTGTATGATctct                                                        >  2:123433/1‑151 (MQ=255)
                                                                                        gatgatGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAACCGTGATATCACCATTTACTGGGGCGGGCGTGAAGAGCAGCATCTGTATGATctct                                                        <  2:1083525/151‑1 (MQ=255)
                                                                                        gatgatGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAACCGTGATATCACCATTTACTGGGGCGGGCGTGAAGAGCAGCATCTGTATGATctct                                                        >  1:973266/1‑151 (MQ=255)
                                                                                         atgatgAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCCGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAACCGTGATATCACCATTTACTGGGGCGGGCGTGAAGAGCAGCATCTGTATGATCTCTg                                                       <  1:1382386/151‑1 (MQ=255)
                                                                                              gaAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCCGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAACCGTGATATCACCATTTACTGGGGCGGGCGTGAAGAGCAGCATCTGTATGATCTCTGCGAGc                                                  >  1:1261918/1‑151 (MQ=255)
                                                                                                 gagCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAACCGTGATATCACCATTTACTGGGGCGGGCGTGAAGAGCAGCATCTGTATGATCTCTGCGAGCTTg                                               <  2:966288/151‑1 (MQ=255)
                                                                                                    cGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAACCGTGATATCACCATTTACTGGGGCGGGCGTGAAGAGCAGCATCTGTATGATCTCTGCGAGCTTGAgg                                            >  1:496192/1‑151 (MQ=255)
                                                                                                    cGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAACCGTGATATCACCATTTACTGGGGCGGGCGTGAAGAGCAGCATCTGTATGATCTCTGCGAGCTTGAgg                                            >  2:1380568/1‑151 (MQ=255)
                                                                                                     gTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAACCGTGATATCACCATTTACTGGGGCGGGCGTGAAGAGCAGCATCTGTATGATCTCTGCGAGCTTGAGgc                                           >  1:780377/1‑151 (MQ=255)
                                                                                                       ccGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAACCGTGATATCACCATTTACTGGGGCGGGCGTGAAGAGCAGCATCTGTATGATCTCTGCGAGCTTGAGgcg                                          <  1:414181/150‑1 (MQ=255)
                                                                                                        cGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCCGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAACCGTGATATCACCATTTACTGGGGCGGGCGTGAAGAGCAGCATCTGTATGATCTCTGCGAGCTTGAGGCGCt                                        >  1:1089318/1‑151 (MQ=255)
                                                                                                         gatgatTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAACCGTGATATCACCATTTACTGGGGCGGGCGTGAAGAGCAGCATCTGTATGATCTCTGCGAGCTTGAGGCGCtt                                       >  1:330514/1‑151 (MQ=255)
                                                                                                          atgatTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAACCGTGATATCACCATTTACTGGGGCGGGCGTGAAGAGCAGCATCTGTATGATCTCTGCGAGCTTGAGGCGCttt                                      >  1:848660/1‑151 (MQ=255)
                                                                                                                 tGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAACCGTGATATCACCATTTACTGGGGCGGGCGTGAAGAGCAGCATCTGTATGATCTCTGCGAGCTTGAGGCGCTTTCGTTGa                                <  1:731997/150‑1 (MQ=255)
                                                                                                                 tGATTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAACCGTGATATCACCATTTACTGGGGCGGGCGTGAAGAGCAGCATCTGTATGATCTCTGCGAGCTTGAGGCGCTTTCGTTGa                                <  2:700339/150‑1 (MQ=255)
                                                                                                                   aTTGCGGGCGGCACCGGGTTCTCTTATGCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAACCGTGATATCACCATTTACTGGGGCGGGCGTGAAGAGCAGCATCTGTATGATCTCTGCGAGCTTGAGGCGCTTTCGTTGAAGc                             >  1:1424275/1‑151 (MQ=255)
                                                                                                                       cgggcggCACCGGGTTCTCTTATGCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAACCGTGATATCACCATTTACTGGGGCGGGCGTGAAGAGCAGCATCTGTATGATCTCTGCGAGCTTGAGGCGCTTTCGTTGAAGCATcc                         >  1:273514/1‑151 (MQ=255)
                                                                                                                       cgggcggCACCGGGTTCTCTTATGCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAACCGTGATATCACCATTTACTGGGGCGGGCGTGAAGAGCAGCATCTGTATGATCTCTGCGAGCTTGAGGCGCTTTCGTTGAAGCATc                          <  2:1286270/150‑1 (MQ=255)
                                                                                                                        gggcggCACCGGGTTCTCTTATGCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAACCGTGATATCACCATTTACTGGGGCGGGCGTGAAGAGCAGCATCTGTATGATCTCTGCGAGCTTGAGGCGCTTTCGTTGAAGCATCCt                        >  2:1710587/1‑151 (MQ=255)
                                                                                                                                ccGGGTTCTCTTATGCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAACCGTGATATCACCATTTACTGGGGCGGGCGTGAAGAGCAGCATCTGTATGATCTCTGCGAGCTTGAGGCGCTTTCGTTGAAGCATCCTGGTCTGc                 >  1:1404101/1‑150 (MQ=255)
                                                                                                                                   ggTTCTCTTATGCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAACCGTGATATCACCATTTACTGGGGCGGGCGTGAAGAGCAGCATCTGTATGATCTCTGCGAGCTTGAGGCGCTTTCGTTGAAGCATCCTGGTCTGCAAgt             >  2:246396/1‑151 (MQ=255)
                                                                                                                                    gTTCTCTTATGCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAACCGTGATATCACCATTTACTGGGGCGGGCGTGAAGAGCAGCATCTGTATGATCTCTGCGAGCTTGAGGCGCTTTCGTTGAAGCATCCTGGTCTGCAAgtg            >  2:548780/1‑151 (MQ=255)
                                                                                                                                      tctcTTATGCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAACCGTGATATCACCATTTACTGGGGCGGGCGTGAAGAGCAGCATCTGTATGATCTCTGCGAGCTTGAGGCGCTTTCGTTGAAGCATCCTGGTCTGCAAgtgg           <  1:810097/150‑1 (MQ=255)
                                                                                                                                         ctTATGCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAACCGTGATATCACCATTTACTGGGGCGGGCGTGAAGAGCAGCATCTGTATGATCTCTGCGAGCTTGAGGCGCTTTCGTTGAAGCATCCTGGTCTGCAAgtggtg         <  1:1340203/149‑1 (MQ=255)
                                                                                                                                          ttATGCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAACCGTGATATCACCATTTACTGGGGCGGGCGTGAAGAGCAGCATCTGTATGATCTCTGCGAGCTTGAGGCGCTTTCGTTGAAGCATCCTGGTCTGCAAGTGGTGCCg      >  2:1060341/1‑151 (MQ=255)
                                                                                                                                            aTGCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAACCGTGATATCACCATTTACTGGGGCGGGCGTGAAGAGCAGCATCTGTATGATCTCTGCGAGCTTGAGGCGCTTTCGTTGAAGCATCCTGGTCTGCAAGTGGTGCCgg     <  1:1009114/150‑1 (MQ=255)
                                                                                                                                             tGCCTGCTCGATTTTGTTGACAGCGTTGGCGCGTAACCCAAACCGTGATATCACCATTTACTGGGGCGGGCGTGAAGAGCAGCATCTGTATGATCTCTGCGAGCTTGAGGCGCTTTCGTTGAAGCATCCTGGTCTGCAAGTGGTGCCggtg   >  1:1167147/1‑151 (MQ=255)
                                                                                                                                              gCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAACCGTGATATCACCATTTACTGGGGCGGGCGTGAAGAGCAGCATCTGTATGATCTCTGCGAGCTTGAGGCGCTTTCGTTGAAGCATCCTGGTCTGCAAGTGGTGCCggtgg  <  2:832013/151‑1 (MQ=255)
                                                                                                                                              gCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAACCGTGATATCACCATTTACTGGGGCGGGCGTGAAGAGCAGCATCTGTATGATCTCTGCGAGCTTGAGGCGCTTTCGTTGAAGCATCCTGGTCTGCAAGTGGTGCCggtgg  >  2:474369/1‑151 (MQ=255)
                                                                                                                                              gCCTGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAACCGTGATATCACCATTTACTGGGGCGGGCGTGAAGAGCAGCATCTGTATGATCTCTGCGAGCTTGAGGCGCTTTCGTTGAAGCATCCTGGTCTGCAAGTGGTGCCggtg   <  1:1635269/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGCCCGCTCGATTTTGCTGACAGCGTTGGCGCGTAACCCAAACCGTGATATCACCATTTACTGGGGCGGGCGTGAAGAGCAGCATCTGTATGATCTCTGCGAGCTTGAGGCGCTTTCGTTGAAGCATCCTGGTCTGCAAGTGGTGCCGGTGG  >  NC_000913/4026736‑4027028

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