Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,871,066 T→G Y228D (TAC→GAC)  yeaI → putative membrane‑anchored diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,871,0660TG97.3% 638.8 / 0.1 186Y228D (TAC→GAC) yeaIputative membrane‑anchored diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (3/2);  new base G (75/105);  total (79/107)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 187 total aligned reads displayed.

GGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCA  >  NC_000913/1870967‑1871162
                                                                                                   |                                                                                                
ggAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTAAGCAGTGATAACGCTGACTATTCTTTGTATGTCGTCGAtt                                                                                                  >  2:3251026/1‑100 (MQ=255)
  aaTCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATGAc                                                                                                <  1:1088800/100‑1 (MQ=255)
  aaTCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATGAc                                                                                                <  2:892161/100‑1 (MQ=255)
    tCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATGACTg                                                                                              <  1:605670/100‑1 (MQ=255)
    tCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATGACTg                                                                                              >  1:3759744/1‑100 (MQ=255)
    tCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATGACTg                                                                                              >  1:1854289/1‑100 (MQ=255)
       tGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATGCCTGTcc                                                                                           >  1:209622/1‑100 (MQ=255)
        gTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATGACTGTCCg                                                                                          <  1:5474453/100‑1 (MQ=255)
        gTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATGACTGTCCg                                                                                          <  1:205021/100‑1 (MQ=255)
           gCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATGACTGTCCGGAc                                                                                       >  1:2899202/1‑100 (MQ=255)
              ggCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATACCGCTGACTATTCTTTGTATGTCGTCGATGACTGTCCGGCCaca                                                                                    >  1:2785149/1‑98 (MQ=255)
              ggCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATGACTGTCCGGACAAc                                                                                    >  1:4547104/1‑100 (MQ=255)
              ggCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATGACTGTCCGGACAAc                                                                                    <  1:2184126/100‑1 (MQ=255)
                cgcTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATGACTGTCCGGACAACCa                                                                                  <  1:4078230/100‑1 (MQ=255)
                 gcTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATGACTGTCCGGACAACCAt                                                                                 <  2:1265289/100‑1 (MQ=255)
                  cTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATGACTGTCCGGACAACCATa                                                                                >  1:503779/1‑100 (MQ=255)
                  cTGAGCCCTTTTTTATTTTTTTTCCTTCTGGCACACAATCTGAGCAGTTTTAACGTTGACTATTCTTTGTATTTCGTCGATGTCTGTCCGGACCACtata                                                                                >  2:4668175/1‑100 (MQ=255)
                    gATCCCTTTTTTATTTTTTTTCCTTCTGGCACACAATCTCAGCAGTTATAACTCTCACTATTCTTTGTATGTCTCCGTTCACTGTCCTGACAACCATTCt                                                                              >  2:3084618/1‑100 (MQ=255)
                    gATCCCTTTTTTAGATTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATGCCTGTCCGGACAACCATACt                                                                              <  2:5028540/100‑1 (MQ=255)
                     aTCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACtn                                                                             <  1:799242/100‑2 (MQ=255)
                      tCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAATGCTGACTATTCTTTGTATGTCGTCGATGACTGTCCGGACAACCATACTGc                                                                            <  2:5411953/100‑1 (MQ=255)
                      tCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATGACTGTCCGGACAACCATACTGc                                                                            >  2:3850004/1‑100 (MQ=255)
                       cccTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATGACTGTCCGGACAACCATACTGCg                                                                           >  1:1393018/1‑100 (MQ=255)
                       cccTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATGACTGTCCGGACAACCATACTGCg                                                                           <  2:4499940/100‑1 (MQ=255)
                          ttttttAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATGACTGTCCGGACAACCATACTGCGAcc                                                                        <  1:1691285/100‑1 (MQ=255)
                          ttttttAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATGACTGTCCGGACAACCATACTGCGAcc                                                                        <  1:2838124/100‑1 (MQ=255)
                          ttttttAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATGACTGTCCGGACAACCATACTGCGAcc                                                                        <  1:1356812/100‑1 (MQ=255)
                             tttAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATGACTGTCCGGACAACCATACTGCGACCTgg                                                                     >  1:2081905/1‑100 (MQ=255)
                              ttAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATGACTGTCCGGACAACCATACTGCGACCTggg                                                                    >  2:246366/1‑100 (MQ=255)
                              ttAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATGACTGTCCGGACAACCATACTGCGACCTggg                                                                    <  2:3508998/100‑1 (MQ=255)
                               tAGTTTTTGCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATGACTGTCCGGACAACCATACTGCGACCTgggg                                                                   <  1:12112/100‑1 (MQ=255)
                               tAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATGACTGTCCGGACAACCATACTGCGACCTgggg                                                                   <  2:4795903/100‑1 (MQ=255)
                                aGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATGACTGTCCGGACAACCATACTGCGACCTGGGGa                                                                  >  2:2567561/1‑100 (MQ=255)
                                 gTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATGACTGTCCGGACAACCATACTGCGACCTGGGGaa                                                                 >  1:1020952/1‑100 (MQ=255)
                                 gTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATGACTGTCCGGACAACCATACTGCGACCTGGGGaa                                                                 >  2:814364/1‑100 (MQ=255)
                                   ttttCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATGACTGTCCGGACAACCATACTGCGACCTGGGGAATc                                                               >  2:2313042/1‑100 (MQ=255)
                                   ttttCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATGACTGTCCGGACAACCATACTGCGACCTGGGGAATc                                                               >  1:1278034/1‑100 (MQ=255)
                                    tttCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATGACTGTCCGGACAACCATACTGCGACCTGGGGAATCa                                                              <  2:2791588/100‑1 (MQ=255)
                                    tttCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATGACTGTCCGGACAACCATACTGCGACCTGGGGAATCa                                                              <  2:4519429/100‑1 (MQ=255)
                                    tttCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATGACTGTCCGGACAACCATACTGCGACCTGGGGAATCa                                                              >  1:1828527/1‑100 (MQ=255)
                                      tCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATGACTGTCCGGACAACCATACTGCGACCTGGGGAATCaan                                                            <  1:4447184/100‑2 (MQ=255)
                                       tccgatCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATGACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACt                                                           <  1:4516636/95‑1 (MQ=255)
                                       ccccTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATGACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACt                                                           <  2:2514890/100‑1 (MQ=255)
                                       ccccTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATGACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACt                                                           <  1:2256530/100‑1 (MQ=255)
                                       ccccTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATGACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACt                                                           >  2:2096759/1‑100 (MQ=255)
                                       ccccTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATGACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACt                                                           <  1:5130918/100‑1 (MQ=255)
                                        cccTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATGACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACta                                                          <  2:1833546/100‑1 (MQ=255)
                                        cccTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATGACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACta                                                          <  1:3230073/100‑1 (MQ=255)
                                         ccTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATGACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACtat                                                         >  2:3498779/1‑100 (MQ=255)
                                          cTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATGACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACtata                                                        <  1:734657/100‑1 (MQ=255)
                                          cTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATGACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACtata                                                        >  1:3611418/1‑100 (MQ=255)
                                           ttCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATGACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATAc                                                       <  2:4118800/100‑1 (MQ=255)
                                           ttCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATGACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATAc                                                       <  2:207601/100‑1 (MQ=255)
                                              tGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATGACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACaaa                                                    <  2:3609182/100‑1 (MQ=255)
                                                gCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATGACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACaaaaa                                                  <  1:271450/100‑1 (MQ=255)
                                                 cacacaATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATGACTGTCCGGACAACCATACTGTGACCTGGGGAATCAACTATACAAAAat                                                 >  1:2368567/1‑100 (MQ=255)
                                                 cacacaATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATGACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACGAAAat                                                 >  1:5271455/1‑100 (MQ=255)
                                                     caATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATGACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTg                                             >  1:1518631/1‑100 (MQ=255)
                                                      aaTCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATGACTGTACGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTgg                                            >  2:4979357/1‑100 (MQ=255)
                                                        tcggAGCAGGTATAGCGCTGACTATTCGTTGTATGTCGGCGATGACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGtt                                          <  1:1243918/97‑1 (MQ=255)
                                                         cTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATGACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGttt                                         <  1:825673/100‑1 (MQ=255)
                                                          tgggCAGTAGTAACGCTGACTACTCGTTGTATGGCGTCGATGACTGTCCGGACACCCATACTGCAACCTGGGGAATCATCTATAGAAAAATATTGGTTTg                                        <  2:4251028/97‑1 (MQ=255)
                                                          tGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATGACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTg                                        >  2:1072053/1‑100 (MQ=255)
                                                          tGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATGACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTg                                        >  1:437928/1‑100 (MQ=255)
                                                             gCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATGCCTGTCCTGACAACCATACTGCGACCTGGTGCATCTAATATGCGAAAATATTGATTTGACt                                     >  2:1469377/1‑100 (MQ=255)
                                                                gTTATAACGCTGACTATTCTTTGTATGTCGTCGATGACTCTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGATTGTCtgtg                                  >  1:4860807/1‑100 (MQ=255)
                                                                 ttATAACGCTGACTATTCTTTGTATGTCGTCGATGACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTgg                                 <  1:2337174/100‑1 (MQ=255)
                                                                 ctataGCGCTGACCATTCTTTGTATGTCGTCGATGACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTgg                                 <  2:5124629/99‑1 (MQ=255)
                                                                  tataACGCTGACTATTCTTTGTATGTCGTCGATGACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTggg                                <  2:4745255/100‑1 (MQ=255)
                                                                   ataACGCTGACTATTCTTTGTATGTCGTCGATGACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGc                               <  2:1742208/100‑1 (MQ=255)
                                                                      ncGCTGACTATTCTTTGTATGTCGTCGATGACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAACTATTGGTTTGTCTGTGGGCAtt                            >  1:5119470/2‑100 (MQ=255)
                                                                      aCGCTGACTATTCTTTGTATGTCGTCGATGACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCAtt                            >  2:114636/1‑100 (MQ=255)
                                                                         cTGACTATTCTTTGTATGTCGTCGATGACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCAttttt                         <  2:1007407/100‑1 (MQ=255)
                                                                         cTGACTATTCTTTGTATGTCGTCGATGACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCAttttt                         >  2:6298/1‑100 (MQ=255)
                                                                           gACTATTCTTTGTATGTCGTCGATGACTGTCCGGACATCCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTAc                       <  2:4998004/100‑1 (MQ=255)
                                                                           gACTATTCTTTGTATGTCGTCGATGACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTAc                       <  2:4789819/100‑1 (MQ=255)
                                                                           gACTATTCTTTGTATGTCGTCGATGACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTAc                       <  1:2603812/100‑1 (MQ=255)
                                                                           gACTATTCTTTGTATGTCGTCGATGACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTAc                       >  1:2157260/1‑100 (MQ=255)
                                                                            aCTATTCTTTGTATGTCGTCGATGACTGTCCGGAAAACCATGCTGAGACATGGGGAATCAACTATCCACCAATCTTGGTTTGTCTGTGGGCCTTTTTACt                      <  2:3026943/100‑1 (MQ=255)
                                                                              tATTCTTTGTATGTCGTCGATGACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGt                    >  1:1294635/1‑100 (MQ=255)
                                                                              tATTCTTTGTATGTCGTCGATGACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGt                    <  2:4328981/100‑1 (MQ=255)
                                                                              tATTCTTTGTATGTCGTCGATGACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGt                    >  1:1294599/1‑100 (MQ=255)
                                                                                  cTTTGTATGTCGTCGATGACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCtt                <  2:3412566/100‑1 (MQ=255)
                                                                                  cTTTGTATGTCGTCGATGACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCtt                >  1:3406501/1‑100 (MQ=255)
                                                                                  cTTTGTATGTCGTCGATGACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCtt                >  1:3196531/1‑100 (MQ=255)
                                                                                  cTTTGTATGTCGTCGATGACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCtt                >  2:824506/1‑100 (MQ=255)
                                                                                  cTTTGTATGTCGTCGATGACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCtt                <  1:2647293/100‑1 (MQ=255)
                                                                                  cTTTGTATGTCGTCGATGACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCtt                >  2:2621533/1‑100 (MQ=255)
                                                                                     tGTATGTCGTCGATGACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTttat             <  2:3578087/100‑1 (MQ=255)
                                                                                       tATGTCGTCGATGACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTttatta           <  1:1241264/100‑1 (MQ=255)
                                                                                        aTGTCGTCGATGACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTtattat          <  2:1499458/100‑1 (MQ=255)
                                                                                         tGTCGTCGATGACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATc         >  2:3364357/1‑100 (MQ=255)
                                                                                          gtcgtcGATGACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCa        <  1:4929361/100‑1 (MQ=255)
                                                                                          gtcgtcGATGACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCa        <  1:1138422/100‑1 (MQ=255)
                                                                                          gtcgtcGATGACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCa        <  2:2951151/100‑1 (MQ=255)
                                                                                           tcgtcgATGACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCAt       <  1:4814483/100‑1 (MQ=255)
                                                                                             gtcgATGACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATgc     <  1:2689548/100‑1 (MQ=255)
                                                                                             gtcgATGACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATgc     <  2:4486131/100‑1 (MQ=255)
                                                                                             gtcgATGACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATgc     >  1:91472/1‑100 (MQ=255)
                                                                                              tcgATGACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATgcg    <  2:2885072/100‑1 (MQ=255)
                                                                                                   |                                                                                                
GGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCA  >  NC_000913/1870967‑1871162

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