Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,478,431 C→T 100% V45I (GTA→ATA)  PP_4815 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,478,4310CT95.6% 213.1 / ‑0.2 68V45I (GTA→ATA) PP_4815membrane protein
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base T (29/36);  total (31/37)
Fisher's exact test for biased strand distribution p-value = 5.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GGTTAAGATAGAGCGGCACAGCCAGGGCCACCGTAGCGGGCCCCAGCAGGATGTTCATGATTTCGGTGCTTTTGCGGTACTCGCTGTAATCGATGCCGCACAGCAGCAACACACCGATCACCAGCAGCATCGATACCAATACAGGTTGCAGGAAGATCCAGCGGGTTTTCTCGTAGGCTGCCAGCACTACCTGATAGGCCCCCAGCGTGATGCCGATACCGAACAATGGGTGATGAATGACCGCATCGAGCGCGCCCTGCCAATCGAGCATCATGGCTGCTCCTCACG  >  NC_002947/5478290‑5478577
                                                                                                                                             |                                                                                                                                                  
ggTTAAGATAGAGCGGCACAGCCAGGGCCACCGTAGCGGGCCCCAGCAGGATGTTCATGATTTCGGTGCTTTTGCGGTACTCGCTGTAATCGATGCCGCACAGCAGCAACACACCGATCACCAGCAGCATCGATACCAATATAGGTTGCa                                                                                                                                            <  1:166877/150‑1 (MQ=255)
      gATAGAGCGGCACAGCCAGGGCCACCGTAGCGGGCCCCAGCAGGATGTTCATGATTTCGGTGCTTTTGCGGTACTCGCTGTAATCGATGCCGCACAGCAGCAACACACCGATCACCAGCAGCATCGATACCAATATAg                                                                                                                                                  <  1:642744/138‑1 (MQ=255)
        tAGAGCGGCACAGCCAGGGCCACCGTAGCGGGCCCCAGCAGGATGTTCATGATTTCGGTGCTTTTGCGGTACTCGCTGTAATCGATGCCGCACAGCAGCAACACACCGATCACCAGCAGCATCGATACCAATATAGGTTGCAGGAAGATc                                                                                                                                    >  1:556026/1‑150 (MQ=255)
             cGGCACAGCCAGGGCCACCGTAGCGGGCCCCAGCAGGATGTTCATGATTTCGGTGCTTTTGCGGTACTCGCTGTAATCGATGCCGCACAGCAGCAACACACCGATCACCAGCAGCATCGATACCAATATAGGTTGCAGGAAGATCCAGCg                                                                                                                               <  1:1002910/150‑1 (MQ=255)
              ggCACAGCCAGGGCCACCGTAGCGGGCCCCAGCAGGGTGTTCATGATTTCGGTGCTTTTGCGGTACTCGCTGTAATCGATGCCGCACAGCAGCAACACACCGATCACCAGCAGCATCGATACCAATATAGGTTGCAGGAAGATCCAGCgg                                                                                                                              <  2:1315435/150‑1 (MQ=255)
                cacaGCCAGGGCCACCGTAGCGGGCCCCAGCAGGATGTTCATGATTTCGGTGCTTTTGCGGTACTCGCTGTAATCGATGCCGCACAGCAGCAACACTCCGATCACCAGCAGCATCGATACCAATATAGGTTGCAGGAAGATCCAGCGGGt                                                                                                                            <  2:2151193/150‑1 (MQ=255)
                    gCCAGGGCCACCGTAGCGGGCCCCAGCAGGATGTTCATGATTTCGGTGCTTTTGCGGTACTCGCTGTAATCGATGCCGCACAGCAGCGACACACCGATCACCAGCAGCATCGATACCAATATAGGTTGCAGGAAGATCCAGCGGGTTTtc                                                                                                                        <  2:1416661/150‑1 (MQ=255)
                     ccAGGGCCACCGTAGCGGGCCCCAGCAGGATGTTCATGATTTCGGTGCTTTTGCGGTACTCGCTGTAATCGATGCCGCACAGCAGCAACACACCGATCACCAGCAGCATCGATACCAATATAGGTTGCAGGAAGATCCAGCGGGTTTtct                                                                                                                       >  1:209436/1‑150 (MQ=255)
                         ggCCACCGTAGCGGGCCCCAGCAGGATGTTCATGATTTCGGTGCTTTTGCGGTACTCGCTGTAATCGATGCCGCACAGCAGCAACACACCGATCACCAGCAGCATCGATACCAATATAGGTTGCAGGAAGATCCAGCGGGTTTTCTCGTa                                                                                                                   <  2:209436/150‑1 (MQ=255)
                          gCCACCGTAGCGGGCCCCAGCAGGATGTTCATGATTTCGGTGCTTTTGCGGTACTCGCTGTAATCGATGCCGCACAGCAGCAACACACCGATCACCAGCAGCATCGATACCAATATAGGTTGCAGGAAGATCCAGCGGGTTTTCTCGTAg                                                                                                                  >  2:1950816/1‑150 (MQ=255)
                          gCCACCGTAGCGGGCCCCAGCAGGATGTTCATGATTTCGGTGCTTTTGCGGTACTCGCTGTAATCGATGCCGCACAGCAGCAACACACCGATCACCAGCAGCATCGATACCAATATAGGTTGCAGGAAGATCCAGCGGGTTTTCTCGTAg                                                                                                                  <  2:556026/150‑1 (MQ=255)
                            cACCGTAGCGGGCCCCAGCAGGATGTTCATGATTTCGGTGCCTTTGCGGTACTCGCTGTAATCGATGCCGCACAGCAGCAACACACCGATCACCAGCAGCATCGATACCAATATAGGTTGCAGGAAGATCCAGCGGGTTTTCTCGTAGGc                                                                                                                >  1:2181938/1‑150 (MQ=255)
                             aCCGTAGCGGGCCCCAGCAGGATGTTCATGATTTCGGTGCTTTTGCGGTACTCGCTGTAATCGATGCCGCACAGCAGCAACACACCGATCACCAGCAGCATCGATACCAATATAGGTTGCAGGAAGATCCAGCGGGTTTTCTCGTAGGCt                                                                                                               >  2:1304473/1‑150 (MQ=255)
                                  aGCGGGCCCCAGCAGGATGTTCATGATTTCGGTGCTTTTGCGGTACTCGCTGTAATCGATGCCGCACAGCAGCAACACACCGATCACCAGCAGCATCGATACCAATATAGGTTGCAGGAAGATCCAGCGGGTTTTCTCGTAGGCTGCCAg                                                                                                          <  1:2163361/150‑1 (MQ=255)
                                   gCGGGCCCCAGCAGGATGTTCATGATTTCGGTGCTTTTGCGGTACTCGCTGTAATCGATGCCGCACAGCAGCAACACACCGATCACCAGCAGCATCGATACCAATATAGGTTGCAGGAAGATCCAGCGGGTTTTCTCGTAGGCTGCCAGc                                                                                                         <  2:189615/150‑1 (MQ=255)
                                         cccAGCAGGATGTTCATGATTTCGGTGCTTTTGCGGTACTCGCTGTAATCGATGCCGCACAGCAGCAACACACCGATCACCAGCAGCATCGATACCAATATAGGTTGCAGGAAGATCCAGCGGGTTTTCTCGTAGGCTGCCAGCACTAcc                                                                                                   <  1:1304473/150‑1 (MQ=255)
                                            agcagGATGTTCATGATTTCGGTGATTTTGCGGTACTCGCTGTAATCGATGCCGCACAGCAGCAACACACCGATCACCAGCAGCATCGATACCAATATAGGTTGCAGGAAGATCCAGCGGGTTTTCTCGTAGGCTGCCAGCACTACCTGa                                                                                                >  1:375558/1‑150 (MQ=255)
                                             gcagGATGTTCATGATTTCGGTGCTTTTGCGGTACTCGCTGTAATCGATGCCGCACAGCAGCAACACACCGATCACCAGCAGCATCGATACCAATATAGGTTGCAGGAAGATCCAGCGGGTTTTCTCGTAGGCTGCCAGCACTACCTGAt                                                                                               <  2:738063/150‑1 (MQ=255)
                                                    gTTCATGATTTCGGTGCTTTTGCGGTACTCGCTGTAATCGATGCCGCACAGCAGCAACACACCGATCACCAGCAGCATCGATACCAATATAGGTTGCAGGAAGATCCAGCGGGTTTTCTCGTAGGCTGCCAGCACTACCTGATAGGcccc                                                                                        <  1:1389898/150‑1 (MQ=255)
                                                      tCATGATTTCGGTGCTTTTGCGGTACTCGCTGTAATCGATGCCGCACAGCAGCAACACACCGATCACCAGCAGCATCGATACCAATATAGGTTGCAGGAAGATCCAGCGGGTTTTCTCGTAGGCTGCCAGCACTa                                                                                                     >  2:2055079/1‑135 (MQ=255)
                                                      tCATGATTTCGGTGCTTTTGCGGTACTCGCTGTAATCGATGCCGCACAGCAGCAACACACCGATCACCAGCAGCATCGATACCAATATAGGTTGCAGGAAGATCCAGCGGGTTTTCTCGTAGGCTGCCAGCACTa                                                                                                     <  1:2055079/135‑1 (MQ=255)
                                                      tCATGATTTCGGTGCTTTTGCGGTACTCGCTGTAATCGATGCCGCACAGCAGCAACACACCGATCACCAGCAGCATCGATACCAATACAGGTTGCAGGAAGATCCAGCGGGTTTTCTCGTAGGCTGCCAGCACTACCTGATAGGCCCCCa                                                                                      >  1:2047548/1‑150 (MQ=255)
                                                         tGATTTCGGTGCTTTTGCGGTACTCGCTGTAATCGATGCCGCACAGCAGCAACACACCGATCACCAGCAGCATCGATACCAATATAGGTTGCAGGAAGATCCAGCGGGTTTTCTCGTAGGCTGCCAGCACTACCTGATAGGc                                                                                           <  2:1769963/142‑1 (MQ=255)
                                                         tGATTTCGGTGCTTTTGCGGTACTCGCTGTAATCGATGCCGCACAGCAGCAACACACCGATCACCAGCAGCATCGATACCAATATAGGTTGCAGGAAGATCCAGCGGGTTTTCTCGTAGGCTGCCAGCACTACCTGATAGGc                                                                                           >  1:1769963/1‑142 (MQ=255)
                                                         tGATTTCGGTGCTTTTGCGGTACTCGCTGTAATCGATGCCGCACAGCAGCAACACACCGATCACCAGCAGCATCGATACCAATACAGGTTGCAGGAAGATCCAGCGGGTTTTCTCGTAGGCTGCCAGCACTACCTGATAGGCCCCCAGCg                                                                                   <  2:2047548/150‑1 (MQ=255)
                                                         tGATTTCGGTGCTTTTGCGGTACTCGCTGTAATCGATACCGCACAGCAGCAACACACCGATCACCAGCAGCATCGATACCAATATAGGTTGCAGGAAGATCCAGCGGGTTTTCTCGTAGGCTGCCAGCACTACCTGATAGGCCCCCAGCg                                                                                   >  2:558078/1‑150 (MQ=255)
                                                          gATTTCGGTGCTTTTGCGGTACTCGCTGTAATCGATGCCGCACAGCAGCAACACACCGATCACCAGCAGCATCGATACCAATATAGGTTGCAGGAAGATCCAGCGGGTTTTCTCGTa                                                                                                                   <  1:1864226/117‑1 (MQ=255)
                                                          gATTTCGGTGCTTTTGCGGTACTCGCTGTAATCGATGCCGCACAGCAGCAACACACCGATCACCAGCAGCATCGATACCAATATAGGTTGCAGGAAGATCCAGCGGGTTTTCTCGTa                                                                                                                   >  2:1864226/1‑117 (MQ=255)
                                                               cGGTGCTTTTGCGGTACTCGCTGTAATCGATGCCGCACAGCAGCAACACACCGATCACCAGCAGCATCGATACCAATATAGGTTGCAGGAAGATCCAGCGGGTTTTCTCGTAGGCTGCCAGCACTACCTGATAGGCCCCCAGCGTGATGc                                                                             <  2:910350/150‑1 (MQ=255)
                                                                ggTGCTTTTGCGGTACTCGCTGTAATCGATGCCGCACAGCAGCAACACACCGATCACCAGCAGCATCGATACCAATATAGGTTGCAGGAAGATCCAGCGGGTTTTCTCGTAGGCTGCCAGCACTACCTGATAGGCCCCCAGCGTGATGcc                                                                            <  1:2265273/150‑1 (MQ=255)
                                                                 gTGCTTTTGCGGTACTCGCTGTAATCGATGCCGCACAGCAGCAACACACCGATCACCAGCAGCATCGATACCAATATAGGTTGCAGGAAGATCCAGCGGGTTTTCTCGTAGGCTGCCAGCACTACCTGATAGGCCCCCAGc                                                                                    >  2:1467168/1‑141 (MQ=255)
                                                                 gTGCTTTTGCGGTACTCGCTGTAATCGATGCCGCACAGCAGCAACACACCGATCACCAGCAGCATCGATACCAATATAGGTTGCAGGAAGATCCAGCGGGTTTTCTCGTAGGCTGCCAGCACTACCTGATAGGCCCCCAGc                                                                                    <  1:1467168/141‑1 (MQ=255)
                                                                  tGCTTTTGCGGTACTCGCTGTAATCGATGCCGCACAGCAGCAACACACCGATCACCAGCAGCATCGATACCAATATAGGTTGCAGGAAGATCCAGCGGGTTTTCTCGTAGGCTGCCAGCACTACCTGATAGGCCCCCAGCGTGATGCCGa                                                                          <  2:1179979/150‑1 (MQ=255)
                                                                     ttttGCGGTACTCGCTGTAATCGATGCCGCACAGCAGCAACACACCGATCACCAGCAGCATCGATACCAATATCGGTTGCAGGAAGATCCAGCGGGTTTTCTCGTAGGCTGCCAGCACTACCTGATAGGCCCCCAGCGTGATGCCGATAc                                                                       >  1:1974042/1‑150 (MQ=255)
                                                                     ttttGCGGTACTCGCTGTAATCGATGCCGCACAGCAGCAACACACCGATCACCAGCAGCATCGATACCAATATAGGTTGCAGGAAGATc                                                                                                                                    <  2:2278884/89‑1 (MQ=255)
                                                                     ttttGCGGTACTCGCTGTAATCGATGCCGCACAGCAGCAACACACCGATCACCAGCAGCATCGATACCAATATAGGTTGCAGGAAGATc                                                                                                                                    >  1:2278884/1‑89 (MQ=255)
                                                                             tACTCGCTGTAATCGATGCCGCACAGCAGCAACACACCGATCACCAGCAGCATCGATACCAATATAGGTTGCAGGAAGATCCAGCGGGTTTTCTCGTAGGCTGCCAGCACTACCTGATAGGCCCCCAGCGTGATGCCGATACCGAACAAt                                                               <  1:1518847/150‑1 (MQ=255)
                                                                              aCTCGCTGTAATCGATGCCGCACAGCAGCAACACACCGATCACCAGCAGCATCGATACCAATATAGGTTGCAGGAAGATCCAGCGGGTTTTCTCGTAg                                                                                                                  >  1:95916/1‑98 (MQ=255)
                                                                              aCTCGCTGTAATCGATGCCGCACAGCAGCAACACACCGATCACCAGCAGCATCGATACCAATATAGGTTGCAGGAAGATCCAGCGGGTTTTCTCGTAg                                                                                                                  <  2:95916/98‑1 (MQ=255)
                                                                                  gCTGTAATCGATGCCGCACAGCAGCAACACACCGCTCTCCAGCAGCATCGATACCAATATAGGTTGCAGGAAGATCAAGCGTGGTTTCTCGTAGGCTGCCAGCACTACCTGATAGGCCCCCAGCGTGATGCCGATACCGAACAATGGGtg                                                          >  1:1355270/1‑150 (MQ=255)
                                                                                  gCTGTAATCGATGCCGCACAGCAGCAACACACCGATCACCAGCAGCATCGATACCAATATAGGTTGCAGGAAGATCCAGCGGGTTTTCTCGTAGGCTGCCAGCACTACCTGATAGGCCCCCAGCGTGATGCCGATACCGAACAATGGGtg                                                          >  1:1281144/1‑150 (MQ=255)
                                                                                      tAATCGATGCCGCACAGCAGCAACACACCGATCACCAGCAGCATCGATACCAATATAGGTTGCAGGAAGATCCAGCGGGTTTTCTCGTAGGCTGCCAGCACTACCTGATAGGCCCCCAGCGTGATGcc                                                                            <  1:413181/128‑1 (MQ=255)
                                                                                      tAATCGATGCCGCACAGCAGCAACACACCGATCACCAGCAGCATCGATACCAATATAGGTTGCAGGAAGATCCAGCGGGTTTTCTCGTAGGCTGCCAGCACTACCTGATAGGCCCCCAGCGTGATGcc                                                                            >  2:413181/1‑128 (MQ=255)
                                                                                      tAATCGATGCCGCACAGCAGCAACACACCGATCACCAGCAGCATCGATACCAATATAGGTTGCAGGAAGATCCAGCGGGTTTTCTCGTAGGCTGCCAGCACTACCTGATAGGCCCCCAGCGTGATGCCGATACCGAACAATGGGtgatga                                                      >  2:751329/1‑150 (MQ=255)
                                                                                      tAATCGATGCCGCACAGCAGCAACACACCGATCACCAGCAGCATCAATACCAATATAGGTTGCAGGAAGATCCAGCGGGTTTTCTCGTAGGCTGCCAGCACTACCTGATAGGCCCCCAGCGTGATGCCGATACCGAACAATGGGtgatga                                                      >  2:1880349/1‑150 (MQ=255)
                                                                                                 gCACAGCAGGAACACACCGATCACCAGCAGCATCGATACCAATATAGGTTGCAGGAAGATCCAGCGGGTTTTCTCGTAGGCTGCCAGCACTACCTGATAGGCCCCCAGCGTGATGCCGATACCGAACAATGGGTGATGAATGACCGCATc                                           <  1:307225/150‑1 (MQ=255)
                                                                                                         gcaACACACCGATCACCAGCAGCATCGATACCAATATAGGTTGCAGGAAGATCCAGCGGGTTTTCTCGTAGGCTGc                                                                                                             >  1:1505444/1‑76 (MQ=255)
                                                                                                         gcaACACACCGATCACCAGCAGCATCGATACCAATATAGGTTGCAGGAAGATCCAGCGGGTTTTCTCGTAGGCTGc                                                                                                             <  2:1505444/76‑1 (MQ=255)
                                                                                                         gcaACACACCGATCACCAGCAGCATCGATACCAATATAGGTTGCAGGAAGATCCAGCGGGTTTTCTCGTAGGCTGCCAGCACTACCTGATAGGCCCCCAGCGTGATGCCGATACCGAACAATGGGTGATGAATGACCGCATCGAgcgcgc                                   <  1:1435654/150‑1 (MQ=255)
                                                                                                            acacacCGATCACCAGCAGCATCGATACCAATATAGGTTGCAGGAAGATCCAGCGGGTTTTCTCGTAGGCTGCCAGCACTACCTGATAGGCCCCCAGCGTGATGCCGATACCGAACAATGGGTGATGAATTACCGCATCGAGCGCGCCCt                                <  2:2095869/150‑1 (MQ=255)
                                                                                                               cacCGATCACCAGCAGCATCGATACCAATATAGGTTGCAGGAAGATCCAGCGGGTTTTCTCGTAGGCTGCCAGCACTACCTGATAGGCCCCCAGCGTGATGCCGATACCGAACAATGGGTGATGAATGACCGCATCGAGCGCGCCCTGcc                             >  1:1107668/1‑150 (MQ=255)
                                                                                                               cacCGATCACCAGCAGCATCGATACCAATATAGGTTGCAGGAAGATCCAGCGGGTTTTCTCGTAGGCTGCCAGCACTACCTGATAGGCCCCCAGCGTGATGCCGATACCGAACAATGGGTGATGAATGACCGCATCGAGCGCGCCCTGcc                             >  1:307937/1‑150 (MQ=255)
                                                                                                                acCGATCACCAGCAGCATCGATACCAATATAGGTTGCAGGAAGATCCAGCGGGTTTTCTCGTAGGCTGCCAGCACTACCTGATAGGCCCCCAGCGTGATGCCGATACCGAACAATGGGTGATGAATGACCGCATCGAGCGCGCCCTGCCa                            >  2:1314220/1‑150 (MQ=255)
                                                                                                                acCGATCACCAGCAGCATCGATACCAATATAGGTTGCAGGAAGATCCAGCGGGTTTTCTCGTAGGCTGCCAGCACTACCTGATAGGCCCCCAGCGTGATGCCGATACCGAACAATGGGTGATGAATGACCGCATCGAGCGCGCCCTGCCa                            >  1:1356848/1‑150 (MQ=255)
                                                                                                                   gATCACCAGCAGCATCGATACCAATATAGGTTGCAGGAAGATCCAGCGGGTTTTCTCGTAGGCTGCCAGCACTACCTGATAGGCCCCCAGCGTGATGCCGATACCGAACAATGGGTGATGAATGACCGCATCGAGCGCGCCCTGCCAATc                         <  1:751329/150‑1 (MQ=255)
                                                                                                                     tCACCAGCAGCATCGATACCAATATAGGTTGCAGGAAGATCCAGCGGGTTTTCTCGTAGGCTGCCAGCACTACCTGATAGGCCCCCAGCGTGATGCCGATACCGAACAATGGGTGATGAATGAcc                                                <  1:833884/125‑1 (MQ=255)
                                                                                                                     tCACCAGCAGCATCGATACCAATATAGGTTGCAGGAAGATCCAGCGGGTTTTCTCGTAGGCTGCCAGCACTACCTGATAGGCCCCCAGCGTGATGCCGATACCGAACAATGGGTGATGAATGAcc                                                >  2:833884/1‑125 (MQ=255)
                                                                                                                     tCACCAGCAGCATCGATACCAATATAGGTTGCAGGAAGATCCAGCGGGTTTTCTCGTAGGCTGCCAGCACTACCTGATAGGCCCCCAGCGTGATGCCGATACCGAACAATGGGTGATGAATGACCGCATCGAGCGCGCCCTGCCAATCGa                       <  2:2181938/150‑1 (MQ=255)
                                                                                                                     tCACCAGCAGCATCGATACCAATATAGGTTGCAGGAAGATCCAGCGGGTTTTCTCGTAGGCTGCCAGCACTACCTGATAGGCCCCCAGCGTGATGCCGATACCGAACAATGGGTGATGAATGACCGCATCGAGCGCGCCCTGCCAATCGa                       <  1:1314220/150‑1 (MQ=255)
                                                                                                                       aCCAGCAGCATCGATACCAATATAGGTTGCAGGAAGATCCAGCGGGTTTTCTCGTAGGCTGCCAGCACTACCTGATAGGCCCCCAGCGTGATGCCGATACCGAACAATGGGTGATGAATGACCGCATCGAGCGCGCCCTGCCAATCGAGc                     <  2:375558/150‑1 (MQ=255)
                                                                                                                        ccAGCAGCATCGATACCAATATAGGTTGCAGGAAGATCCAGCGGGTTTTCTCGTAGGCTGCCAGCACTACCTGATAGGCCCCCAGCGTGATGCCGATACCGAACAATg                                                              >  1:761714/1‑108 (MQ=255)
                                                                                                                        ccAGCAGCATCGATACCAATATAGGTTGCAGGAAGATCCAGCGGGTTTTCTCGTAGGCTGCCAGCACTACCTGATAGGCCCCCAGCGTGATGCCGATACCGAACAATg                                                              <  2:761714/108‑1 (MQ=255)
                                                                                                                             agcaTCGATACCAATATAGGTTGCAGGAAGATCCAGCGGGTTTTCTCGTAGGCTGCCAGCACTACCTGATAGGCCCCCAGCGTGATGCCGATACCGAACAATGGGTGATGAATGACCGCATCGAGCGCGCCCTGCCAATCGAGCATCATg               <  1:125591/150‑1 (MQ=255)
                                                                                                                                 tCGATACCAATATAGGTTGCAGGAAGATCCAGCGGGTTTTCTCGTAGGCTGCCAGCACTACCTGATAGGCCCCCAGCGTGATGCCGATACCGAACAATGGGTGATGAATGACCGCATCGAGCGCGCCCTGCCAATCGAGCATCATGgctg           >  1:245775/1‑150 (MQ=255)
                                                                                                                                        cAATATAGGTTGCAGGAAGATCCAGCGGGTTTTCTCGTAGGCTGCCAGCACTACCTGATAGGCCCCCAGCGTGATGCCGATACCGAACAATGGGTGATGAATGACCGCATCGAGCGCGCCCTGCCAATCGAGCATCATGGCTGCTCCTCa    >  1:1155497/1‑150 (MQ=255)
                                                                                                                                         aaTATAGGTTGCAGGAAGATCCAGCGGGTTTTCTCGTAGGCTGCCAGCACTACCTGATAGGCCCCCAGCGTGATGCCGATAc                                                                       >  1:2063470/1‑82 (MQ=255)
                                                                                                                                         aaTATAGGTTGCAGGAAGATCCAGCGGGTTTTCTCGTAGGCTGCCAGCACTACCTGATAGGCCCCCAGCGTGATGCCGATAc                                                                       <  2:2063470/82‑1 (MQ=255)
                                                                                                                                          aTACAGGTTGCAGGAAGATCCAGCGGGTTTTCTCGTAGGCTGCCAGCACTACCTGATAGGCCCCCAGCGTGATGCCGATACCGAACAATGGGTGATGAATGACCGCATCGAGCGCGCCCTGCCAATCGAGCATCATGGCTGCTCCTCACg  >  2:751374/1‑150 (MQ=255)
                                                                                                                                             |                                                                                                                                                  
GGTTAAGATAGAGCGGCACAGCCAGGGCCACCGTAGCGGGCCCCAGCAGGATGTTCATGATTTCGGTGCTTTTGCGGTACTCGCTGTAATCGATGCCGCACAGCAGCAACACACCGATCACCAGCAGCATCGATACCAATACAGGTTGCAGGAAGATCCAGCGGGTTTTCTCGTAGGCTGCCAGCACTACCTGATAGGCCCCCAGCGTGATGCCGATACCGAACAATGGGTGATGAATGACCGCATCGAGCGCGCCCTGCCAATCGAGCATCATGGCTGCTCCTCACG  >  NC_002947/5478290‑5478577

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